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Economies  of  scale  provide  the  main  incentive  for  regionalization 
of  v/astewater  treatment  facilities.   The  cost  of  collecting  wastewater 
increases  v/ith  the  size  of  plant  and  serves  to  offset  within  plant  scale 
economies.   Two  types  of  secondary  treatment  and  three  tertiary  treat- 
ment processes  were  considered  as  feasible  alternatives.   They  were 
high-rate  trickling  filter  system,  activated  sludge  system,  and  two 
clarifier  lime  clarification  processes,  line  recalcination,  and  ammonia 
stripping  process,  respectively. 

The  space  and  time  dimensions  of  regionalization  were  considered. 
In  order  to  find  the  potential  sites  of  regional  treatment  plants,  a 
cost  minimization  model  was  used.   For  the  consolidation  over  time  an 
exponential  objective  function  was  minimized  subject  to  discrete  ca- 
pacity requirements. 

Estimates  of  construction  and  operation  and  maintenance  costs  were 
based  upon  Environmental  Protection  Agency  (EPA)  studies.   The  rates  of 


cost  increases  over  time  vzere  also  estimated  from  EPA  cost  indexes  as 
.12,  .10,  and  .07  per  year  for  force  main  construction  cost,  pumping 
station  and  sewage  treatment  plant  construction  costs,  and  operation 
and  maintenance  costs  for  both  force  mains  and  sewage  treatment  plants, 
respectively.   Service  lives  for  pumping  station  and  wastewater  treat- 
ment were  taken  as  25  years,  and  for  pipelines  as  50  years.   An  inter- 
est rate  of  5%  per  year  was  assumed. 

A  least  cost  staging  policy  to  satisfy  the  wastewater  treatment 
requirements  of  the  nine  cities  in  Alachua  County,  Florida  for  the 
years  1975,  1980,  and  1990  was  determined  with  the  model.   The  results 
were  sensitive  to  the  required  quality  of  secondary  effluents. 

If  the  required  quality  of  secondary  effluent  is  set  at  above  80% 
removal,  activated  sludge  plants  V7ill  be  required  to  provide  the  least 
cost  method  of  treating  wastewater.   The  results  show  that  secondary 
treatment  cost  will  be  minimized  if  High  Springs,  Alachua,  Archer,  and 
Newberry  cooperate  to  build  a  regional  treatment  plant  in  1975  to 
satisfy  their  1980  dem.ands,  and  then  construct  individual  plants  in 
1990  to  satisfy  the  additional  demands  at  that  time.   The  other  cities 
would  satisfy  their  demands  by  building  individual  activated  sludge 
plants  in  1975  to  satisfy  their  1980  demands  with  additions  being  made 
in  1990.   The  minimum  tertiary  treatment  cost  combination  would  require 
the  construction  of  two  regional  treatment  plants,  one  for  High  Springs, 
Alachua,  Archer,  and  Newberry,  and  one  for  Waldo  and  Hawthorne,  in  1975 
to  satisfy  the  1980  demands  of  these  six  cities.   Their  additional  de- 
mands in  1990  would  be  met  with  plants  added  in  each  city.   Gainesville. 
LaCrosse,  and  Micanopy  would  satisfy  their  tertiary  treatment  dem.ands 
with  individual  plants  constructed  in  1975  with  additions  in  1990. 


If  the  required  quality  of  secondary  effluent  is  set  at  80a  re- 
noval,  the  results  shovj  that  secondary  treatment  cost  for  the  county 
will  be  minimized  if  each  city  bui]ds  a  trickling  filter  pl.ant  which 
satisfies  1980 's  wastc-ater  treatment  demands  in  1975  with  additions  to 
capacity  being  made  in  1990.   The  minimum  tertiary  treatment  cost  com- 
bination suggests  that  there  v/ill  be  regionalization  of  treatment  in 
1975  to  satisfy  the  19S0's  demands  of  three  citie3--Newberry,  Archer, 
and  Alachua.   Additional  demands  in  1990  were  satisfied  by  building  in- 
dividual treatment  plants  in  these  three  cities.   The  other  six  cities 
would  satisfy  their  tertiary  treatment  with  individual  plants  con- 
structed on  the  same  schedule  as  the  trickling  filter  plants. 

The  possibility   of  substituting  land  spreading  for  tertiary 
treatment  Is  discussed,  and  the  factors  which  should  be  considered  in 
land  treatment  are  analyzed  with  the  aid  of  an  optimization  model.   The 
break-even  distances  of  transmission  pipelines  for  land  treatment  as  a 
substitute  for  tertiary  treatment  are  calculated  and  presented. 


CHAPTER  I 
INTRODUCTION 


The  Problem  and  Objectives 

The  environmental  thrust  of  the  1970 's  has  resulted  in  regulations 
to  clean  up  wastewater  from  a  number  of  sources  including  industry, 
municipalities,  and  nuclear  power  generating  stations.   Communities 
throughout  the  country  faced  with  meeting  the  1985  goal  of  zero  dis- 
charge of  pollutants  into  streams  encounter  both  economic  and  environ- 
mental constraints  [Public  Law  92-500,  1972].   The  high  costs  of 
providing  primary,  secondary,  and  especially  tertiary  treatment  by  con- 
ventional methods  place  an  economic  burden  on  communities,  especially 
on  the  smaller  ones  typically  found  in  rural  areas. 

Several  wastewater  management  alternatives  exist,  including 
(1)  advanced  biological  treatment,  (2)  physical-chemical  treatment,  and 
(3)  land  treatment.   Comjnunities  upgrading  their  treatment  systems  must 
determine  which  individual  or  combination  of  these  treatment  methods 
has  the  greatest  potential  for  meeting  both  economic  and  environmental 
constraints . 

The  concept  of  land  treatment  is  receiving  increased  attention  as 
an  alternative  wastewater  management  approach  [Dalton  and  Murphy,  1973; 
Egeland,  1973;  Parizek  et_  aJL.  ,  1967;  Sopper  and  Kardos,  1973;  Thomas, 
1973;  Williams  et  al. ,  1969].   Briefly  stated,  land  treatment  involves 
the  use  of  agricultural  land  and  crops  or  forest  products  to  absorb  and 
filter  nitrates,  phosphates,  and  other  elem.ents  from  wastewater  that 


has  undergone  primary  and,  usually,  secondary  treatment.   Water  filtered 
by  the  soil  is  then  returned  to  the  underground  water  supply. 

The  many  issues  surrounding  the  use  of  land  for  treatment  of  muni- 
cipal wastewater,  and  subsequent  land  reclamation  through  soil  infil- 
tration and  plant  growth,  encourages  a  m.ultidisciplinary  planning 
approach.   Systematic  investigation  of  institutional,  economic,  and 
social  issues  should  be  made  concurrently  with  technical  and  engineering 
studies.   To  date,  most  attention  has  focused  on  technical  and  engineer- 
ing questions.   An  objective  of  this  study  is  to  articulate  some  of 
these  issues  and  to  suggest  a  framework  which  social  scientists,  and 
agricultural  economists  in  particular,  might  use  as  a  starting  point 
for  further  inquiry  into  the  socioeconomic  aspects  of  land  treatment. 

In  the  planning  and  design  of  a  wastewater  treatment  facility  for 
a  city,  town,  or  an  area  vrhere  the  requirements  for  treatm.ent  are  ex- 
pected to  increase  with  time,  the  initial  size  of  the  treatment  plant 
and  the  timing  of  capacity  additions  and/or  replacements  over  some  time 
horizon  need  to  be  answered  in  the  context  of  an  optimal  staging  policy. 
Such  a  policy  is  affected  by  the  v;astewater  treatment  requirem.ents ,  the 
rates  of  interest  and  inflation,  construction  cost,  operating  costs, 
maintenance  and  repair,  service  life,  and  the  staging  efficiency  of  the 
system  to  be  designed. 

Recently,  the  concept  of  regionalization  of  wastewater  treatment 
has  been  suggested  as  an  effective  means  of  meeting  the  water  quality 
goals  at  a  minimum  cost.   Although  economy  of  scale  provides  the  pri- 
mary incentive  for  regionalization,  a  host  of  other  advantages  such  as 
more  qualified  operating  personnel   and  higher  degrees  of  automation 
may  be  gained  through  the  utilization  of  this  concept. 


The  recent  UteraU,re  of  citv  planning  and  regional  science  con- 
tains few  analytical  studies  concerned  with  facility  planning  for 
pnhlic  urhan  systems.   However,  the  .anage.ent  of  these  systems  needs  a 
coherent  framework  for  planning  which  deals  with  facility  si.a.  loca- 

tion,  and  tiaiing. 

This  study  develops  a  decision  .odel  to  use  as  a  guide  in  planning 
facilities  for  a  specific  wastewater  treatment  syste..   An  optimisation 
technique  will  be  ntill.ed  in  selecting  the  opti.u»  program.   The  final 
plan  will  include  the  nunber  of  treatment  plants  needed,  treatment 
plant  capacities,  waste  sources  to  be  served  by  each  plant,  disposal 
methods,  and  plant  locations. 

Meth£d_of_Study   |/ 

The  selection  of  an  optimal  regional  plan  over  time  is  a  complex 
problem,  the  solution  of  which  requires  a  considerable  variety  and  a 
quantity  of  information,  a  great  expenditure  of  time,  money,  and  a  wide 
diversity  of  expertise.   The  optimum  plan  provides  the  least  cost 
system  of  collection,  treatment,  and  disposal  of  wastewater  subject  to 
the  physical,  financial,  and  sociopolitical  constraints  of  the  region. 
Choosing  an  optimal  plan  from  among  a  large  number  of  alternatives 
(which  involve  different  construction  costs,  plant  location,  facilities 

o>  -Jc  n.n-1-p  difficult  without  a  formal 
and  capacity,  and  farming  programs)  is  quite  dxtricux 

,  ^   ,  i^-r'.tP  nnmber  of  variables  that  are  deter- 
model  that  can  accomodate  a  l^rge  number  oi 

Mn.d  simultaneously  in  a  dynamic  frameworlc.   The  steps  to  develop  such 
a  model  and  the  lin.s  between  them  are  shown  in  Figure  1.1.   The  steps 
consist  of  two  sets:   (D  the  estimation  of  the  form  of  cost  functions 
and  the  values  of  uncontrolled  variables;  and  (2)  the  choice  of  values 
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of  decision  variables.   Decisions  deal  v/ith  sizes,  types,  locations, 
and  timing  of  plants  and  facilities  to  be  constructed. 

In  this  study,  a  mathematical  optimization  model  is  proposed  \7hich 
will  minimize  the  present  value  of  total  expenses  over  the  planning 
horizon  subject  to  the  requirements  for  wastewater  treatment  and  land 
availability  in  the  region  of  interest.   Treatment  plant  facility  costs 
as  functions  of  capacity,  land  requirements,  wastewater  production  and 
land  values  as  functions  of  time  and  population,  and  information  on 
zoning  and  population  growth  estimates  were  incorporated  in  the  model 
to  determine  a  minimum  cost  expansion  program. 

Area  of  Study  and  Source  of  Data 

Alachua  County,  Florida  w^as  chosen  as  the  study  area.   The  county 
has  a  total  land  area  of  568,320  acres  or  888  square  miles.   These 
figures  do  not  include  the  large  bodies  of  water  representing  about 
76.4  additional  square  miles  [North  Central  Florida  Regional  Planning 
Council,  1973].   Figure  1.2  shov;s  the  location  of  Alachua  County. 
There  are  nine  incorporated  cities  or  population  centers  in  Alachua 
County.   Only  one  of  these  nine  incorporated  communities  has  adequate 
wastewater  collection,  treatment,  and  disposal;  the  remaining  eight 
incorporated  municipalities  have  either  privies  or  septic  tanks  w^ith 
the  large  flow  users  such  as  hotels,  restaurants,  and  service  stations, 
primarily  located  along  major  thoroughfares,  being  served  by  small 
package  plants. 

Population  estimates  for  this  study  were  obtained  from  the  esti- 
mates provided  by  North  Central  Florida  Regional  Planning  Council 
[1973].   The  relationships  between  capacity  and  costs,  both  construction 
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Figure  1.2   Location  of  Alachua  County 


and  operation  and  maintainence,  were  obtained  from  EPA  studies  [Michel 
and  Johnson,  1970;  Smith  and  McMichael,  1969;  Smith  and  Eilers,  1971]. 
The  information  about  cost  increases  over  time  was  provided  by  EPA  cost 
indexes.  Other  knowledge  required  to  complete  this  study  was  obtained 
from  journals,  textbooks,  and  miscellaneous  publications  on  wastewater 
treatment. 

This  study  developed  a  procedure  to  find  the  optimum  number,  sizes, 
and  locations  of  wastewater  treatment  facilities  in  Alachua  County. 
The  time  period  involved  is  from  1975  to  1990,  and  wastewater  sources 
are  only  the  nine  incorporated  cities  in  the  county.   The  concept  of 
regionalization  used  in  this  study  was  that  more  than  one  city  might 
transfer  its  waste  to  a  regional  treatment  facility  in  some  selected 
location. 


CHAPTER  II 
RELEVANT  ECONOMIC  THEORY 


This  chapter  develops  a  theoretical  framework  for  structuring 
the  mathematical  optimization  m.odel.   Cost  minimization  theory',  loca- 
tion theory,  and  multidimensional  utility  analysis  are  combined  to  de- 
termine the  optimum  number,  size,  and  location  of  v;astewater  treatment 
plants. 

Cost  I? jnim-ization  with  Respect  to  Output 

Corsider  a  system  of  equations  consisting  of  a  production  function 

(2.1),  a  cost  equation  (2.2),  and  an  expansion  path  function  (2.3): 

q   =   f(x^,  X2)  (2.1) 

c   =   r   X   +  r   X  +  b  (2.2) 

o   =   g(:-:^,  x^)  (2.3) 

where  q  is  output,  x^ ,  x„  are  variable  inputs,  r^  and  r„  are  the  re- 
spective prices  of  x   and  x  ,  and  b  is  the  cost  of  the  fixed  inputs. 
This  system  of  three  equations  in  four  variables  can  be  reduced  to  a 
single  equation  in  which  cost  is  stated  as  an  explicit  function  of  the 
level  of  output  plus  the  cost  of  the  fixed  inputs: 

c   =   ^(q)  +  b  (2.4) 

The  fixed  cost  must  be  paid  regardless  of  how  much  the  firm  produces, 
or  whether  it  produces  at  all.   The  cost  function  gives  the  minimum 

This  system  of  equations  can  be  expanded  to  include  more  than  two 
variable  inputs.   [See  Ferguson,  1971,  pp.  154-168] 


cost  of  producing  each  output  and  is  derived  on  the  assumption  that  the 
entrepreneur  acts  rationally. 

Let  the  levels  of  the  entrepreneur's  fixed  inputs  be  represented 
by  a  parameter  k,  which  gives  the  "size  of  his  plant"' — the  greater  the 
value  of  k,  the  greater  the  size  of  his  plant.   The  entrepreneur's 
short-run  problems  concern  the  optimal  utilization  of  a  plant  of  given 
size.   In  the  long-run  he  is  free  to  vary  k  to  select  a  plant  of  optimum 
size.   The  shapes  of  the  entrepreneur's  production  and  cost  functions 
are  given  in  the  short-run  and  depend  upon  fixed  plant  size.   In  the 
long-run  he  can  choose  a  production  function.   Once  he  has  selected  this 
function,  i.e.,  selected  a  value  for  k,  he  is  faced  v/ith  the  conven- 
tional short-run  optimization  problems. 

Assume  that  k  is  continuously  variable  and  introduce  it  explicitly 
into  the  production  function,  cost  equation,  and  expansion  path 
function: 

q   =   f(x^,  X2,  k)  (2.1a) 

c   =  r^  x^  +  r^  x^  +  if'(k)  (2.2a) 

o  =  g(x^,  x^,  k)  (2.3a) 

Fixed  cost  is  an  increasing  function  of  plant  size:  —7; >  0.   As 

dk 

before,  total  cost  may  be  expressed  as  a  function  of  output  level  and 
plant  size: 

c  =   <{>(q,  k)  +  ^^(k)  (2.4a) 

The  entrepreneur's  long-run  total  cost  function  gives  the  minimum. 
cost  of  producing  each  level  of  output  as  a  function  of  plant  size. 
This  function  is  the  envelope  of  the  short-run  functions:   it  touches 
each  and  intersects  none.   Write  the  equation  for  the  family  of  short- 
run  cost  functions  (2.4a)  in  implicit  form: 
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C  -  i<(q,    k)  -  i|^(k)   -   G(C,  q,  k)   =   0  (2.4h) 

and  set  the  partial  derivative  of  (2.4b)  v/ith  respect  to  k  equal  to 
zero: 

G,  (C,  q,  k)   =   0  (2.4c) 

The  equation  of  the  envelope  curve  (the  long-run  cost  curve)  is  obtained 
by  eliminating  k  from  (2.4b)  and  (2.4c)  by  solving  for  C  as  a  function 
of  q,  i.e.: 

C  =   (j>(q)  (2.4d) 

Long-run  total  cost  is  a  function  of  output  level,  given  the  condition 
that  each  output  level  is  produced  in  a  plant  of  optimum  size. 

Since  fverage  cost  equals  total  cost  divided  by  output  level,  the 
minimum  average  cost  of  producing  a  particular  output  level  is  attained 
at  the  same  plant  size  as  the  minimum  total  cost  of  producing  that  out- 
put level.   The  long-run  average  cost  curve  can  be  derived  by  dividing 
long-run  total  cost  by  output  level,  or  by  constructing  the  envelope  of 
the  short-run  average  cost  curves.   The  two  constructions  are  equivalent. 

Generally,  long-run  average  cost  (LAC)  is  considered  as  a  planning 
device,  which  given  the  state  of  the  arts  would  tell  an  entrepreneur 
the  plant  size  which  will  enable  him  to  produce  a  desired  output  level 
at  the  least  possible  cost  per  unit.   The  significance  of  plant  size  in 
relation  to  the  plant  costs  results  from  the  nature  of  the  economies  of 
scale. 

The  possible  existence  of  substantial  scale  economies  v/hich  would 
permit  relatively  large  treatment  plants  to  process  their  wastewater  at 
lower  average  cost  per  unit  than  relatively  smaller  treatment  plants  is 
a  key  factor  in  stimulating  interest  in  regional  wastewater  treatment. 
The  principal  basis  of  scale  economies  is  specialization,  or  the 
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division  of  labor — a  phenomenon  Adam  Smith  deemed  so  central  that  he  de- 
voted the  first  three  chapters  of  Wealth  of  Nations  to  it.   In  Smith's 
view,  great  increases  in  the  productivity  of  labor  (and  hence  great  re- 
duction in  product  cost)  v/ere  due  to  three  repercussions  of  the  division 
of  labor:   an  increase  in  worker  dexterity,  the  saving  of  time  coranionly 
lost  in  "passing  from  one  species  of  work  to  another,"  and  "the  inves- 
tion  of  a  great  num.ber  of  machines  which  facilitate  and  abridge  labor, 
and  enable  one  man  to  do  the  work  of  many." 

A  somewhat  different  basis  of  scale  economies  is  found  in  the  pro- 
cess industries  such  as  petroleum  refining,  chemica]  production,  cement 
making,  glass  manufacturing,  and  steam  generation.   The  output  of  a  pro- 
cessing unit  tends  within  certain  physical  limits  to  be  roughly  propor- 
tional to  the  volume  of  the  unit,  other  things  being  equal;  while  the 
amount  of  materials  and  fabrication  effort  (and  hence  investment  cost) 
required  to  construct  the  unit  is  more  apt  to  be  proportional  to  the 
surface  area  of  the  units  reaction  chambers,  storage  tanks,  connecting 
pipes,  etc.   Since  the  area  of  a  sphere  or  cylinder  of  constant  propor- 
tions varies  as  the  two-thirds  power  of  volume,  the  cost  of  constructing 
process  industry  plants  can  be  expected  to  rise  as  the  two-thirds  power 
of  their  output  capacity,  at  least  up  to  the  point  where  they  become  so 
large  that  extra  structural  reinforcement  and  special  fabrication  tech- 
niques are  required.   There  is  considerable  empirical  support  for  the 
existence  of  this  "two-thirds  rule,"  which  is  used  by  engineers  in  esti- 
mating the  cost  of  new  process  equipment  [Haldi  and  Whitecomb,  1967; 
Koorc,  1959;  Schuman  and  Alpert,  I960]. 

Still  another  benefit  of  size  arises  from  what  E.  A.  G.  Robinson 
calls  "the  economies  of  massed  reserves"  [Robinson,  1958,  pp.  26-27]. 
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A  firm  anxious  to  maintain  continuity  of  production  must  hold  equipir.ent 
in  reserve  against  machine  breakdowns.   A  firm  large  enough  to  use  only 
one  specialized  machine  may  be  forced  to  double  its  capacity  if  it 
insists  on  hedging  against  breakdown;  the  larger  firm  with  numerous  ma- 
chines can  obtain  virtually  the  same  degree  of  protection  by  holding 
only  a  small  proportion  of  its  capacity  in  reserve.   Likewise,  the  num- 
ber of  repairmen  a  company  must  employ  to  provide  any  stipulated  amount 
of  service  in  the  event  of  random  machine  failure  rises  less  than  pro- 
portionately with  the  number  of  machines  in  operation  [VJhitin  and 
Peston,  1954]. 

Two  other  common  arguments  are  the  theory  of  "perfect  divisibility" 
promoted  by  Kaldor  [l943]  and  Lerner  [1946,  pp.  186-199],  and  the 
theory  of  "proportionality"  developed  by  Chamberlin  [l948].   The  theory 
of  perfect  divisibility  states  that  the  inperfect  divisibility  of 
factors  explains  economies  of  scale  and  that  v;ith  perfect  divisibility 
the  economies  of  scale  would  be  absent.   In  criticizing  the  theory  of 
perfect  divisibility,  Chamberlin  argued  that  it  neglects  the  effects 
of  divisibility  of  the  efficiency  of  factors.   Rather,  he  claimed  that 
all  economies  of  scale  villi   be  explained  by  the  proportion  of  factors. 
The  theory  of  proportionality  indicates  that  there  is  a  certin  optimum 
proportion  of  factors,  and  because  factors  are  obtainable  only  in  dis- 
crete units,  this  optimum  proportion  can  be  closely  approximated  only 
when  the  aggregate  of  factors  is  large.   In  other  words,  small  scale 
production  will  be  relatively  inefficient  in  achieving  these  optimum 
proportions. 

It  is  quite  clear  that  econom.ies  of  scale  do  exist,  and  that  unit 
costs  decline  with  increases  in  plant  and  firm  size,  at  least 
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within  limits.   However,  there  are  several  reasons  for  believing  that 
scale  economies  are  limited.   First,  in  nearly  all  production  and  dis- 
tribution operations  the  realization  of  scale  economies  appears  to  be 
subject  to  diminishing  returns.   Sooner  or  later  a  point  is  reached  at 
which  all  opportunities  for  making  further  cost  reductions  through  in- 
creased size  are  exhausted.   Second,  it  is  possible  that  rising  unit 
costs  related  to  the  difficulty  of  managing  an  enterprise  of  increas- 
ingly large  scale  will  offset  and  eventually  over^rhelm  the  savings  at- 
tributable to  high-volume  production  and  distributions  [Robinson,  1958, 
Chapters  3,  10,  12;  Coase,  1937;  Williamson,  1967].   A  third  major  in- 
fluence which  may  prevent  economies  of  scale  from  being  realized  in- 
definitely is  the  cost  of  transportation. 

The  theory  of  the  firm  implies  that,  under  pure  competition, 
individual  firms  tend  to  operate  at  equilibrium  in  the  long  run,  seek- 
ing to  minimize  the  long-run  cost  with  respect  to  output  [Leftwich, 
1956,  pp.  167-177].   The  number  and  size  of  plants  of  a  given  firm  are 
assumed  to  be  determined  within  the  equilibrium  conditions.   However, 
wastewater  treatment  in  Alachua  County  is  a  regulated  utility  subject 
to  the  control  of  a  single  authority.   It  is  a  local  monopoly  and  its 
size  is  determined  largely  by  the  extend  of  local  demand.   The  distinc- 
tion between  perfect  competition  industries  on  the  one  hand  and  regu- 
lated and  imperfectly  competitive  industries  on  the  other  is  that,  in 
the  former,  the  size  of  the  market  plays  no  part  in  determining  the 
size  of  firm  but  merely  determ.ines  the  number  that  can  survive,  while 
in  the  public  utility  type  of  industry  the  spatial  distribution  of  de- 
mand determines  both  the  number  and  relative  scale  of  the  firms  in  the 
optimal  (cost-minimizing)  setup. 
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The  Theory  of  Location 

The  basic  concept  of  location  theory  is  to  incorporate  assembly 
costs  into  production  costs  as  the  basis  for  determining  the  cost  mini- 
mizing conditions.   This  combination  overco:?.es  the  shortcomings  of 
neglecting  the  selection  of  location  as  a  production  activity  in  the 
pure  theory  of  the  firm..   Since  the  production  origins  and  plants  are 
generally  not  at  the  same  location,  the  assembly  costs  are  incurred 
v-jhile  transferring  the  raw  material  from  production  origins  to  plants. 
Furthemore,  the  marketing  costs  incurred  while  transporting  final  pro- 
ducts from  plant  site  to  market  should  also  be  included  in  the  analysis. 
Assembly  costs  consist  of  two  parts  in  this  study.   The  collection 
costs,  v;hich  relate  to  the  process  of  removing  the  sewage  from  the  site 
of  production  by  small  sewer;  and  the  transmission  costs,  v/hich  relate 
to  transporting  the  sewage  from  the  small  sewers  by  larger  sevjers  to 
the  treatment  site.   Those  costs  are  directly  related  to  the  distance 
from  the  site  of  production  to  treatment  site,  the  amount  of  flow  trans- 
mitted, and  the  difference  in  elevation  between  the  tv7o  sites  to  be 
joined  by  sewers. 

The  pure  theory  of  the  firm  m.ust  be  combined  with  the  theory  of 
location  in  order  to  interpret  the  joint  problem  with  which  this  study 
is  concerned.   The  incorporation  of  assembly  activities  into  the  pro- 
duction process  makes  the  goal  of  the  firm  one  of  minimizing  the  total 
combined  cost  of  assembly  and  production  rather  than  those  associated 
with  production  alone. 

\iJhile  assembly  is  clearly  external  to  the  operations  of  a  particu- 
lar plant,  both  plant  operations  and  transportation  are  integral  parts 
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of  the  joint  function.   In  order  to  minimize  the  production  cost,  the 
firm  seeks  to  gain  the  economies  of  scale  by  increasing  plant  size. 
However,  as  the  scale  of  plant  increases,  there  will  be  accompanying 

changes  in  the  magnitude  of  the  required  assembly  costs  of  raw  materi- 

2 

als  as  well  as  the  required  costs  of  disposing  of  its  final  products. 

More  specifically,  increases  in  plant  size  will  be  associated  with  ex- 
panded supply  requirerp.ents,  and  hence,  with  increased  assembly  costs 
and  disposal  costs.   Therefore,  potential  economies  of  scale  will 
gradually  be  offset  by  increases  in  these  costs. 

For  instance,  consider  two  communities,  A  and  B,  which  both  have  a 
need  for  v^astewater  treatment  by  means  of  conventional  activated  sludge 
and  associated  processes.   If  neither  community  has  an  existing  vraste- 
water  treatment  plant  and  each  community  needs  a  one  mgd  (million  gal- 
lon per  day)  plant,  the  two  alternatives  are  for  each  community  to 
build  and  operate  its  own  plant  or  for  some  joint  effort  with  one  plant. 
Assume  the  costs  for  a  single  one  mgd  plant  is  25.5  cents/kgal,  22 
cents/kgal  for  a  two  mgd  plant  and  pipeline  cost  of  1.71  cents/kgal/ 

mile.   Under  these  assumptions,  if  two  communities  are  more  than  4.1 

3 

m.iles  apart,  each  community  would  build  a  one  mgd  plant.    The  same 

kind  of  argument  can  be  expanded  to  include  more  than  t\-io   communities 
with  different  wastewater  treatment  requirements. 


Disposal  costs  for  an  effluent  plant  are  somewhat  analogous  to 
marketing  costs  for  a  firm  producing  consumer  goods. 

3 

The  breakeven  length  pipeline  (L)  is  calculated  by  equating  the 
total  daily  cost  of  the  two  alternatives  and  solving  for  the  length  of 
the  pipeline  at  the  breakeven  point,  i.e.,  to  solve  (25.5)  x  1000  + 
25.5  X  1000  =  1.71  L  X  1000  +  22.0  x  2000  for  L,  which  gives  L  =  4.1 
miles. 
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However,  the  potential  location  for  a  plant  need  not  be  situated 
in  a  site  of  demand  or  supply.   The  question  is  how  location  decisions 
are  made.   In  the  ideal  case  they  are  made  with  perfect  knowledge  of 
costs  and  benefits  at  alternative  locations,  and  they  are  made  so  as  to 
maximize  net  benefits  from  operations  at  the  optimal  location.   Often 
in  the  real  v;orld  case,  perfect  knowledge  is  replaced  by  error  and  an 
adverse  reaction  to  risk,  and  the  search  for  optimal  locations  is  re- 
duced to  a  search  for  a  satisfactory  location. 

Transport  systems  v^ere  built  to  promote  interactions  over  space, 
and  prevailing  location  theory  models  are  largely  transportation 
oriented.   To  study  the  effect  of  cost  of  geographic  space  on  the  loca- 
tion of  human  activity  let  us  examine  the  location  problem  in  general 
terms  rather  than  in  terms  of  specific  industry  case  studies.   We  shall 
consider  only  manufacturing  activity  and  will  assume  that  wage,  inter- 
est, and  other  prices  are  equal  ever>n-7here .   Only  transportation  charges 
which  vary  with  distance  are  considered  in  the  analytical  approaches 
presented  below. 

Problems  on  location  analysis  can  be  classified  into  two  major 
structural  categories. 

A.   Location  on  a  plane,  which  is  characterized  by 

1.  an  infinite  solution  space.   That  is,  central  facilities 
may  be  located  anywhere  on  the  plane  and  are  confined  neither  to  nodes 
of  the  network  nor  to  points  on  the  links  between  these  nodes. 

2.  distance  measurement  according  to  a  particular  metric. 
One  example  is  the  Euclidean  metric  where 

d..   =   ((x.  -  x.)^  +  (y.  -  y.)^)''^ 

d..   =   the  distance  between  points  i  and  j. 
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X.,  y.   =   the  coordinates  in  a  rectangular  system  of  the 
^'   ^     ith  point. 

Another  example  is  the  metropolitan  metric  v;here 

d..   =   ixi-x.l  +  iy.  -  y.|. 

B.   Location  on  a  network  which  is  characterized  by 

1.  a  solution  space  consisting  of  the  points  on  the  netxrork 
(both  nodes  and  points  on  the  arcs  which  join  the  nodes), 

2.  distance  measurement  or  time  measurement  along  the  network, 

d..   =   the  length  (time)  of  the  shortest  path  [Cooper, 
"'"-'     1967]  from  node  i  to  node  j. 

Alfred  Weber  [1909]  pioneered  the  Location  Theory.   He  considered 
the  location  of  an  industry  between  two  resources  and  a  single  market 
where  the  criterion  was  minimization  of  transportation  cost.   However, 
it  is  only  two  decades  since  the  advent  of  mathematical  programming  and 
computer-aided  computation  that  the  problem  has  received  real  attention 
from  the  standpoint  of  research. 

Werson  et  al.  [1962]  presented  one  of  the  earliest  modern  consid- 
erations of  location  problems.   They  located  solid  waste  disposal  sites 
to  minimize  hauling  costs  in  a  metropolitan  environment.   Travel  was 
assumed  to  be  along  the  rectangular  grid  which  typifies  many  American 
cities.   An  optimal  single  disposal  site  was  determined  by  linear  pro- 
gramming to  be  at  the  median  of  the  generating  sites. 

Beginning  with  the  formulation  of  Cooper  [1963]  and  Kuhn  and 
Kuenne  [1962],  interest  in  location  analysis  has  quickened.   Their 
work  which  appeared  independently  described  an  iterative  process  for 
solving  the  generalized  Weber  problem.   The  problem  is  to  find  the 
single  point  which  minimizes  the  sum  of  the  v/eighted  Euclidean  distances 
to  that  point.   The  objective  is 
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minimize  z   =   ^   v; .  d 
i=3   ^   ^ 
where : 

V7.   =   the  weight  attached  to  the  i   point  (goods  de- 
manded, resources  sent,  population,  etc.); 

2  2  ^ 

d.   =   ((x.  -  x)   +  (y.  -  y)  )',  which  is  the  Euclidean 

distance  from  point  i  to  central  point; 

X.,  y.   =   the  location  of  the  i   point  relative  to  some 
i    i 

fixed  cartesian  coordinate  system; 

X,  y   =   the  unknown  coordinates  of  the  central  point; 

n   =   the  number  of  points  served. 

Partial  differentiation  with  respect  to  x  and  y  yields  a  pair  of 

equations  that  are  the  first  order  conditions  for  a  minimum: 

9z 


3x     -     -^  ^  (2.5) 

8z        w.  (y.  -  y) 

-97  =  h-^—f. =  °  (2.6) 

1 

Somev/hac  more  attention  has  been  devoted  to  problems  of  locating  central 

points  on  a  network.   The  problem  of  warehouse  or  plant  location  has  the 

following  general  characteristics.   Given  a  number  of  demand  areas  of  a 

certain  product,  each  with  a  demand  D.,  and  a  number  of  alternative 

X 

sites  where  facilities  may  be  built  to  satisfy  these  demands,  determine 
where  the  facilities  should  be  placed  and  which  demand  areas  are  to  be 
served  by  a  given  facility.   The  objective  is  that  the  sum  of  the  trans- 
portation cost  and  the  amortized  facility  cost  is  minimized. 

A  tradeoff  exists  betv/een  facility  and  transport  costs;  clearly, 
the  greater  the  number  of  well-placed  facilities,  the  lower  will  be  the 
cost  o£  distribution.   But  as  shipment  cost  decreases,  the  investment 
in  facilities  must  rise.   At  some  number  of  facilities  the  total  cost 
should  be  a  minimum.   Beyond  that  point,  the  cost  of  adding  a  facility 
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exceeds  the  savings  in  distribution  cost.   It  is  assun:ied  that  both  the 
number  of  demand  areas  for  products  and  the  number  of  plant  or  ware- 
house sites  are  finite. 

Maranzana  [1954]  considered  warehouse  location  \-7here  distances  are 
not  Euclidean  but  are  measured  on  the  road  network.   He  also  utilized 
an  heuristic  procedure  to  locate  a  specified  number  of  warehouses  to 
serve  a  region  of  known  demands.   His  criterion  was  the  minimization  of 
transport  costs.   Kuehn  and  Hamburger  [1963]  also  presented  a  heuristic 
method  for  warehouse  location;  their  approach  was  to  minimize  the  sum 
of  transport  and  v/arehousing  costs. 

ReVelle  and  Sv7ain  [1970]  consider  a  location  problem;  their  objec- 
tive was  the  m.inimization  of  average  time  or  distance  vjhich  people  must 
travel  to  the  facilities.   The  problem  is  constrained  by  a  fixed  number 
of  facilities.   The  model  was  structured  as  a  0-1  mix  integer  program- 
ming problem  with  facilities  restricted  to  nodes  of  the  network. 

In  this  study,  equations  (2.5)  and  (2.6)  were  used  to  find  the 
potential  sites  for  treatment  plants.   Euclidean  metric  distance  meas- 
urement was  used  to  determine  the  optimal  locations  for  wastewater 
treatment  facilities.   The  solution  space  consists  of  potential  sites 
and  sources  of  wastewater  treatment  demands. 

Multidiiiientional  Utility  Analysis 

Since  their  inception,  cost  minimizing  models  for  firni  behavior 
have  certainly  not  been  v;ithout  their  critics.   But  for  the  most  part 
the  critics  have  failed  to  propose  substantive  alternatives  for  theo- 
retical use.   An  approach  not  directly  involving  simply  cost  minimiza- 
tion was  suggested  by  Scitovsky  [1943].   In  particular,  he  showed  that 
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if  an  entrepi^eneur  attempts  to  maximize  satisfaction,  he  will  expend 
the  same  amount  of  effort  as  if  he  v.'ere  attempting  to  maximize  profit 
only  in  the  special  case  in  which  the  marginal  rate  of  substitution 
between  entrepreneruial  activity  and  money  income  is  ind-';pendent  of 
the  level  of  money  income.   The  point  of  interest  is  the  applicability 
of  multidimensional  utility  theory  to  a  certain  class  of  raicroeconomic 
problems. 

The  theory  of  multidimensional  vector  ordering,  or  what  is  now 
more  generally  called  lexicographic  ordering,  has  the  following  meaning. 
Consider  two  alternatives,  which  may  be  bundles  of  commodities,  combina- 
tions of  lottery  tickets,  business  objectives,  etc.,  i.e.,  x'^  =  (x^ , 

x2,  ••.,  y^^)    and  x^  =  (x| ,  x]- ,  ....  x^).   Let  u  be  a  preference  index 
/-        n  1    2        n 

function.   A  regular  ordering  ranks  u(x'^)  >  u(x^)  if,  and  only  if,  x9  2. 

x.  for  all  i  and  the  strict  inequality  holds  for  at  least  one  component. 

In  a  lexicographic  ordering  a  hierachy  of  wants  is  recognized;  the 

components  of  the  vector  x  are  not  regarded  as  equally  important.   For 

convenience,  let  the  elements  of  each  vector  be  numbered  so  that  x   is 

more  important  than  x  ,  x   is  more  important  than  x  ,  etc.   Then  uCx'^)  > 

u(xM  if  x?  >  x} ,  irrespective  of  the  relationships  between  x?  and  xt , 
11  11 

for  i  =  2,  3,  ...,  n.   If  x^  =  x  ,  comparison  is  based  upon  the  second 
component.   Thus,  u(x^)  >  u(x^)  if  x^  =  x   and  x^  >  x^ ,  etc.   Proceeding 
in  this  manner,  vector  elements  associated  with  variables  lower  in  the 
hierachy  of  wants  are  considered  only  after  the  higher  order  wants  are 
satisfied  [Encarnacion,  1964]. 

Let  us  consider  a  wastewater  treatment  system  in  a  basin,  in  which 
it  is  assumed  that  the  policy  makers  attempt  to  minimize  wastev/ater 
treatment  cost  for  the  entire  basin  subject  to  maintaining  a 
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satisfactory  level  of  vrater  quality.   This  model  actually  has  two  possi- 
ble outcomes.   First,  if  the  level  of  v/ater  quality  is  less  than  the 
satisfactory  level,  the  policy  makers  disregard  the  cost  minimization 
goal  and  behave  as  though  they  were  purely  radical  environmentalists. 
Second,  if  the  satisfactory  level  of  water  quality  is  achieved,  the 
policy  makers  would  seek  to  minimize  wastewater  treatment  costs  for  the 
entire  basin.   In  the  terminology  of  multidimensional  utility  analysis, 
level  of  water  quality  is  the  dominant  component  and  cost  minimizing  is 
the  subordinate  one. 

Denote  the  level  of  water  quality  by  x   and  costs  by  x  .   Thus  any 
situation  is  represented  by  the  vector  x  =  (x  ,  x  ) .   Suppose  the  satis- 
factory level  of  v/ater  quality  is  x  .   According  to  the  description  of 
the  model,  further  v.'ater  quality  does  not  enhance  the  utility  or  satis- 
faction of  the  policy  makers.   In  other  words,  if  u(x)  is  the  utility 
function, 

^^  ]      ^  =   0 
3x   X   >  X 

Thus  the  optimum  vector  x   is  found  by  selecting  the  policy  that  mini- 
mizes X   (costs)  subject  to  x  ^  >^-|  •   In  particular,  let  us  compare  tvro 
situations  x'^  =  (x^,  x^)  and  x^  =  (x| ,  x^).   The  former  is  preferred  to 

the  latter  if,  and  only  if,  x  >^   x   for  i  =  0,  1,  and  x°  <  xL   If  this 

i    ^ 

problem  is  not  feasible,  i.e.,  x   <  x   for  i  =  0,  1,  the  optimum  vector 

X   is  the  one  whose  first  component  is  greater  (x   >  x::  for  i  7^  j). 

This  simple  two-variable  model  can  easily  be  generalized.   Suppose 

there  are  M  goals  x  ,  x  ,  ...,  x   arranged  in  order  of  descending 
i    ^        m 

importance.   Further,  let  each  goal  be  defined  so  that 

•^7—  J       *  >  0  and  ^r ]       *   =   0 

3x.   x.  <  x.         9x.   X.  >  X. 
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Every  possible  situation  is  described  by  a  vector  x  =  (x  ,  x  ,  ...,  x  ); 

-L    2        n 

and  the  optimal  vector  x   is  found  by  selecting  the  set  of  policies 
which  solves  the  following  constrained  maximization  (or  minimization) 

A 

problem:   maxlTiize  x   subject  to  x.  >  x.  for  i  =  1,  ? m-1.   If 

m  1—1  >   J     5 

this  problem  is  unfeasible,  the  least  important  goal  (x  )  is  dropped 

from  consideration.   The  new  problem  accordingly  becomes:   maximize 

x  -,    subject  to  X.  >  X.  for  i  =  1,  2,  ,  m-2.   One  thus  works  in  se- 

m-1  1—1 

quence  until  a  feasible  problem  is  determined,  all  lower  ordered  goals 
being  discarded  in  route.  If  goal  k(<ra)  is  the  least  important  objec- 
tive for  which  a  feasible  problem  exists,  the  vector  x*^  is  the  optimal 

vector  if,  a'ld  only  if,  x.  >  x.  for  i  =  1,  2,  ...,  k-1  and  x?  >  yJ}    for 

1—1  k    k 

all  j  /  0. 

In  this  study,  treating  V7astewater  is  the  dominant  component  and 
cost  minimizing  is  the  subordinate  one.   Possible  wastewater  treatment 
methods  will  be  investigated,  and  their  costs  will  be  estimated;  a 
decision  will  be  made  according  to  the  framework  of  multidimensional 
utility  analysis. 


CHAPTER  III 

A  GENEEIAL  MODEL  FOR 

REGIONAL  WASTEWATER  TREATMENT 


Regionalization  may  consist  of  the  centralization  of  administra- 
tion (separate  sewerage  systems  operated  by  a  single  authority)  or  one 
central  plant  serving  several  municipalities.   There  are  tvjo  ways  to 
centralize  the  wasteti/ater  treatment  for  the  second  case.   One  sugges- 
tion which  has  been  made  to  control  pollution  along  reaches  of  a  system 
is  to  collect  the  sewage  from  all  the  cities  and  firms  along  the  stream 
and  treat  it  at  one  central  downstream  location.   Another  suggestion 
often  made  is  that  small  outlying  cities  or  villages  transport  their 
waste  to  a  regional  treatment  facility  in  a  large  city.   The  feasi- 
bility of  such  proposals  can  be  judged  by  comparing  the  cost  of  trans- 
porting waste  with  the  treatment  savings  which  can  be  realized. 

In  order  to  meet  the  condition  that  the  available  treatment  plant 
capacity  at  any  time  during  the  study  period  is  sufficient  to  meet  the 
requirement,  there  are  many  expansion  policies  that  can  be  chosen. 
Planning  an  expansion  scheme  involves  the  choices  of  choosing  plant 
capacity,  building  site,  and  treatment  facilities  (such  as  land  spread- 
ing versus  integrated  biological-chemical  systems  and  the  like).   Also, 
different  rates  of  interest  and  inflation  can  change  the  preferred 
ordering  of  alternative  schemes. 

In  this  chapter,  a  mathematical  model  is  developed  for  the  purpose 
of  choosing  an  expansion  scheme  v;hich  minimizes  the  present  value  of 


23 


24 


investment  given  construction  costs  and  wastev/ater  treatment  require- 
ments over  time. 

The  Model 

The  primary  factors  which  will  be  considered  in  selecting  a 
regional  plan  in  this  study  include: 

a)  waste  sources  to  be  served; 

b)  number  and  capacities  of  treatment  plants  needed; 

c)  treatment  plant  locations; 

d)  desired  degrees  of  treatment; 

e)  collection  system;  and 

f)  cost  functions. 

Arrangemient  of  outfall  structure,  cost  allocation  among  partici- 
pants, methods  of  financing,  and  methods  of  im.plementation  v;ill  not  be 
considered  in  the  model.   The  influence  of  factors  (a)  -  (f)  on  the 
selection  of  an  optimal  regional  wastevzater  management  system  is  de- 
scribed more  fully  below. 

The  development  of  a  sound  regional  wastewater  treatment  plan 
requires  a  full  knowledge  of  the  major  existing  and  potential  sources 
of  wastes.   The  type  of  the  treatment  facility,  physical,  chemical,  or 
biological,  as  well  as  the  design  of  the  collection  system,  pumping 
stations,  and  treatment  depends  on  the  quantity  and  characteristics  of 
the  wastes.   Therefore,  characterization  and  quantification  of  all  the 
major  sources,  both  municipal  and  industrial,  within  the  region  may  be 
considered  as  one  of  the  most  important  factors  in  the  development  of 
an  effective  wastev/ater  management  plan. 
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The  locations  of  v^astewater  treatment  plants  constitute  a  major 
factor  in  the  overall  plan  for  optimization  of  v;astewater  management 
in  a  region.   Such  factors  as  the  water  quality  criteria,  the  distribu- 
tions of  ^^?aste  sources,  cost  of  treatment  facility,  cost  of  pumping 
stations,  and  interceptor  sewers  must  be  considered. 

A  general  model  is  to  supply  answers  regarding  the  location  of 
treatment  facility,  capacities,  and  the  number  of  treatment  plants. 
Assume  that  we  have  I  wastewater  sources,  J  potential  regional  and/or 
individual  plant  locations,  K  types  of  treatment  plant,  and  a  planning 
horizon  from  time  1  to  T.   In  the  following  relations: 

i  refers  to  waste  source;  i   =  1,  ...,  I 

j  to  plant  location;  j   =   1,  ...,  J 

k  to  type  of  treatment  plant;       k   =   1,  ...,  K 
t  to  time  period;  t   =   1,  .  .  .  ,  T 

d   to  amount  of  wastewater  generated  at  source  i  at  time  t; 

CC    to  present  value  of  wastewater  treatment  costs  at  loca- 
tion j  for  type  k  facility  at  time  t; 

CS.    to  present  value  of  wastewater  transmission  cost  from 
ij  t         .     -     .    . 

source  i  to  location  j  at  time  t; 

RC    to  required  wastewater  treatment  capacity  at  location  j 
for  type  k  facility  at  time  t; 

.    to  required  wastewater  transmission  capacity  from 
source  i  to  location  j  at  time  t; 

C.,   to  existing  wastewater  treatment  capacity  at  location  j 
for  type  k  facility  at  time  t; 

S     to  existing  wastewater  transmission  capacity  from 
source  i  to  location  j  at  time  t; 

CX    to  new  wastewater  treatment  capacity  to  be  built  at 
location  j  of  type  k  at  time  t; 

SX    to  new  wastewater  transmission  capacity  to  be  built 
from  source  i  to  location  j  at  time  t;  and 

Z..   to  the  amount  of  vraste  flow  from  i  to  j  at  time  t. 


RS 
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Specifically  the  model  can  be  stated  as  follows: 

Minimize     Z      (CC.,   +  CS .  .  )  (3.1) 

j,k,t      J^^     ^J^ 

Subject  to  the  constraints: 


Z   Z.  .   <  d.  for  all  i  and  t      (3.2) 

j=l  -^'-     -^ 

(v7aste   at   each   source  must   be   satisfied) 

I  K 

E      Z..      <      Z      RC.,  for   all   j    and   t  (3.3) 

i=l      ^J^~k=l        J^^ 

(the  amount  of  waste  flow  into  location  j  must  be  less  than  or  equal  to 

treatment  capacity'  at  location  j) 

Z..   <  RS..  for  all  i,  i,  and  t   (3.4) 

ijt  -   ijt 

(the  amount  of  waste  flow  from  source  i  to  location  j  at  tim.e  t  must  be 

less  than  or  equal  to  its  required  transmission  capacity) 

RS    -  S..   <  SX..  for  all  i,  j,  and  t   (3.5) 

-.jt     ijt  -   ijt 

(the  difference  between  required  capacity  and  existing  capacity  for 

transmission  must  be  less  than  or  equal  to  the  nev;  capacity  to  be  added 

for  the  v;aste  flow  from  source  i  to  location  j  at  time  t) 

RC  ,   -  C  ,   <  ex.,  for  all  j,  k,  and  t   (3.6) 

jkt    3kt  -   3kt 

(the  difference  between  required  capacity  and  existing  capacity  for 

treatment  must  be  less  than  or  equal  to  the  new  capacity  to  be  added 

for  type  k  facility  at  location  j  at  time  t) 

where  CC  ,   and  CS     are  functions  of  CX..   and  SX._,  respectively, 
jkt       ijt  jkt       ijt 

This  model  involves  three  dimensions.  These  are  time,  space,  and 
types  of  facility.  In  order  to  solve  this  model  one  needs  information 
about  the  amount  of  waste  generated  at  each  source,  the  cost  functions 
of  treating  and  transporting  the  waste,  and  the  locations  for  regional 
wastewater  treatment  plants.   In  the  follovvring  sections  we  will  discuss 
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briefly  how  this  infornation  v;as  obtained  for  this  study.   A  brief  dis- 
cussion will  be  given  on  the  basic  ideas  v/hich  have  been  built  into 
this  general  model;  then  a  way  of  solving  this  model  will  be  introduced 
as  a  preview  for  the  next  three  chapters. 

Wastewater  Production 

Uastev/ater  comes  from  four  primary  sources:   municipal  sewage,  in- 
dustrial wastewater,  agricultural  runoff,  and  stom  water  and  urban 
runoff.   Estimation  of  municipal  and  industrial  wastewater  flows  and 
loadings  can  be  done  in  one  of  several  v;ays ,  based  on  knowledge  of  past 
and  future  growth  plans  for  the  community,  sociological  patterns,  and 
land-use  planning. 

Municipal  and  industrial  (M6I)  water  use  rates  are  affected  by 
variables  that  may  be  grouped  into  tvjo  broad  categories  [Department  of 
Water  Resources,  California,  1958]: 

1.  Climatic  factors  such  as  temperature,  rainfall,  wind  speed, 
and  so  on. 

2.  Kan-made  factors  v/hich  further  divided  into  two  groups: 

a)  Residential-related  factors  such  as  economic  level,  educa- 
tion, price  of  water,  family  size  and  age,  metering  and 
sewering. 

b)  Other  urban  related  factors  such  as  greenery,  kind  of  com- 
munity, changing  industrial  water  requirements,  water  pro- 
duction and  use  measurements,  population  served,  and  other 
factors  such  as  worn  flow  meters  and  inadequate  distribution 
systems. 
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Watking  [1968]  points  out  about  the  same  results  in  a  study  of 
sociological  perspective  of  water  consumers  in  South  Florida  households. 

An  aggregate  wastewater  production  function  can  be  specified  by 
the  following  relationship: 

d.^   =   f(w    PN.  ,  M   )  (",  7N 

It        It'    It'   it  u./; 

where: 

d_.^  is  the  amount  of  V7astewater  at  zone  i  at  time  t; 
^it  ^^  ^^^   price  of  water  at  zone  i  at  time  t; 

PN_^^  is  the  population  served  by  the  treatment  p]ant  at 
zone  i  at  time  t;  and 

\t  ^^  ^^^   income  level  at  zone  i  at  time  t. 
In  an  urban  water  demand  study,  Hanke  [1968]  pointed  out  that 
average  household  water  uses  differ  littJe  between  metered  and  unmetered 
areas,  but  sprinkling  uses  and  peak  demand  differ  considerably.   In 
other  words,  household  water  uses  are  price  inelastic  and  sprinkling 
uses  are  price  elastic.   Since  in  most  cases  only  household  wastex^ater 
which  is  not  sensitive  to  the  price  of  water  goes  into  sewage  system, 
equation  (3.7)  may  be  simplified  as  follows: 

^it   =   ^'^P\t'  \t^  (3.8) 

Cost  Functions 


The  cost  of  a  regional  wastewater  treatment  netvjork  is  a  function 
of  a  number  of  variables.  Some  of  these  variables  are  internal  to  the 
technology  of  the  service  in  question.  Others  are  a  function  of  condi- 
tions which  are  unique  to  a  certain  site  and  life  style  which  the  plan- 
ner usually  accepts  as  given,  while  a  third  set  involves  factors  which 
relate  to  the  urban  morphology  (form  and  structure).  This  last  set  of 
factors  will  be  regarded  as  choice  variables  from  the  point  of  view  of 
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the  city  planner  who  is  Involved  in  studying  the  cost  implications  of 
alternative  zoning  and  development  policies  but  will  be  considered  as 
given  in  this  study. 

Emphasis  on  regional  planning  usually  centers  around  optimization 
of  the  collection,  treatment  and  disposal  of  wastewater  in  the  region 
to  meet  a  set  water  quality  goal  at  a  minimum  cost  to  the  region.   A 
host  of  alternative  solutions  must  be  entertained  and  each  alternative 
has  an  associated  cost  function.   These  cost  functions  include  cost  of 
facilities,  construction,  operation  and  maintenance  of  collection,  and 
fi.al  disposal.   Although  economies  of  scale  provide  the  main  incentive 
to  a  regional  wastewater  treatment  system,  it  is  not  always  true  that  a 
single  large  treatment  plant  for  the  entire  region  will  provide  the 
least  cost  system,  nor  is  it  always  possible  to  have  a  single  plant. 
This  is,  of  course,  due  to  the  fact  that  the  location  and  the  capacity 
of  the  treatment  facility  can  not  be  independent  of  other  factors  such 
as  waste  sources  to  be  served,  the  topography  of  the  region,  length  of 
interceptors,  and  locations  and  number  of  existing  sewage  facilities. 
The  final  analysis  in  the  regional  plan  must,  therefore,  consider  trade- 
offs between  combining  treatment  facilities,  combining  and/or  abandoning 
transmission  facilities,  and  final  disposal  methods. 

In  this  study,  cost  functions  of  the  wastewater  treatment  facili- 
ties will  be  divided  into  two  categories-cost  functions  of  collect: 
systems,  cost  functions  of  treatment  plant  facilities  (secondary  and 
advanced).   Furthermore,  each  category  of  cost  functions  will  be  sub- 
divided into  two  subcategories,  i.e.,  construction  cost  functions  and 
operation  and  maintenance  cost  functions. 
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Since  an  optimal  staging  policy  is  also  affected  by  the  rate  of 
interest  and  inflation,  one  dollar  spent  today  will  not  be  considered 
equivalent  to  one  dollar  spent  at  some  future  date.   The  decision  pro- 
blem at  hand  is  a  multipoint-input,  multipoint-output  allocative  pro- 
blem [Henderson  and  Quandt,  1958,  p.  244]  with  the  demand  given  for 
output  treatment  capacity  at  each  point  of  time. 

To  project  year-by-year  operating  costs  for  each  segment  of  the 
proposed  system,  a  present  value  of  future  expenditures  scheme  was  used 
in  the  optimization  model.   A  basic  operating  cost  for  the  present  year 
was  set  for  each  treatment  plant  and  transmission  line.   The  value  was 
then  inflated  at  a  parameterized  rate  to  reflect  increased  labor  and 
maintenance  expenditures  in  the  future.   The  present  value  of  the  re- 
sulting series  of  expenditures  was  then  added  to  the  initial  cost  of 
each  system  component  to  obtain  a  present  value  of  costs  over  the  plan- 
ning horijcon. 

Let  i   be  the  market  rate  of  interest  and  y  be  the  inflation  rate 
t  t 

connecting  marketing  date  t-1  and  t.   Then  the  present  value  of  one 
dollar  payable  at  the  end  of  t   marketing  period  is 

t  _1 

[  n   (1  +  Y  )  (1  +  i  )   ]   =  V^   t  =  1,  ...,  T     (3.9) 

,  T  T  t 

T  =  l 

For  example,  the  total  present  value  of  costs  for  a  type  k  treatment 

facility  at  location  j  at  the  end  of  time  t  is 

CC.,  ,   =   (FC,  ^  +  FM.,  J  V^  (3.10) 

jkt        jkt     jkt    t 

where  FC.,   and  FM  ,   are  the  construction  cost  and  operation  and  main- 
jkt       jkt 

tenance  costs  for  type  k  facility  at  location  j  at  time  t,  respectively. 
This  same  rule  will  be  used  to  calculate  the  present  value  of  a  trans- 
mission system. 
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In  order  to  solve  this  general  model,  we  need  information  about 
the  potential  sites  for  individual  and/or  regional  treatment  plants. 
Theoretically,  a  regional  treatment  plant  can  serve  from  tv;o  to  nine 

cities;  therefore,  for  a  nine-city  problem,  such  as  this  study,  we  end 

9  9 

up  with   E   (  .  )  potential  sites,  where  (  .  )  represents  the  number  of 

i^l 
combinations  of  choosing  i  cities  out  of  nine  cities.   This  study  will 

be  started  from  the  estimation  of  wastewater  production,  finding  po- 
tential sites  for  regional  treatment  plants,  collection  cost  data,  and 
eventually  reaching  a  solution  for  this  general  model. 

The  first  dimension  to  be  considered  in  this  study  is  the  con- 
solidation of  treatment  plants  over  space.   This  dimension  is  included 
because  the  existence  of  scale  economies  may  permit  relatively  large 
treatment  plants  to  process  their  wastewater  at  lower  average  cost  per 
unit.   Thus,  the  cost  for  treating  wastes  from  several  sources  at  a 
regional  plant  may  be  lower  than  the  total  cost  if  wastes  were  treated 
at  each  source.   There  is  a  tradeoff  between  the  gains  from  scale 
economies  and  the  losses  incurred  by  constructing  transmission 
pipelines . 

The  second  dimension  to  be  considered  in  this  study  is  the  con- 
solidation of  treatment  plants  over  time.   Since  the  demands  for  waste- 
water treatment  increase  with  time  and  the  inflation  rate  exceeds  the 
interest  rate  there  may  be  a  net  gain  by  building  a  treatment  plant 
larger  than  necessary.   This  is  the  tradeoff  between  the  gains  in  econo- 
mies of  scale,  savings  in  inflation,  and  the  losses  in  interests  and  in 
operating  and  maintaining  excess  capacity. 

The  third  dimension  to  be  considered  is  the  effect  of  different 
types  of  treatment  facilities.   In  this  study,  we  consider  only  two 
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types  of  secondary  treatment  plant,  three  tertiary  processes,  and  the 
possibility  of  substituting  tertiary  treatment  by  land  treatment. 

The  problem  of  consolidation  of  treatment  plants  over  space  V7as 
first  solved  for  this  problem;  the  central  location  model  mentioned  in 
Chapter  II  was  used.   Then  the  consolidation  of  treatment  plants  over 
time  was  incorporated  with  the  consolidation  over  space,  and  solved 
simultaneously.   Finally,  a  feasibility  analysis  of  substituting  land 
treatment  for  tertiary  treatment  is  presented. 


CHAPTER  IV 
WASTEWATER  TREATMENT  SITES 

In  order  to  develop  a  working  model,  the  following  information  is 
required: 

1.  potential  wastewater  sources  and  potential  sites  for  v/aste- 
water  treatment  plant, 

2.  the  amount  of  v/astewater  produced  at  each  source,  and, 

3.  cost  functions  of  wastevrater  treatm.ent  facilities  as  well 
as  cost  functions  of  transmission  pipeline  from  each  waste- 
water source  to  the  potential  sites  for  wastewater  treat- 
ment facilities. 

This  chapter  will  provide  the  above  information  and  then  compare 
the  potential  treatment  sites  found  in  this  study  with  those  provided  by 
North  Central  Florida  Regional  Planning  Council  [1973,  pp.  11.1-11.27]. 

General  Description  of  the  Area 

Alachua  County  is  located  in  the  North  Central  section  of  the 
Florida  peninsula.   The  Santa  Fe  River,  on  the  north,  separates  Alachua 
County  from  Columbia,  Bradford,  and  Union  Counties.   Putnam  County  lies 
along  most  of  the  Eastern  boundary,  with  Clay  County  at  the  north  end. 
Marion  County  is  at  the  south  end  of  the  east  boundary  and  runs  along 
roughly  two-thirds  of  the  south  boundary.   Levy  County  borders  Alachua 
County  to  the  south  and  west,  and  Gilchrist  to  the  west.   All  of  Alachua 
County's  territory  falls  within  the  29°  25'  and  29°  57'  parallels  of 
latitude  north  and  the  82°  03'  and  82°  40'  meridans  of  longitude  west. 
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The  county  has  a  total  land  area  of  558,320  acres  or  888  square 
miles.   These  figures  do  not  include  the  large  bodies  of  water  repre- 
senting about  76.4  additional  square  miles  [North  Central  Florida  P.e- 
gional  Planning  Council,  1973].   Alachua  County  lies  within  the  central 
highlands  of  the  state.   The  area  generally  ranges  from  level  or  nearly 
level  to  gently  sloping.   Most  ground  elevations  range  betvzeen  50  and 
210  feet  above  sea  level.   Climate  in  Alachua  County  is  mild,  with  a 
mean  annual  temperature  of  70.2  degrees  and  an  average  yearly  rainfall 
of  approximately  52  inches.   Precipitation  varies  from  year  to  year, 
particularly  in  the  summer  m.onths .   Alternating  v;et  and  dry  cycles  of 
several  years'  duration  are  observed.   Most  of  the  rain,  an  average  of 
about  36  inches  per  year,  is  tropical  in  nature,  falling  during  the 
sunmier  months  as  thundershowers .   About  16  inches  of  the  yearly  average 
come  as  winter  rain  of  the  cyclonic  or  frontal  type,  usually  slow  and 
drizzly  and  followed  by  a  drop  in  temperature. 

As  a  result  of  the  increase  in  population  and  its  changing  charac- 
teristics, the  way  in  v;hich  land  is  used  in  the  county  has  tv7o  basic 
land  use  patterns.   Agricultural  land  use  has  retained  its  dominance 
and  is  still  the  primary  land  use  category  as  it  was  in  the  early  1900's. 
The  Gainesville  Urban  Area  has  developed  as  a  major  urban  center  at  the 
approximate  geographic  center  of  the  county.   Land  use  in  the  urban 
center  is  primarily  residential  and  institutional  in  nature.   Commercial 
land  uses,  although  occupying  a  relatively  smaller  number  of  acres,  have 
developed  to  support  the  agricultural,  residential,  and  institutional 
uses.   Table  4.1  summarizes  the  estimated  major  land  use  categories  for 
Alachua  County  in  1972.   The  most  significant  changes  expected  are 


Table  4.1    Estimated  Land  Uses  in  Alachua  County 


Nonagriculture  Agriculture 

Public  Developed  Water  „    ^   ,  .,   ^    „    ^  ^^, 

,  Croplano  Pasture  Forest  Otnar 
Land    Land    Areas 


11.17  20.81      19.16      27.20 


Source:   North  Central  Florida  Regional  Planning  Council  [1973]. 

increase  in  forest  lauds  and  urban  area,  and  a  decrease,  in  croplands  or 
cultivated  areas. 

Alachua,  Archer,  Hawthorne.  High  Springs,  LaCrosse,  Micanopy, 
Newberry,  and  Waldo  are  the  eight  incorporated  centers  of  population  in 
Alachua  County  exclusive  of  Gainesville.   Theii"  location  is  shovjn  on 
Figure  4.1.   Melrose,  an  unincorporated  community  is  situated  on  the 
boundary  line  between  Alachua  and  Putnam  Counties. 

There  are  many  unincorporated  towns  within  Alachua  County.   The 
most  important  are  Arredondo,  Campville,  Cross  Creek,  Earleton, 
Evinston,  Fairbanks,  Forest  Grove,  Grove  Park,  Hague,  Island  Grove, 
Lochloosa,  Monteocha,  Orange  Heights,  Rochelle,  Santa  Fe,  and  Windsor. 
Most  of  the  rural  dwellings  are  found  v;ithin  the  western  half  of  the 
county.   Table  4.2  provides  a  brief  description  of  each  of  the  incor- 
porated cities. 

The  city  of  Gainesville  is  the  only  incorporated  community  with 
near  adequate  wastewater  collection,  treatment,  and  disposal.   Present 
treatment  facilities  in  Gainesville  consist  of  a  5-ngd  high-rate  trick- 
ling filter  plant  and  the  recently  completed  4.5-mgd  contact  stabiliza- 
tion plant.   In  addition,  the  city  ov;ns  a  third  treatment  plant;  this 
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High  Springs 
(1.62,  22.45) 


LaCrosse 
(13.19,  23.67) 


Alachua 
(7-59,  20.24) 


Waldo 
(27.62,  20.05) 


Newberry 
(0.86,  10.05) 


Gainesville 
(18.05,  10.48) 


Hawthorne 
(32.57,  6.62) 


Archer 
(6.14,  2.09) 


Micanopy 
(20.95,  0.83) 


Figure  4 . 1 


Locations  of  the  Incorporated  Cities  in  Alachua  County, 
Florida 
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xs  an  extended  aeration  plant  V7ith  a  nominal  capacity  of  55,000  gallons 
per  day. 

The  city  of  Waldo  is  the  only  other  city  providing  public  sanitary 
facilities;  the  present  v;astewater  treatnent  facility  is  a  40,000  gallon 
septic  tank,  which  is  grossly  inadequate.   The  remaining  eight  incorpo- 
rated municipalities  have  either  privies  or  septic  tanks  with  the  large 
floi-:  users  such  as  hotels,  restaurants,  and  service  stations,  primarily 
located  along  major  thoroughfares,  being  served  by  small  package  plants. 
Table  4.3  gives  an  inventory  of  the  treatment  plants  in  Alachua  County 
by  their  capacities. 


Table  4.3    ISiumber  and  Capacities  of  the  Treatment  Plants  in 
Alachua  County,  December,  1973 


Capacity  (million  gallons  per  day) 


0.01    0.01-0.05    0.05-0.1    0.1-1.0    1.0 

Number         8         14  2  3        3 

of  Plants 

Source:   A  computer  output  provided  by  Alachua  County  Pollu- 
tion Control  Board. 


As  mentioned  in  Chapter  III,  the  amount  of  wastewater  produced  may 
be  estimated  as  a  function  of  the  size  of  population  served,  and  the  in- 
come level.   However,  due  to  the  lack  of  accurate  data  on  the  size  of 
population  served,   this  relationship  v.'as  not  estimated.   Instead,  an 


The  monthly  sewage  treatment  operation  report  submitted  by  each 
treatment  plant  should  provide  this  information.   However,  a  review  of 
those  reports  submitted  by  the  treatment  plants  in  Alachua  County  indi- 
cated incompletenesses  as  well  as  inaccuracies  which  prohibit  meaningful 
estimation  of  wastewater  production  function.   Therefore,  an  average 
daily  flow  per  capita  estimate  was  used. 
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average  daily  flow  per  capita  was  used  to  estimate  the  wastev/ater  pro- 

2 
duced  at  each  wastewater  source.   Assuming  a  sewage  flow  of  135   gallon 

per  capita  per  day,  Table  4.4  shows  the  estimates  of  sewage  flow  of  each 

city  in  Alachua  County  in  year  1975,  1980,  and  1990,  respectively. 


Table  4.4   Estimated  Wastewater  Flows  in  Alachua  County, 

Florida,  by  Incorporated  Cities  for  1975,  1980, 
and  1990 


Wastewater  Flows  (gal/day) 

Name 

1975 

1980 

1990 

Alachua 

456,975 

479,250 

600,750 

(3,385) 

(3,550) 

(4,450) 

Archer 

133,650 

149,850 

182,250 

(990) 

(1,110) 

(1,350) 

Gainesville 

13,729,500 

16,044,615 

20 

,632,185 

Urban  Area 

(101,700) 

(118,849) 

(152,831) 

Hawthorne 

168,750 

177,525 

199,800 

(1,250) 

(1,315) 

(1,480) 

High  Springs 

409,050 

445,500 

532,575 

(3,030) 

(3,300) 

(3,945) 

LaCrosse 

58,050 

65,475 

81,000 

(430) 

(485) 

(600) 

Micanopy 

108,540 

116,100 

120,825 

(804) 

(860) 

(895) 

Newberry 

184,275 

201,150 

239,625 

(1,365) 

(1,490) 

(1,775) 

Waldo 

114,075 

120,825 

134,325 

(845) 

(895) 

(995) 

Total 

15,362,865 

17,800,290 

22 

,723,875 

(113,799) 

(131,854) 

(168,325) 

Remark:   The  figures  in  parentheses  are  estimated  population 
sizes  from  North  Central  Florida  Regional  Planning 
Council  [1973]. 


Miich  is  an  average  daily  flow  of  73  cities  in  27  states  in  th« 
United  States  [Loehr,  1968]. 


AO 


After  potential  v/astewater  sources  and  the  amount  of  wastev;ater 
produced  by  each  source  were  determined,  the  potential  sites  for  v/aste- 
water treatment  facilities  were  considered.   In  order  to  chose  the  po- 
tential sites,  the  central  location  model  mentioned  in  Chapter  II  was 

3 
used.   The  number  of  potential  sites  is  quite  large.    As  a  consequence, 

the  large  amount  of  computer  time  makes  it  impractical  to  find  a  solu- 
tion for  each  potential  site  using  tlie  model  presented  in  the  previous 
chapter.   To  overcome  this  problem,  a  computer  program  v/as  developed 
(see  Appendix  I  and  next  chapter  for  the  details  of  this  program)  to 
scan  through  all  the  potential  sites,  and  eliminate  those  vmich  have 
higher  sum.  of  regional  treatment  and  transmission  costs  than  the  sum  of 
individual  treatment  costs.   This  computer  program  used  a  central  loca- 
tion model  to  determine  potential  sites  and  to  calculate  the  transmis- 
sion costs  and  the  combined  treatment  costs.   The  program  then  compares 
the  sum  of  transmission  and  treatment  costs  v/ith  the  sum  of  the  treat- 
ment costs  if  each  city  involved  built  its  own  treatment  facilities. 

Wast ev7a ter  Transmission  Costs 

The  cost  functions  for  transmission  are  as  follov;s: 

1.   The  relationship  between  flow  and  the  diameter  of  sewers,  which 
minimizes  the  sum  of  the  cost  ot  frictional  pov/er  and  the  debt  service 

for  the  pipe,  estimated  in  [Smith  and  Eilers,  1971]  as: 

.463 
Economic  Diameter  (inches)   =   8.55Q'  (4.1) 

where: 

Q   =   designed  flow  in  mgd. 


3  ^    9 

The  total  number  of  potential  site  equals  to   Z   (  .  )  =  491,  where 

(  .  )  is  defined  as  9J .  i=l 

'  (9-i)!  (i)! 
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2.  The  cost  associated  with  sewers  which  includes  the  amortization 
cost  and  a  small  maintenance  cost  in  terms  of  October,  1970  dollars  v;as 
estimated  [Smith  tmd  Eilers,  1971]  as  follows: 

Construction  costs  in  dollars  per  mile   = 

1540. 7(ID  +  2.0436)^-^''^^^  (4.2) 

where: 

ID   =   inside  pipe  diam.eter  in  inches. 

The  cost  of  maintaining  the  sev7er  was  taken  as  $40/yr/m.ile. 

3.  The  construction  cost  for  pumping  stations  has  been  estimated 
[Smith  and  Eilers,  1971]  as: 

Construction  cost,  $  =   76,:-i00Q'  (4.3) 

where: 

Q  =  average  flow  in  mgd. 
The  operation  and  maintenance  cost  for  pum.ping  stations  has  been  esti- 
mated [Smith  and  Eilers,  1971]  as: 

cents/lOQO  gal.   =   1.59Q~'^^"^  (4.4) 

v^7here  electrical  power  was  assumed  to  cost  one  cent  per  kw-hr.   The 
hydraulic  efficiency  of  the  pump  was  taken  as  60%  and  the  electrical 
efficiency  of  the  driving  motor  v/as  taken  as  80%.   All  costs  were  keyed 
to  October,  1970. 

The  Ten  State  Standards  [Great  Lakes — Upper  Mississippi  River  Board 
of  State  Sanitary  Engineers,  1971]  require  interceptor  sewers  to  be 
sized  to  carry  3.5  times  the  design  flov;  for  the  plant.   This  factor  was 
used  for  sizing  both  the  pipelines  and  the  pum.ping  stations,  since  the 
flow  of  sanitary  sewage  reaching  the  treatment  plant  varies  over  time, 
and  the  capacity  of  the  sewage  facility  is  fixed  over  the  short  run. 
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Excess  capacity  (capacity  above  the  average  daily  flow)  is  necessary  to 
collect  and  treat  the  sewage. 

The  collection  and  transmission  of  peak  flov;s  may  be  handled  in  one 
of  the  three  ways:   (1)  sewers  could  be  built  to  handle  the  expected 
peak  or  some  greater  flow,  (2)  a  holding  tank  could  be  installed  to 
average  flows,  and  (3)  some  method  could  be  employed  which  V70uld  average 
the  release  of  flov;s.   Of  these  three,  engineers  consider  the  first  as 
the  only  feasible  alternative.   Holding  tanks  have  not  been  used  because 
there  are  design  problems  v/hich  engineers  feel  prevent  their  use. 
Averaging  at  the  source  of  the  flovjs  is  an  alternative  which  engineers 
avoid  completely. 

Wastewater  Treatment  Costs 

Treatuent  of  sewage  is  the  process  of  removing  undesirable  materi- 
als from  the  water  and/or  changing  them  into  less  objectionable  forms. 
In  the  typical  plant  this  is  done  by  a  combination  of  physical  and  bio- 
logical processes.   In  essence,  the  treatment  process  uses  a  more  rapid 
version  of  the  same  process  v.'hich  would  occur  if  the  v/astes  were  re- 
leased to  the  receiving  water  without  treatment. 

The  selection  of  a  wastewater  treatment  process  or  a  combination  of 
processes  v/ill  depend  upon:   (1)  the  characteristics  of  the  wastewater, 
(2)  the  required  effluent  quality,  (3)  the  costs  and  availability  of 
land,  and  (4)  the  future  upgrading  of  water  quality  standards.   In  this 
study,  two  conventional  treatment  system.s  will  be  considered,  i.e.. 
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high-rate  trickling  filter  system   and  activated  sludge  system.   Since 
the  high-rate  trickling  filter  systems  cannot  remove  more  than  about 
80%  of  BOD  and  SS,  most  new  plants  for  domestic  wastewater  treatment  are 
activated  sludge  plants. 

Of  course,  there  are  minor  variations  in  the  designs  of  treatment 
plants.   For  example,  there  are  different  ways  to  circulate  the  primary 
as  well  as  secondary  sludge,  and  there  are  different  designs  for  re- 
turned sludge  entering  the  primary  settling  tank  and  the  like.   Hov;ever, 
the  cost  variations  incurred  by  these  different  designs  are  not  large, 
hence  these  differences  in  cost  vjill  not  be  considered  in  this  study. 

There  are  several  ways  one  can  estimate  the  construction  and  opera- 
tion and  maintenance  costs  of  a  treatment  plant.   A  cost  estimation  com- 
puter program  was  developed  by  Smith  [Environmental  Protection  Agency, 
1973b]  which  allows  the  consideration  of  influent  waste  strength  as  well 
as  the  different  combinations  of  treatment  processes.   Such  a  program 
was  purchased  from  National  Technical  Information  Service  and  adapted  to 
the  IBM  370  system  at  the  University  of  Florida.   However,  there  were 
errors  in  the  program  which  precluded  its  use. 

EPA  [1973a]  has  estimated  cost  equations  with  regression  analysis 
using  data  on  accepted  bids  for  construction  of  new  niunicipal  wastewater 


4 
The  differences  of  design  between  a  high-rate  trickling  filter  sys- 
tem and  a  standard-rate  (or  low-rate)  trickling  filter  system  depend  on 
their  hydraulic  loading,  200-1000  gpd  per  sq.  ft.  for  high-rate  opera- 
tion, and  25-100  gpd  per  sq.  ft.  for  standard-rate  operation  [Fair  et  al. 
1971]. 

Treatment  plant  design  examples  provided  by  National  Technical  In- 
formation Service  were  used  to  test  the  computer  program,  however,  the 
computer  output  showed  that  during  the  calculation  some  logarithm  and 
exponential  functions  had  illegal  arguments.   Several  attempts  have  been 
made  without  success  to  correct  these  errors,  hence  the  program  v."as  not 
used  in  this  study. 
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treatment  plants  in  several  states  between  1967  and  1969  as  a  function 
of  design  flow.   Costs  were  updated  to  September  1972  dollars  by  using 
the  EPA  Sev/age  Treatment  Plant  Construction  Cost  Index.   EPA's  estimates 
provide  another  means  of  estimating  costs.   Their  estimates  were  as 
follows : 

1.  For  high-rate  trickling  filter  system: 

AC  =   852049.58  m'"^^'^^"^  (4.5) 

2.  For  activated  sludge  system: 

AC  =   699835.09  m"""^^^^  (4.6) 

where: 

AC   =   average  construction  cost,  and 
M   =   amount  of  design  flow  in  mgd . 
And  estimates  of  operation  and  maintenance  cost  collected  from  re- 
ports on  600  plants  betv/een  1968  and  19  70  in  January  1968  dollars  have 
the  following  relationships  [Michel  and  Johnson,  1970] : 

1.  For  high-rate  trickling  filter  system: 

TC  =   31959.50  M'^^^^  (4.7) 

2.  For  activated  sludge  system: 

TC  =   46989.41  M'^^^^  (4.8) 

where: 

TC   =   total  annual  operation  and  maintenance  cost,  and 
M   =   amound  of  design  floxs?  in  mgd. 
A  third  way  of  estimating  wastewater  treatment  costs  is  the  appli- 
cation of  the  six-tenths  rule  [Berthauex,  1972;  Williams,  1947;  Chilton, 
1950;  Eckenfelder  and  Ford,  1969].   Under  this  rule,  the  cost-capacity 
relationship  has  the  following  form: 


^5 


Cost  of  Plant  A  _   , Capacity  of  Plant  A  ,M 
Cost  of  Plant  B   ~    Capacity  of  Plant  B 


Dab        a 


where: 


C   =   cost  of  an  item  of  capacity  Q  , 
a  '     ^  ^a' 

M 
K  =   base  cost  factor  and  equals  C,  /Q,  ,  and 

b   b 

M  =  measure  of  the  economy  of  scale. 
The  name  "Six-tenths  rule"  arose  because  the  exponent  M  is  often  about 
0.6  for  chemical  processing  plants.   The  rule  not  only  applies  to 
scaling-up  the  costs  of  an  entire  plant  but  shows  that  the  exponent 
varies  vjith  the  type  of  plant  under  construction.   Some  of  the  estimates 
of  factor  K  are  shown  in  Table  4.5. 


Table  4.5   Cost-Capacity  Factors  for  Municipal  Wastewater 
Treatment  Plants 


Size  Factor  (M) 
Type  of  Plant 


(1)  (2) 


Activated  Sludge  0.77  0.65 

Trickling  Filter  0.60  0.62 

Source:   (1)  from  Berthauex  [1972].   (2)  from  EPA  [1973a]. 

In  this  study,  the  estimated  equations  (4.1)  through  (4.8)  by  EPA 
were  used.   All  the  cost  figures  v.'ere  adjusted  to  September  1974 
dollars.   Equations  (3.9)  and  (3.10)  were  used  to  calculate  the  present 
values  of  1975,  1980,  and  1990  cost  figures  in  terms  of  September  1974 
dollars. 
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Tertiary  Treatme.nts 

Tertiai-y  treatment  may  be  categorically  defined  ar.  treatnient  for 
tlie  removal  csf  pollutants  not  removed  by  conyentional  treatment  pro- 
cesses (activated  sludge,  trickling  fil'^ers,  etc.)-   These  pollutants 
v;ill  include  suspended  solids  (SS) ,  biochemical  oxygen  denictnd  (HOT)), 
refractory  organics  (u;-ually  reported  as  chemical  oxygen  demand  (COD) 
or  total  organic  carbor'.  (TOC)),  nutrients  (nitrogen  and  phospLcrun)  ^  and 
inorganic  salts.   In  the  United  States  today,  increasing  emphiir.is  is 
being  placed  on  the  rem.oval  of  phosphorus  and  the  removal  of  unoxi.dilzed 
nitrogen  vhlch  will  exhibit  a  long-term,  oxygen  demand  in  the  receiving 
waters. 

The  lime  clarification  process  is  used  primarily  for  the  removal 
of  phosphcras  and  suspended  organic  matter.   An  additional  benefit  is 
the  incrcc-.sing  pH  resulting  from  lime  addition  which  makes  amiiionia 
nitrogen  available  for  removal  by  air  stripping.   Even  thougn  m.any  pro- 
blems associated  with  the  use  of  ammonia  stripping  processes  reoiain  to 
be  solved,  it  is  presently  viev/ed  as  the  most  promising  process  for  re- 
moving ammonia  nitrogen  from  wastewater.   A  process  employed  in  T.ake 
Tahoe  is  shown  in  Figure  4.2.   The  effluent  characteristics  are  shoi<ra  in 
lower  part  of  the  same  figure. 

The  estim.ated  cost  equations  for  three  tertiary  treatment  processes 
to  be  u;;ed  in  this  study  from  the  data  provided  by  [Smith  and  HcMichacl, 
1969]  are  presented  below: 

1.   Tv70  clarifier  lime  clarification  process  without  chemicals: 

There  is  some  evidence  that  a  coagulant  aid  such  as  iron  might  be 
required  in  the  second  clarifier.   Since  the  need  for  this  chemic:il  is 
not  clearly  established,  it  was  not  included  in  the  cost. 
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Const  rue  Lion  cost  (niillioa  dollars) 

. 77610^ 
0.128371Q  (4.9) 

Om   cost  (cic-llars/day)   =   62. 851SQ"  ^^"^^-^^^        (4.10) 

2.  Lime  recal cination  plus  make  up  lime  for  use  with  lime 
clarification: 

Construction  cost  (million  dollars)   = 

0.196977Q-'°^^°''  (4.11) 

O&M  cost  (dollars/day)   =  49. 37cSQ"  ^^^^"^"^        (4.12) 

3.  Ammonia  stripping  of  lime  clarified  v;astewater: 

Construction  cost  (million  dollars)   = 

0.490535Q^-^'^'^'^  (4.13) 

7  7  R^ R7 
O&M  cost  (dollars/day)   =  40.9442Q*  (4.14) 

where  Q  represents  design  capacity  measured  by  miillions  of  gallons  per 

day  (ragd) . 

As  pointed  out  in  Chapter  II,  the  long-run  cost  curve  is  the 

envelope  of  the  short-run  cost  curves.   Cost  equations  (4.5)  through 

(4.8)  were  estimated  by  regression  analysis  from  cross-section  data,  and 

these  cost  observations  v/ere  subject  to  variations  in  the  age,  type,  and 

cost  of  equipment,  the  quality  of  executives  as  well  as  the  effects  from 

spatial  variation  in  factor  prices.   Cost  equations  (4.2)  through  (4.4) 

and  (4.9)  through  (4.14)  were  estimated  as  the  loci  of  minimum  costs  by 

the  engineering  approach.   Therefore,  equations  (^.  1)  through  (4.14)  as 

used  in  this  study  to  estimate  the  wastewater  treatment  costs  give 

higher  cost  estimates  than  the  minimum  theoretical  costs  from  eqi:ation 

(2.4d). 
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Assumptions  about  Cost  Increases 

In  order  to  estimate  the  rates  of  cost  increases  over  time  to  be 
used  in  the  calculation  of  present  values,  regression  analysis  was 
used.   In  this  analysis,  we  assume  that  cost  index  is  a  function  of  time, 
i.e., 

Cost  index  =  a  +  bt  +  u^  t  =  1966,  ...,  1973 
v/here  a  and  b  are  parameters  to  be  estimated,  and  u  is  the  disturbance 
term  at  time  t. 

Three  different  kinds  of  cost  indexes  were  used  in  this  analysis, 

they  are:   sewer  construction  cost  index  (SVJ) ,  sewage  treatment  plant 

construction  cost  index  (ST),  and  consumer  price  index  for  residential 

water  and  sewerag;e  services  (SR)  .    All  the  September  cost  index  figures 

from  1966  to  1973  were  used  for  the  estimation  of  parameters  a  and  b, 

and  the  results  are  as  follows: 

SW   =  -23567.0  +  12.0444t 
^  (1990.00)  (1.0104) 

ST   =  -20378.8  +  10.4214t 
^      (1585.40)  (0.8050) 

SR   =  -14604.9  +  7.4762t 
(973.204) (0.4941) 

where  the  "hats"  on  top  of  SW  ,  ST  ,  and  SR   represent  estimated  values, 
^      t'    t        t    ^  ' 

the  figures  in  the  parenthesis  are  estimated  standard  errors  of  the 

2 
parameters,  and  R  's  are  multiple  determination  coefficients. 

Parameter  b  represents  the  effect  of  one  unit  change  in  t  on  the 

corresponding  cost  index,  hence  it  can  be  explained  as  the  annual  rate 

of  increase  in  costs  and  they  were   used  in  equations  (3.9)  and  (3.10). 


r2  . 

.9595 

(4.15) 

r2  = 

.9654 

(4.16) 

9 
R"   = 

.9745 

(4.17) 

These  three  cost  indexes  vjere  provided  by  Advanced  Waste  Treatment 
Research  Laboratory,  Cincinnati,  Ohio. 
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h  in  (4.15)  was  used  in  the  calculation  of  sewer  construction  costs,  b 
in  (4.16)  was  used  in  the  claculation  of  sewage  treatvaent  plant  con- 
struction cost,  and  b  in  (4.17)  was  used  in  the  calculation  of  Of-M  costs 
for  both  sewer  and  sewage  treatment  plant.   Service  lives  for  pumping 
stations  and  wastewater  treatment  plants  were  assumed  to  be  25  years,  and 
for  transmission  pipelines  50  years.   Interest  rate  v;as  assumed  to  be 
5%. 

Potential  Sites  for  Wastewater  Treatment  Plants 

In  order  to  reduce  the  possible  number  of  potential  sites,  a  com- 
puter program  vzritten  in  FORTRAl.'  IV  was  developed  (Appendix  I).   This 
computer  program  calculates 

1.  the  cordinates  in  miles  from  the  origin  (29°  30'  parallel 
of  latitude  north  and  82°  37'  30"  meridian  of  longitude 
west)  for  each  of  the  cities  of  this  study  (Table  4.6). 

2.  construction  costs  and  operation  and  maintenance  costs  for 
wastewater  treatment  plants,  transmission  pipelines,  and 
pumping  stations  for  each  city  using  equations  (4.1) 
through  (4.14),  and 

3.  the  transmission  cost  from  each  source  to  the  potential 
sites  obtained  using  equations  (2.5)  and  (2.6).    Finally, 
the  program  compares  the  sum  of  pumping  station  costs, 
transmission  costs,  and  regional  wastewater  treatment 
costs  for  each  site  with  the  wastewater  treatment  costs 
for  a  completely  disaggregated  system  (i.e.,  one  in  v;hich 
each  of  the  nine  cities  build  their  ovm  facilities).   If 
the  regional  cost  is  less  than  or  equal  to  the  sum  of  in- 
dividual costs,  the  program  prints  out  the  related  cost 
figures. 


w-  represents  transmission  pipeline  cost  per  mile  from  city  i  as 
calculated  from  equations  (4.2)  or  (4.4). 
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Table  4.6   Cordinates  of  Municipalities  in  Alachua  County  in 
Miles^ 


^,   1,   b 

City  Number 


Cordinate 


East(x)  North (y) 


High  Springs  1  1.6190  22.4526 

Newberry  2  0.8572  10.0476 

Archer  3  6.1428  2.0952 

Alachua  4  7.9524  20.2380 

Gainesville 

Urban  Area  5  18.0476  10.4762 

Lacrosse  6  13.1904  23.6665 

Micanopy  7  20.9524  0.3810 

Waldo  8  27.6190  20.0476 

Hawthorne  9  32.5714  6.6190 


These  cordinates  are  calculated  from  the  origin  80°  37' 
30"  longitude  West  and  29°  30'  latitude  North  in  miles. 

These  numbers  will  be  used  later  to  represent  the  cor- 
responding cities. 


9 
It  is  assumed  that  force  mains  were  used  as  the  connecting  pipe- 
lines between  contributing  cities  and  receiving  regional  plants,  and  a 
raw  sewage  pumping  station  has  been  sized  to  handle  the  flow  at  each 


9 
Smith  and  Eilero  [1971]  indicate  that  the  cost  of  constructing 

force  mains  is  not  significantly  different  from  the  cost  of  constructing 
gravity  sewers,  hence,  equation  (4.2)  may  be  used  to  find  the  construc- 
tion cost  for  force  mains.   However,  by  a  conversation  with  Mr.  Baldwin 
(Agricultural  Engineer),  this  approach  may  overestimate  the  cost  of 
force  mains. 
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The  Newton-Raphson  method  was  first  used  to  solve  equations  (2.5) 
and  (2.6)  simultaneously,  but  this  method  required  too  many  operations 
to  finish  one  iteration,  hence,  it  v/as  not  used  in  this  study.   Instead, 
an  iterative  procedure  suggested  by  Kuhn  and  Kuenne  [1962]  and  by 
Cooper  [1963]  V7as  used  to  solve  these  two  equations  for  x  and  y  in  terras 
of  w. ,  X . ,  y . ,  and  d . : 

1     X     1  X 

w .  X .      w . 

X   =   E.  -^-  /  T..  ~   ,  (4.18) 

X    d .      X  d . 

X  X 

W.  V .       w. 

y  =   Z.  -V^  /    T..  -^  ,  (4.19) 

X   d.      X  d. 

X  X 

d.   =   ((x^  -  x)^  +  (y.  -  y)^)^  (4.20) 

The  value  of  d.  is  recalculated  via  (4.20)  and  the  procedure  repeated 

X 

until  successive  differences  betvjeen  values  of  x  and  between  values  of 
y  are  negligible.   In  this  study,  the  initial  values  of  x  and  y  for 
the  calculation  of  d.  are  the  weighted  average  of  x.'s  and  y.'s, 
respectively.   As  mentioned  by  ReVelle  et  al.  [1970]  this  ro.ethod  has 
been  found  to  be  extremely  fast,  usually  terminating  at  the  global  mini- 
mum in  less  than  ten  iterations,  and  no  problems  of  lack  of  convergence 
have  been  reported.   Therefore,  a  maximum  number  of  fifty  iterations  was 
set  and  a  four-digit  accuracy  was  required  for  the  solution  of  (x,  y) . 
The  last  (x,  y)  in  the  iteration  procedure  was  taken  as  the  solution  for 
(2.5)  and  (2.6). 

A  regional  treatment  plant  in  this  study  is  defined  as  one  serving 
more  than  one  city.   Let  R(.)  represent  a  regional  treatment  plant  which 
serves  the  cities  represented  by  the  numbers  within  the  parentheses;  for 
example,  R(l,  2)  is  a  regional  treatment  plant  that  serves  cities  1  and 
2  (the  number  for  each  city  has  been  assigned  in  Table  4.6).   If  there 
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is  one  number  within  the  parentheset; ,  it  indicates  that  a  treatment 
plant  serves  only  one  city. 

Assuming  that  ^^^575,  '^iggQ^    ^^^   ^1990  ''^"^^  satisfied  in  years  1975, 
1980,  and  1990,  respectively.  Table  4.7  shov/s  the  possible  sites,  capa- 
cities, and  costs  of  regional  treatment  plants,  the  cities  they  will 
serve,  and  the  years  they  are  to  be  built.   There  are  regional  treatment 
plant  sites  in  the  southeast  corner  of  Alachua  County  other  than  those 
presented  in  Table  4.7  v/hich  were  not  considered  because  of  the  lakes 
and  swamps  in  that  area.   These  sites  involve  the  following  regional 
treatment  plants:   R(6,  7,  8),  R(6,  7,  8,  9),  R(7,  8,  9),  R(2,  3,  4,  5, 
6,  7,  8,  9),  and  R(l,  2,  3,  4,  5,  6,  7,  8,  9). 

Based  on  the  population  estimates  and  wastewater  production  of  135 
gallons  per  person  per  day  [Loehr,  1968],  demands  for  wastewater  treat- 
ment and  the  corresponding  treatment  costs  are  calculated.   The  treat- 
ment cost  figures  in  column  five  in  Table  4.7  are  the  future  service 
values  (see  next  chapter  for  definition)  up  to  1990  for  the  correspond- 
ing treatment  plants.   From  these  estimates,  one  finds  that  (1)  high- 
rate  trickling  filter  plants  are  cheaper  than  activated  sludge  plants, 
(2)  the  costs  for  tertiary  treatment  processes  are  higher  than  the  costs 
of  both  high- rate  trickling  filter  plants  and  activated  sludge  plants. 

Assuming  that  secondary  treatment  is  adequate  to  meet  the  require- 
ments, there  are  several  vjays  to  satisfy  the  demands  in  the  county.  If 
secondary  treatment  is  adequate  to  meet  the  requirements  the  least  cost 
method  of  treatment  is  for  each  city  to  build  a  high-rate  trickling 


Federal  regulation  requires  a  90%  removal  of  BOD's  and  SS  which 
high-rate  trickling  filter  plant  cannot  achieve.   Hov;ever,  if  tertiary 
treatment  processes  were  added  to  this  system,  the  results  of  trickling 
filter  plant  from  this  study  are  still  valuable.   Therefore,  the  results 
of  trickling  filter  system  V7ere  presented  and  discussed. 
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filter  plant  to  meet  its  deinand  (Table  4.8).   Ho'./ever,  if  tertiary 
treatment  processes  are  required  to  meet  the  requirements,  the  results 
show  that  High  Springs,  Nev7berry,  and  Archer  should  cooperate  in  develop- 
ing a  regional  facility  if  their  objective  is  to  minimize  cost.   The 
potential  site  has  a  cordinate  of  :■:  =  0.88  and  y  =  10.07,  which  is  in 
Newberry,  and  an  estimated  regional  treatment  cost  of  4,085,362 
September  1974  dollars  to  satisfy  their  1975  demands  and  to  meet  their 
demands  in  1980  and  1990  by  building  individual  plants.   The  other  six 
municipalities  should  build  separate  treatment  plants  to  meet  their 
individual  demands. 

Concluding  Remarks 

As  mentioned  in  Chapter  III,  this  is  a  study  over  time  and  space. 
This  chapter  discussed  the  consolidation  of  V7aste  treatment  facilities 
over  space  for  specific  points  in  time,  i.e.,  1975,  1980,  and  1990, 
respectively.   There  are  other  ways  that  wastewater  treatment  demands 
can  be  satisfied  over  time.   For  example,  the  estimated  treatment  de- 
mands of  High  Springs  are  .4090  ragd,  .4455  mgd,  and  .5236  mgd  for  years 
1975,  1980,  and  1990,  respectively;  these  demands  could  be  satisfied  by 
building  a  .5326  mgd  treatment  plant  in  1975.   Of  course,  there  are 
other  options.   Chapter  V  brings  the  tim.e  dimension  into  the  optim.iza- 
tion  process  and  considers  the  simultaneous  determination  of  costs  over 
both  time  and  space  dimensions. 

Since  interest  is  in  minimizing  costs,  combinations  over  time  and 
space  which  lead  to  greater  costs  than  those  presented  in  Table  4.8  need 
not  be  considered.   Therefore,  the  results  of  Chapter  V  help  to  narrow 
the  search  for  an  optimal  county  solution. 
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Table  4.S    Possible  Op  Lions  for  Satisfying  County's  Wastewater  Treat- 
ment Demands  and  Associated  Costs,  1975  to  1990^ 


Option  Total  Cost  for  County 


I.   Without  Tertiary  Processes 

i.   high-rate  trickling  filter  system 

9 

Z  R(i)  12,276,820 

i=l 
ii.   activated  sludge  system 

9 
R(l,  2,  3,  A)  +  Z  R(i)  14,464,670 

i=5 
R(l,  2,  3,  4,  5,  6,  7,  8)  +  R(9)  14,550,152 

II.   VJith  Tertiary  Processes 

i.   high-rate  trickling  filter  system 

9 
R(l,  2,  3)  -I-   E  R(i)  44,589,425 

1=4 

9 
R(l)  +  R(2,  3,  4)  +  Z  R(i)  44,590,757 

i=5 
9 
R(l)  +  R(2,  3)  +  Z   R(i)  44,765,688 

i=4 

ii.   activated  sludge  system 

7 
R(l,  2,  3)  +  Z  R(i)  +  R(8,  9)  46,719,711 

i=4 

7 
R(l)  +  R(2,  3)  +  Z  R(i)  +  R(8,  9)  46,986,390 

i=4 
7 
R(l,  2)  +  Z   R(i)  +  R(8,  9)  47,015,999 

i=3 

7 
R(l,  2,  3,  4)  +  Z   R(i)  +  R(8,  9)  45,698,220 

i=5 
R(l,  2,  3,  4,  5,  6,  7,  8)  +  R(9)  46,454,733 

7 
R(l)  +  R(2,  3,4)+  Z   R(i)  +  R(8,  9)         46,698,720 
i=5 


"" Assuming  that  d^   3,  d^^g^,  and  d^^^Q  were  satisfied  in  years  1975^ 
1980,  and  1990,  respectively. 


CHAPTER  V 
CAPACITY  EXPANSION  SCHEME 


Investment  Over  Time 


When  a  firm  experiences  a  growing  demand  for  its  output  and  makes 
investments  aimed  at  increasing  its  capacity  in  existing  lines  of  ac- 
tivity, the  outlays  are  called  induced  investments.   When,  on  the  other 
hand,  the  firm  diversifies  into  new   lines  of  activity,  the  capital  in- 
vestments involved  are  referred  to  as  autonomous  Investments  [Dernburg 
and  McDougall,  1968].   This  study  deals  only  v.'ith  induced  investments  of 
wastewater  treatment  facilities. 

The  simplest  feedback  theory  of  how  induced  investment  takes  place 
is  called  the  acceleration  principle.   This  principle  holds  that  the 
rate  at  which  a  firm  invests  in  new,  as  opposed  to  replacement,  equip- 
ment is  a  linear  function  of  the  rate  of  change  of  output.   Suppose  a 
wastevjater  treatment  plant  has  a  capacity  of  1.0  mgd,  and  is  operating 
at  that  capacity.   Assuming  that  the  only  way  to  increase  treatment 
capacity  is  to  add  more  capacity,  then  as  demand  increases  capacity  must 
be  increased  in  direct  proportion. 

The  acceleration  principle  assumes  that  capacity  is  v;ell  defined 
and  that  when  it  is  reached  such  alternatives  as  overtime,  subcontract- 
ing, and  back  ordering  are  not  available.   The  principle  purports  to 
explain  only  new  investment  or  net  investment  and  also  assumes  that  the 
firm  is  able  to  obtain  the  funds  to  finance  the  indicated  expansion. 
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Wlien  demand  is  constant  the  firm's  investment  raay  be  confined  to 
replacement  of  its  existing  capacity,  although  if  this  was  not  being 
fully  utilized  the  firm  might  actually  neglect  replacement.   VHiile  the 
acceleration  principl;-  provides  a  first  indication,  effective  implemen- 
tation of  an  expansion  program  requires  considerably  greater  insight. 

Capacity  is  not  a  particularly  well-defined  concept  in  most  firms. 
There  is  a  graduaD-  growth  in  unfilled  orders  and  overtime,  shifts  are 
added,  and  any  slack  or  hidden  capacity  comes  into  use.   Thus  management 
has  several  alternatives  to  capacity  expansion.   It  can  use  inventories 
as  a  demand  buffer  to  insulate  to  some  extent  the  steady  operation  of 
production  facilities  from,  changes  in  the  market.   Negative  inventories 
or  back  orders  are  a  useful  means  of  meeting  v;hat  appear  to  be  temporary 
increases  in  demand.   The  capacity  of  existing  facilities  can  also  be 
extended  by  v7orking  overtime  or  expanding  the  working  force.   Sometimes 
capacity  can  be  expanded  by  subcontracting  part  of  the  production  pro- 
cess to  others  aiid  effectively  buying  rather  than  making  some  of  the 
required  output. 

In  the  planning  and  designing  of  a  wastewater  treatment  facility 
for  a  town,  city,  or  an  area  v.-'here  wastev/ater  treatment  requirements 
are  expected  to  increase  with  time,  the  questions  as  to  the  initial  size 
of  the  treatm.ent  plant  and  the  timing  of  capacity  additions  and/or  re- 
placements over  the  period  of  study  or  time  horizon  have  to  be  answered 
in  the  context  of  an  optimal  staging  policy.   Such  a  policy  is  affected 
by  the  v/astewater  treatment  requirement  and  its  growth  rate;  the  rates 
of  interest  and  inflation;  the  cost  of  the  treatm.ent  plant  and  its 
operation  and  maintenance;  the  load  factor,  service  life,  and  expected 
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salvage  value  of  the  treatme-at  plant;  and  the  staging  efficiency  of  the 
system  to  be  designed. 

Incremental  savings  in  construction  and  operation  and  maintenance 
costs  due  to  economies  of  scale  make  it  desirable  to  bear  the  cost  of 
overcapacity  until  demand  catches  up.   A  major  decision  variable  in 
public  and  private  wastewater  treatment  plant  staging  policy  is  the 
amount  of  excess  capacity  to  be  built  initially  into  a  new  system  and 
the  staging  of  capacity  additions  and/or  replacem.ent  (as  the  old  plants 
become  uneconomical  to  run,  being  past  their  useful  service  life)  to 
meet  demands  increasing  with  time. 

For  example,  the  vertical  axis  in  Figure  5.1  represents  the  long- 
run  average  construction  cost  for  vjastewater  treatment,  and  the  hori- 
zontal axis  represents  the  demand  for  X'jastewater  treatment  capacity. 
Q-iqjr    and  Qn  goQ  represent  the  required  capacities  in  1975  and  1980  re- 
spectively; LL  is  the  long-run  average  construction  cost  curve.   Assume 
there  are  two  choices  to  chose  from,  i.e.,  to  build  a  single  plant  v%Tith 
capacity  Q-,qoq  in  1975,  or  to  build  one  plant  with  capacity  Q-,q-,c  in 

1975  and  another  one  with  capacity  Q,r,Q„  -  Q-,mc  in  1980.   The  differ- 

lyou    ly  /  D 

ence  in  construction  costs  represented  by  present  value  is 

D   =  ocdj  -  obeh  -  oafg  ^^   '^  ^\    ,  (5.1) 

""  (1  +  i)^ 

where  oa  =  oc  --  ob,  A  represents  rate  of  inflation  in  construction  cost, 

i  represents  interest  rate,  and  D   represents  the  difference  in  construc- 

c 

tion  costs  between  these  two  options.   If  D   >  0,  then  there  is  a  saving 

in  construction  cost,  if  D   <  0,  then  there  is  a  loss  in  the  two-plant 

c  ^ 

option.   However,  the  choice  between  a  single-plant  option  and  two-plant 
option  depends  not  only  on  construction  costs  but  on  operation  and 
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Cost  ($) 


Capacity  (mgd) 


'1975 


1980 


Figure  5.1   Long-run  Average  Construction  Cost  Curve 


^1975   ^1980 


Capacity  (mgd) 


Figure  5.2   Average  Operation  and  Maintenance  Cost  Curve 
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maintenance  (05;M)  costs  also.   In  Figure  5.2,  L'L'  represents  long-run 
average  O&M  cost  curve,  which  is  the  envelop  of  short-run  cost  curves 
^a'  ^b'  ^c'  ^°^  capacities  Q^^^^  -  Q^^^^,  q^,^^.    and  Q^^g^,  respectively, 
od  equals  the  difference  between  Qiq-,c  31""^  ^iQfin"   ''•'"'^  present  values  of 
O&M  cost  for  the  two  options  mentioned  above  follow. 


f  1980 
Single-plant  option:   D   =   present  value  of  L     Q"f   dQ      (5.2) 


TXvzo-plant  option:   D   =   present  value  of 


1979 
+   Odeg  '^  (5.3) 


where  f  and  f,  are  functions  of  operating  capacity  Q.   Let 
c      b 

D    =   D^  -  D„  (5.4) 

om      12 

then  if  D  +  D    is  negative,  the  single-plant  option  will  be  chosen, 
c    om 

other\'.'ise   the   tvjo-plant    option  v/ill  be   chosen. 

For  each  potential  site  the  model  may  be  stated  as 
Model  1 

''^^  ^^1975^^1975^  -"   ^^980^^1980^  +  ^^990^^1990^ 

Ql975  ^^1975 

^1975  "^  ^1980         -  ^^1980 

^1975  "^  ^1980  "^  ^1990  -  ^1990 

^1975'  ^1980'  ^1990   -  ° 

where  FT  (Q  )'s  are  total  cost  functions  of  treatment  facilities  in 
t^t 

terms  of  present  value  of  plants  to  be  built  at  time  t  v/ith  design  capa- 
city Q  ,  and  d  's  are  scalars  as  demand  for  wastewater  treatment  at  time 
t,  where  t  =  1975,  1980,  and  1990. 

As  mentioned  above,  FT  (Q. )  for  t  =  1975,  1980,  1990,  includes  con- 
struction costs  as  v:ell  as  O&M  cost.  The  O&M  cost  as  mentioned  in  equa- 
tions (5.2)  and  (5.3)  involves  short-run  cost  functions.   There  are 
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several  points  about  short-run  cost  functions  that  should  be  discussed. 

First,  as  shown  in  Figure  5.3,  L'L'  represents  a  long-run  cost  curve, 

and  f   represents  a  conventional  short-run  cost  curve.   To  the  left  of 
c 

point  g,  short-run  average  cost  increases  as  Q  decreases.   This  may  be 
explained  as  the  indivisibility  of  management,  production  factors,  and 
the  like.   The  decrease  in  Q  will  not  decrease  the  quality  of  effluent. 
However,  to  the  right  of  point  g,  can  one  expand  Q  without  reducing  the 
quality  of  effluent?   Or,  can  one  expand  Q  by  working  overtime,  adding 
more  shifts,  or  increasing  the  inventory  of  raw  sevjage?   The  capacity  of 
a  treatment  plant  is  measured  by  a  rate  (e.g.,  million  gallons  per  day) 
and  the  possibility  of  expanding  treatment  capacity  Q  beyond  the  design 
capacity  without  reducing  the  quality  of  effluent  may  not  exist.   Hence, 
the  segment  of  the  short-run  cost  curve  to  the  right  of  point  g  may  not 
exist  if  the  quality  of  effluent  has  to  be  maintained. 

Second,  there  are  very  fcv;  studies  about  short-run  costs  of  vraste- 
water  treatment,  and  therefore  it  is  difficult,  if  not  impossible,  to 
find  short-run  cost  functions  with  the  capacities  needed  in  this  study. 
As  a  proxy,  long-run  O&M  cost  functions  were  used.   There  are  shortcom- 
ings with  this  approach.   For  example,  if  a  treatment  plant  is  to  be 

build  with  capacity  Q,^„^  in  1975,  one  should  use  f   to  estimate  the  O&M 
^  ^    ^1980  c 

costs  from  1975  through  1980:  however,  there  is  no  information  abour  f  . 

CD  1  ■>  Q 

Therefore,  if  L'L'  is  used  to  estimate  O&M  cost  with  capacity  Q-|QgQ5  it 
is  difficult  to  determine  if  the  cost  has  been  over-  or  underestimated 
without  knowing  f   since  there  is  no  information  about  the  difference 
between  ojgh  and  oacd. 

In  this  study,  each  treatment,  transmission,  or  pumping  facility 
has  been  assumed  to  operate  at  its  design  capacity.   The  reasons  are 


Cost  ($) 
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Capacity  (nigd) 


^1975 


^1980 


Figure  5.3   Opei-ation  and  Maintenance  Cost  over  Time 
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that  1)  there  is  not  enough  information  about  short-run  cost  functions, 
especially  for  advanced  treatment  processes,  2)  this  study  v;as  for  plan- 
ning purposes  and  emphasis  was  placed  on  the  long-run  effect,  and  3)  in 
most  cities  the  differences  in  demands  for  wastewater  treatment  over 
time  are  not  large;  therefore,  the  differences  betv/een  long-run  cost 
estimates  under  design  capacity  and  short-run  cost  estimates  under  capa- 
city in  use  may  not  be  substantial. 

In  this  study,  interest  is  in  satisfying  the  demands  for  wastewater 
treatment  in  1975,  1980,  and  1990.   The  demand  functions  for  wastewater 
treatment  over  time  can  be  considered  as  step  functions  rather  than  as 
continuous  functions.   One  example  is  shown  in  Figure  5.4,  v/here  d(t) 
represents  a  continuous  function,  and  d'  „-,  ,  ^^'igorv^  snd  d'      are 
segments  of  a  step  function.   One  of  the  disadvantages  of  using  step 
functions  for  wastev%fater  treatment  demands  is  that  one  may  underestimate 
the  quantity  of  demand  toward  the  end  of  each  time  period,  but  if  each 
time  period  is  not  large,  this  underestimation  will  not  be  serious.   The 
advantage  of  using  step  functions  is  that  it  simplifies  the  staging 
process.   For  example,  instead  of  considering  Model  2  in  this  study,  we 
consider  Model  1  which  when  compared  to  Model  2  is  easier  to  solve. 
Model  2 

Min        1990 

Z    FT  (Q  ) 
t=1975 

s.t.      ^1975-^^1975 


t" 

t-1975 


E    Q  >  d  „     for  all  t"  <  1990 
^t  -  t"  - 


Q_  >^  0,  1990  <_   t  <  1975 
v;here  FT,  Q,  and  d  are  as  defined  in  Model  1. 
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Capacity  (mgd) 


d(t) 

/ 

'                    i9yo 

^ 

^'1975 

1980 

^\~ 

1975   1980        1990 


—  time  (ye-.ar) 


Figure  5.4    Demand  for  Wastev/ater  Treatment  over  Time 
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The  difference  betv/een  Model  1  and  Model  2   Is  mainly  in  the  time 
variable,  t.   In  Model  2,  t  is  a  continuous  variable;  therefore  the  two 
constraints  in  the  model  are  not  v.'ell  defined.   Unless  one  can  split  the 
time  variable  into  soma  finite  segments,  that  is  to  say  t  is  a  step 
function  rather  than  a  continuous  function,  one  may  not  be  able  to  solve 
Model  2. 

^'^-^   "^1975  ^  ^1975  ^  "^1980'  ^""^  ^^1975^''^'  ^'^1980^''^  ^''  functions  in 
present  value  for  years  1975  and  1980  respectively,  and  x  be  the  design 

capacity.   If  one  can  show  that 


^'^1975^^1975^  "^  ^^1980^*^1980   ^1975^ 


>  ^^'1975(^1980^  ^'-'^ 

then  the  capacities  between  d.„^^  and  d,,^^-  for  year  1975  do  not  have  to 

iy/5       1980 

be  considered.   This  is  to  say,  it  is  more  expensive  to  build  two  treat- 
ment plants  with  capacities  Q,nTE:  and  d     -  Qimi;  ^^  1975  and  1980, 

19/j       19oU     ly/j 

respectively.   If  this  is  true,  the  problem  is  simplified  to  choosing 
the  minimum  cost  plan  from  the  following  options. 


Table  5.1   Reference  Numbers  for  Time  Options  in  Alachua 
County,  Florida,  1975  to  1990 

Demand  for  V/astewater  To  Be  Satisfied 

Optxons                 rn    ^    ^  ■      V         /■    \ 

^                                                Treatment  m  Year(s) 


^1975  ^^" 

^1980  -  ^1975  "^° 

^^1990  -  ^1980  '''' 

^1980  '''' 

^^1990  -  ^1980  "^° 

^1980  '''' 

^1990  -  ^^1980  '''' 
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Table  5,1    (Cont.) 


Demand  for  Wastevzater         To  Be  Satisfied 
°P^^°"^  Treatment  in  Year(s) 


^1975  '''' 

^1990  -  ^1975  '''' 

^^1990  '''' 


Now,  split  FT  into  construction  cost,  CC,  and  Operation  and  mainte- 
nance cost,  OC,  and  consider  construction  cost  only  in  the  following 
inequality: 


CC^975(Ql975)  +  CC^c)gQ(d_^ggQ   ^1975) 


>^^1975^980>  ^'-'^ 

where  d     <  Qm-^n  <  ^moo-   ^^  inequality  (5.6)  holds,  and  if  only  the 

ly/j    ly/j   lybj 

construction  cost  is  considered,  it  is  more  economical  to  build  a  plant 

with  capacity  d,„„„  in  1975  than  any  other  plan  which  can  meet  the  treat- 
1^    -^   1980 

ment  requirement,  d     ,  in  1980.   If  inequality  (5.6)  holds,  one  should 
19oU 

only  consider  the  tradeoff  between  the  savings  in  capital  by  meeting 

d   n   in  1975  with  the  losses  of  O&M  cost  in  excess  capacity.   In  other 
19o(J 

words,  one  only  has  to  consider  three  capacities,  d     ,  d     ,  and 
d   ^_,  in  year  1975,  not  a  continuum  on  the  interval  [d     ,  d     ]. 

iyyu  iy/5   xyyu 

The  following  lemma  is  used  to  prove  that  (5.6)  holds 

0  p 

Lemma   Suppose  f  (x)   =   (1  +  x  )  -  (1  +  x)   is  continuous  on  (O,"')} 

o 

f'(x)  exists  at  some  point  xe(0,«'),  and  0  <  B  <  1,  then  x   >  (1  +  x)  -  1, 

Proof:     f(0)   =   (1  +  0  )  ^  -  (1  -1-  0)^   =   0  =   lim  f(x) 

B-1  B-]  ^^° 

f'(x)   =  f  x^    -  P(l  +  x)^  ^^ 
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Now  (— )     -  (  ■:|-^7 — )   "  >  0,  for  all  x  >  0,  since  x  <  1  +  x.   There- 
fore, -  >  -——  ,    and  since  1  -  6  >  0,  (-)^"^  >  (~—)-^~^.      Hence, 

Xl  +  X  X  1-iX 

f'(x)  >  0,  which  implies  f(x)  is  an  increasing  function  and  lim  f (x)  =  0. 

xiO 
In  other  words, 

f (x)   =   (1  +  x^)  -  (1  +  x)^  >  0        for  X  >  0 

=0  for  X  =  0 

where  0  <  S  <  1  and  x  ^  0.   Therefore,  1  +  x  >  (1  +  x)  ,  or  x  > 

(1  +  x)^  -  1. 

All  the  construction  cost  curves  except  the  one  for  transmission 

pipelines  are  of  form 

CC  (Q  )   =  a  Q^,  Q  >  0,  0  <  B  <  1,  and  te[1975,  1990] 
t^t        t   t  *-     >        L     , 

and  a   >  a  ,  if  t  >  t '  because  the  rate  of  inflation  is  greater  than  the 
rate  of  interest.   If  (5.6)  holds,  the  following  will  hold: 

^1975^1975  "^  "l98o'-'^1980  ~  ^1975'^ 

Divide  both  sides  of  (5.7)  by  o.  to  obtain, 

6  S 

^1975  "^  ^'^"^1980  "  ^1975^  "      ^1980 

,        al980  ^  ^ 
where  y  =  ~~  >  1>  or 

"1975 

^(^^1980  -  ^1975^'  >  ^980  "  ^975  ^'''^ 

where  d     >  Q-,q-,^  >   0,  and  0  <  6  <  1.   Inequality  (5.8)  can  be  re- 
written as: 

,[!i9_80_l]B,  (!l280/.i  (5.9) 

^1975  ^1975 

Therefore,  if  [  ~^^   -  1]"  >  (7:-^)   -  1,  inequality  (5.6)  holds. 
^1975  ^1975 

a       3 

By  lemma,  V7e  have  x  >  (1  +  x)   -  1  if  x  >  0,  and  0  <  6  <  1, 

I  Qfin  1 980 

let  X  =  (      )  -  1  >  0,  then  ^     =  1  +  x.   Hence,  by  the  lemma 

^1975  ^975 
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inequality  (5.9)  holds.   Therefore,  given  the  estimated  cost  functions 
all  one  has  to  do  is  to  choose  the  minimum  cost  plan  from  the  five 
options  mentioned  above. 

Future  Service  Value  of  Treatment  Facilities 

In  this  study,  the  service  lives  for  treatment  plants  and  pumping 
stations  are  assumed  to  be  25  years,  and  for  transmission  pipelines  50 
years.   The  five  options  mentioned  in  the  previous  paragraph  involve 
building  treatment  plants  with  different  capacities  in  different  years. 
In  order  to  make  a  comparison  among  the  costs  of  these  five  options, 
the  cost  of  unused  services  of  each  option  after  1990  should  be  sub- 
tracted from  the  total  cost. 

From  an  economic  point  of  view,  a  capital  asset  is  nothing  but  a 
store  or  reservior  of  valuable  future  services,  from  which  alone  the 
value  of  the  asset  derives.   We  say  "future  services"  because  they  can 
be  used  up  only  over  a  period  of  time.   They  are  deferred  benefits.   If 
these  services  were  valued  at  par,  this  is  to  say,  at  what  they  will 
eventually  be  worth  when  realized,  the  capital  value  of  the  asset  would, 
of  course,  equal  their  sum.   But  a  business  enterprise  does  not  know- 
ingly buy  future  services  at  par  if  profits  are  to  exist.   To  earn  pro- 
fits a  firm  must  buy  the  asset  for  less  than  its  services  will 
eventually  be  worth,  the  profit  over  the  asset  life  as  a  v7hole  being  the 
excess  of  these  service  values  over  the  acquisition  and  use  costs.   The 
capital  value  of  the  asset,  in  other  words,  is  the  discounted  value,  or 
present  worth,  of  these  service  values,  the  latter  in  turn  being  the 
amount  available  for  capital  recovery  and  profit  combined. 


70 


It  is  a  matter  of  co'.--~:on  observation  that  the  services  of  capital 
assets  tend  to  become  less  and  less  valuable  as  time  goes  on.   The  ma- 
jority of  the  assets  require  during  their  service  lives  a  flow  of  main- 
tenance expenditures,  v/hich  as  a  rule  rises  irregularly  with  age  and 
use.   Most  of  them  suffer  a  progressive  deterioration  in  the  quality  or 
the  adequacy  of  their  service.   Moreover,  in  dynamic  technology  they  are 
subject  to  the  competition  of  improved  substitutes,  so  that  the  quality 
of  their  service  declines  relative  to  the  available  alternatives  even 
when  it  does  not  deteriorate  absolutely.    All  of  these  factors — rising 
operation  costs,  impaired  service  quality  or  adequacy,  and  improved 
alternatives — combine  to  reduce  the  value  of  the  service  as  the  asset 
ages. 

In  practice,  the  future  service  value  of  an  asset  has  to  be 
estimated.   Such  estimates  are  quite  common  in  the  valuation  of  real 
estate,  where  annual  service  values  are  projected  and  discounted  into 
capital  value.   It  must  be  admitted  that  they  are  not  much  employed  in 
the  valuation  of  other  types  of  property. 

With  the  exception  of  land,  the  productive  facilities  of  industry 
are  wasting  assets  and  are  depreciated  with  time  and  use  until  the  capa- 
tal  embodied  in  the  facilities  is  exhausted  over  their  productive  ser- 
vice lives.  The  worth  of  the  asset  to  the  owner  at  any  given  point  is 
what  he  could  afford  to  pay  for  it  in  competition  with  the  various  al- 
ternatives then  available,  if  he  did  not  already  have  it.  So  long  as  he 
elects  to  hold  the  asset,  it  may  be  presumed  that  this  worth  is  above 


This  deterioration  in  the  quality  of  service,  both  absolutely  and 
relative  to  current  alternatives,  is  reflected  in  a  general  tendency  to 
reduce  the  intensity  or  continuity  of  use  of  the  asset  as  time  goes  on. 
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resale  value.   It  may  also  be  presumed  that  the  decline  therein  since 
the  asset  was  acquired  has  been  less  than  the  decline  from  original  to 
resale  value. 

Costs  make  up  the  depreciation  base  for  purchased  assets.   In  set- 
ting up  a  depreciation  schedule  a  firm  must  establish  a  useful  life  and 
salvage  value  of  the  asset  in  question.   The  most  common  methods  of  de- 
preciation are  straight-line,  declining-balance,  sum  of  the  digits,  and 
unit-of-production  or  service  [Dougall,  1973,  p.  474].   Since  the  proba- 
ble loss  of  capital  value  is  decidedly  concentrated  in  the  early  part  of 
life,  the  straight-line  v/riteoff  is  not  a  completely  satisfactory  method 
of  depreciation  for  productive  equipment.   Any  realistic  allocation  pro- 
cedure should  get  rid  of  at  least  one  half  of  the  initial  value  over  the 
first  third  of  the  service  life  and  at  least  tv;o-thirds  over  the  first 
half.   The  straight-line  method  is  perhaps  less  objectionable  for  build- 
ings and  structures  than  for  equipment  [Terborgh,  1954,  pp.  37,  47]. 

In  this  study,  wastewater  treatment  facilities  are  long-term  invest- 
ments, and  a  straight-line  method  was  used  to  calculate  the  depreciation 
of  facilities.   The  undepreciated  balance  of  the  assets  was  discounted 
back  to  present  value  for  each  of  the  options  mentioned  above  and  sub- 
tracted from  the  present  value  of  costs. 

Construction  costs,  less  discounted  salvage  value  divided  by  the 
estimated  years  of  life,  produced  (for  each  type  of  fixed  asset)  a  uni- 
form depreciation  e;:pense  each  year.   The  salvage  values  for  treatment 
facilities  in  this  study  were  assumed  negligible.   The  present  value  of 
future  costs  included  both  construction  and  O&M  costs. 
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As  mentioned  in  Chapter  III,  this  study  tries  to  find  the  least 
cost  option  to  satisfy  the  waste-water  treatment  demands  for  the  nine 
cities  in  Alachua  County,  subject  to  five  constraints. 

By  inequality  (5.6),  it  is  obvious  that  if  cost  curves  are  of  form 
CC^(Q^)   =  c.^  Q^   Q  >  0,  0  <  S  <  1,  and  t£[1975,  1990] 
there  is  no  need  to  transport  v/aste  to  more  than  one  location  for 
treatment.   All  the  cost  curves  used  in  this  study  but  two — one  for 
transmission  pipelines  and  one  for  ammonia  stripping — are  of  this  form, 
and  a  study  by  Smith  and  Eilers  [1971]  shows  that  there  are  economies  of 
scale  existing  in  transmission  pipeline  construction  cost,  v/hich  sug- 
gests this  cost  curve  has  the  same  form.   Ammonia  stripping  is  a  very 
small  part  of  the  total  cost  for  the  facilities  studied  here  so  no  v/aste 
was  transported  to  more  than  one  location  for  treatment. 

A  computer  program  v;as  developed  to  scan  through  all  the  possible 

vrays  of  satisfying  wastewater  treatment  demands  of  the  nine  cities  in 

Alachua  County.   This  computer  program  consists  of  one  main  program, 

four  subroutines,  and  two  functions,  they  are  MAIN,  GRAVTY,  TRT,  COST, 

ZERO,  TRI-1,  and  PMG,  respectively  (see  Appendix  I).   Their  individual 

functions  follov;: 

MAIN:   1.   To  read  city  cordinates,  population  figures,  and  calcu- 
late v/astewater  treatment  demands  over  time; 

2.  To  calculate  construction  costs  as  well  as  operation  and 
maintenance  costs  of  pipelines  per  mile,  pumping  sta- 
tions, and  treatment  costs  (both  secondary  and  tertiary 
treatment)  for  each  city  and  time  of  interest  according 
to  the  needs  of  each  city,  and  to  print  out  the  results. 
The  printout  includes  wastewater  treatment  dem.ands, 
population  figures,  and  all  the  needed  cost  figures  over 
time; 

3.  To  generate  a  combination  of  cities  and  send  it  to  sub- 
routine GRAVTY  for  regionalization  analysis. 
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GRAVTY:  1.  Using  nhe  combination  provided  by  MA.IN  to  find  out  the 
potential  site  for  regional  treatment  and  to  calculate 
transmission  costs  for  regional  treatment; 

2,  To  sum  up  the  individual  treatment  costs  of  the  cities 
in  the  combination  and  to  calculate  regional  total  cost; 

3.  To  compare  the  sum  of  Individual  treatment  cost  and 
regional  total  cost.   If  regional  total  cost  is  lower, 
print  out  the  cordinates  of  potential  treatment  site, 
regional  treatment  costs,  transmission  costs,  pumping 
station  costs,  sum  of  individual  treatm.ent  costs,  and 
the  distance  from  each  city  in  the  combination  to  the 
potential  site.   If  the  sum  of  individual  costs  is 
lower,  return  to  step  3  in  MAIN  to  get  a  new  combination, 
until  all  the  possible  combinations  are  exhausted. 

TRT,  TRI-I,  PMG:   To  provide  cost  figures  of  treatment,  transmission  pipe- 
line, and  pumping  station,  respectively,  which  are  re- 
quired by  MAIN,  GRAVTY,  and  COST. 

COST:   Using  the  combination  provided  by  MAIN  to  screen  out  the 
required  information  about  costs  and  send  it  to  GRAVTY 
for  further  calculation. 

ZERO:   To  initialize  some  of  the  variables  in  COST. 
There  are  491  potential  sites  for  each  time  option,  five  time  op- 
tions, two  kinds  of  secondary  treatment,  and  three  tertiary  treatment 
processes.   The  program  is  instructed  to  choose  one  secondary  treatment 
with  and  without  tertiary  treatm.ent  processes,  one  time  option,  and  one 
potential  site  at  a  time  to  calculate  cost  figures,  and  give  answers  to 
the  feasibility  of  regional  treatment  until  all  the  possible  combinations 
of  treatment  types,  time  options,  and  potential  sites  are  exhausted. 
This  is  a  brute-force  search  for  the  solution  of  the  optimization  model 
mentioned  in  Chapter  III. 

Minimum  Cost  Option  for  Wastewater  Treatment 

In  the  last  chapter  it  was  shown  that  if  only  secondary  treatment 
facilities  are  required  a  regional  wastewater  treatment  plant  is  not 
needed.   In  this  case  individual  high-rate  trickling  treatment  plants 
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v;il]  give  minimuni  cost  option  provided  that  d,^-,^,  d-,^„^,  and  d,,,_„  are 
satisfied  in  years  1975,  1980,  and  1990,  respectively.   However,  under 
the  same  scheme,  if  tertiary  treatment  processes  are  required,  there  is 
one  regional  treatm.ent  plant  site  vhich  has  a  higher  gain  in  economies 
of  scale  than  the  losses  in  pumping  and  transmission  pipelines,  and  has 
a  minimum  cost  among  the  nine  options  under  considerations,  and  this 
regional  treatment  plant  has  a  capacity  which  satisfies  the  1975  demands 
for  wastev;ater  treatment  of  three  cities  in  1975. 

The  same  procedures  in  the  previous  chapter  were  used  to  scan 
through  all  the  combinations  of  tim.e,  treatment  types,  and  combinations 
of  cities  to  meet  the  demand  for  each  of  the  time  options  mentioned  in 
Table  5.1. 

To  simplify  the  discussion,  Table  5.2  shows  all  the  possible  re- 
gionalization  options,  and  a  num.ber  is  assigned  to  each.   In  the  follow- 
ing discussion  or  representation,  instead  of  repeating  the  constituents 
of  these  options,  reference  numbers  v/ill  be  used.   For  example,  regional 
option  11  represents  the  plan  which  consists  of  a  regional  treatment 
plant  which  serves   cities  1  through  8  and  a  single  plant  serving  city  9. 

Table  5.3  shows  the  seven  ways  of  satisfying  demands  over  1975, 

1980,  and  1990.   Present  value  of  future  cost  of  possible  regional 

treatment  plants,  cities  to  be  served,  and  capacities  to  be  built  by 

types  are  also  presented.   For  example,  if  d   ^   is  satisfied  in  1975, 

19  /-' 

there  are  four  possible  regional  high-rate  trickling  filter  regional 
plants  with  tertiary  processes,  seven  possible  regional  activated  sludge 
plants  with  tertiary  processes,  and  txv^o  possible  regional  activated 
sludge  plants  v/ithout  tertiary  treatment  processes.   Present  value  of 
costs  up  to  1990  in  terms  of  September  1974  dollars  are  shov^/n  in  columns 
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Table  5.2   Reference  Numbers  for  Regional  Wastewater  Treatment  Options 
in  Alachua  County,  Florida,  1975  to  1990 

Regional  Wastewater  Treatment  Options         Reference  Number 


9 

Z 

1=1 

R(i) 

R(l) 

+ 

R(2, 

.3) 

9 
+      Z 
1=4 

R(l) 

R(l, 

2, 

,    3) 

+ 

9 

Z      R(: 
L=4 

L) 

R(l) 

+ 

R(2. 

,    3 

,    4)   + 

9 

Z      R(j 
1=5 

0 

R(l, 

2, 

.    3, 

A) 

9 
+      Z 
1=5 

R(i) 

7 
R(l,  2)  +  Z   R(i)  +  R(8,  9)  6 

1=3 

7 
R(l)  +  R(2,  3)  +  Z   R(i)  +  R(8,  9)  7 

1=4 

7 
R(l,  2,3)+   Z   R(l)  +  R(8,  9)  8 

i=4 

7 
R(l)  +  R(2,  3,4)+   Z   R(i)  +  R(8,  9)  9 

1=5 

7 

R(l,  2,  3,  4)  +  Z  R(i)  +  R(8,  9)  10 

1=5 

R(l,  2,  3,  4,  5,  6,  7,  8)  +  R(9)  11 
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Table  5.3    Present  Value  of  Costs  of  Alternative  Regional  Treatment^ 
Plants  for  Alachua  County,  Florida,  1975,  1980,  and  1990^ 


Activated  Sludge     High-rate  Trickling 
System  Filter  System 


Regional  Treatment 


Plant  tertiary  processes     tertiary  processes 


with      without     with      v;ithout 


1.  d      is  satisfied  in  1975 

R(2,  3)  2,232,163     2,147,296 

R(l,  2)  3,603,283 

R(8,  9)  2,216,676 

R(l,  2,  3)  4,336,789     4,086,862 

R(2,  3,  4)  4,550,842     4,261,926 

R(l,  2,  3,  4)  6,171,936   3,422,759   5,816,155 

R(l,  2,  3,  4,  5,  6,  7,  8)   35,550,720  9,605,506 

2.  d^gg^  -  d^g^^  is  satisfied  in  1980 

R(l,  2,  3)  758,340 

R(2,  3,  4)  657,067 

R(l,  2,  3,  4)  946,246     880,555 

^-    ^1990  "  ^1980  ^^  satisfied  in  1990 

R(l,    2,    3,    4)  226,315            —                 235,555 

""■    '^1990  ~  ^1980  ^^  satisfied  in  1980 

R(2,  3)  618,276 

R(l,  2,  3)  1,260,849     —     1,273,937 

R(2,  3,  4)  1,426,352 

R(l,  2,  3,  4)  1,859,144     1,742,912 

5-  d^-o_  is  satisfied  in  1975 
1980 

R(l,  2)  3,824,863 

R(2,  3)  2,495,646     2,301,272 

R(8,  9)  2,297,853 

R(l,  2,  3)  4,650,832     4,357,962 

R(2,  3,  4)  4,773,810     —     4,474,878 

R(l,  2,  3,  4)  6,506,944   3,558,569   6,137,994 

R(l,  2,  3,  4,  5,  6,  7,  8)   40,885,104 


Table  5.3    (Cont.) 
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Regional  Treatment 
Plant 


Activated  Sludge 
System 


tertiary  processes 


with 


without 


High-rate  Trickling 
Filter  System 


tertiary  processes 


with 


without 


6.  d  is  satisfied  in  1975 

R(l,  2) 

R(2,  3) 

R(8,  9) 

R(l,  2,  3) 

R(2,  3,  4) 

R(l,  2,  3,  4) 


4,327,626 

2,830,965     2,617,682 

2,478,977 

5,272,472     -—  4,952,246 

5,516,895     — -  5,185,699 

7,474,954  4,052,585  7,068,858 


R(l,  2,  3,  4,  5,  6,  7,  8)   53,032,960 


"^1990  ''l')75 

R(l,  2) 
R(2,  3) 

R(l,  2,  3) 

R(2,  3,  4) 

R(l,  2,  3,  4) 


is  satisfied  in  1980 

1,256,987 
872,681 
1,594,391 
1,681,106 
2,244,183 


1,489,681 
1,571,759 

2,110,057 


All  costs  are  in  September  1974  dollars, 


78 


2  to  5  in  Table  5.3.   Other  possibilities  and  their  costs  up  to  year 
1990  are  also  shov/n  in  Table  5.3. 

In  order  to  find  the  least  cost  combination  these  cost  figures 
(Table  5.3)  v;ere  assembled  V7ith  other  individual  treatment  costs  to 
estimate  county  total  costs.   All  the  regional  plants  shown  in  Table 
5.3  have  lower  costs  than  the  sun  of  the  costs  of  individual  plants  for 
the  cities  each  regional  plant  serves.   Therefore,  for  the  five  time 
options  mentioned  in  Table  5.1,  regional  treatment  plants  were  chosen 
first,  and  the  discrepancies  over  time  and  over  space  were  filled  by  in- 
dividual plants. 

Table  5.4  shows  all  the  county  total  cost  figures  according  to  time 
options  and  regional  wastewater  treatment  schemes.   Cost  figures  are 
discounted  present  values  of  all  costs  up  to  1990.   From  these  figures, 
one  finds  that  time  option  2,  i.e.,  the  option  to  satisfy  d     in  1975 
and  d     -  d  qq.^  in  1990,  has  the  lowest  cost  for  both  systems  and  for 
all  the  treatment  schemes.   And  as  before,  the  activated  sludge  system 
has  a  higher  cost  than  the  high-rate  trickling  filter  system. 

Though  simpler  and  cheaper  to  operate,  the  typical  trickling  filter 
plant  is  being  used  less  and  less  on  domestic  sewage  in  North  America 
due  to  the  fact  that  it  does  not  achieve  the  percent  removal  of  organic 
material  that  an  activated  sludge  plant  does.   The  least  cost  method 
thus  depends  on  the  required  quality  of  the  effluent.   Typical  effluents 
from  trickling  filter  plants  treating  domestic  wastes  have  BOD's  and  SS 
usually  greater  than  20  mg  per  liter.   Typical  removal  efficiencies  for 
both  BOD  and  SS  are  in  the  area  of  80%.   Activated  sludge  plant  efflu- 
ents have  BOD's  and  SS  between  10  and  20  mg  per  liter  and  removal  ef- 
ficiencies in  the  area  of  90%. 
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Hence,  if  the  required  quality  of  secondary  effluent  is  set  at  90% 
removal,  a  trickling  plant  will  not  achieve  this  level  of  removal,  and 
an  activated  sludge  plant  would  be  required  to  provide  the  least  cost 
method  of  treating  wastewater.   The  results  in  Table  5.4  show  that  if 
only  secondary  treatment  is  required,  the  least  cost  method  suggests 
that  the  cities  of  High  Springs,  Alachua,  Archer,  and  Newberry  should 
cooperate  to  build  a  regional  treatment  plant  v/ith  cordinate  x=5.02, 
y=18.77  in  1975  to  satisfy  their  1980  demands.   Added  demands  in  1990 
were  satisfied  in  the  solution  by  the  construction  of  individual  plants. 
For  the  other  five  cities,  tim.e  option  2  vjas  used  to  satisfy  their  de- 
mands by  building  individual  activated  sludge  plants  for  each.   Hov/ever, 
if  tertiary  treatment  processes  are  required,  one  V70uld  build  two 
regional  activated  sludge  plants,  one  for  the  cities  of  High  Springs, 
Alachua,  Archer,  and  Newberry  with  cordinate  x=5.02,  y=18.77,  and  one 
for  Kaldo  and  Hawthorne  w^ith  cordinate  x=32.57,  y=6.62,  in  1975  to  sat- 
isfy the  1980  demands  of  the  six  cities.   Added  demands  in  1990  v/ere 
satisfied  in  the  solution  by  the  construction  of  individual  plants.   For 
the  other  cities — Gainesville  Urban  Area,  LaCrosse,  and  Micanopy--time 
option  2  and  individual  activated  sludge  plants  were  used  to  satisfy 
demands . 

If  the  required  quality  of  secondary  effluent  is  set  at  80%  re- 
moval, a  trickling  filter  plant  vrill  provide  the  least  cost  method 
(Table  5.4)  of  treating  wastewater.   The  solution  suggests  that  each 
city  should  built  its  own  high-rate  trickling  filter  plant  according  to 
time  option  2.   Hovjever,  if  tertiary  treatment  processes  are  required, 
one  would  build  a  regional  high-rate  trickling  filter  plant  for  the 
cities  of  Alachua,  Archer,  and  Nev.berry  v.'ith  cordinate  x=4.02,  y=12.03 
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in  1975  to  satisfy  the  19S0  demands  of  the  three  cities.   Added  demands 
for  these  three  cities  in  1990  were  satisfied  in  the  solution  by  the 
construction  of  individual  plants.   For  the  other  six  cities — High 
Springs,  the  Gainesville  Urban  Area,  LaCrosse,  Micanopy,  Waldo,  and 
Hawthorne — time  option  2  and  high-rate  trickling  filter  plants  should 
be  used  to  satisfy  the  demands. 

Figure  5.5  shows  the  locations  of  the  regional  activated  sludge 
treatment  plants  and  the  lengths  of  pipelines  from  the  cities  to  these 
two  sites.   Figure  5.6  does  the  same  for  the  regional  trickling  filter 
treatment  plant.   Tables  5.5  and  5.6  show  the  details  of  the  minimum  cost 
schemes  for  secondary  treatment  and  advanced  treatment,  respectively. 
The  figures  in  parentheses  are  the  corresponding  treatment  capacities  in 
mgd  to  be  built.   The  difference  in  treatment  capacities  between  second- 
ary treatment  and  tertiary  treatment  in  the  Gainesville  Urban  Area  is 
due  to  the  fact  that  the  Gainesville  Urban  Area  has  secondary  treatment 
plants  with  capacity  9.5  mdg,  but  has  no  tertiary  treatment  facilities. 
Therefore,  if  tertiary  treatment  is  required,  the  facilities  should  be 
large  enough  to  handle  all  the  wastewater  generated  in  this  area. 

The  results  shown  in  Tables  5.5  and  5.6  were  calculated  under  the 
assumption  that  one  can  build  a  treatment  plant  as  small  as  one  likes. 
However,  this  assumption  may  not  be  realistic.   If  one  suspects  that  the 
capacity,  for  instance,  required  for  Micanopy  in  1990  under  time  option 
2  is  too  small  to  be  practical,  then  time  option  5  would  provide  second 
best  solutions.   The  present  values  of  cost  for  each  city  or  regional 
plant  for  both  secondary  treatment  and  secondary  treatment  v.'ith  tertiary 
processes  for  time  option  5  are  presented  in  Table  5.7.   The  conclusions 
are  still  the  same  as  those  from,  time  option  2,  except  the  cost  figures 
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Table  5.5    Present  Value  of  Secondary  Treatment  Costs  for  the  Minimum 
Cost  Option  for  Alachua  County,  Florida,  1975  to  1990^ 


Activated  SI 

.udge 

High-rate 

Trickling 

City 

System 

Filter 

System 

1975 

1990 

1975 

1990 

High  Springs 



39,059 
(.0871) 

969,928 
(.4455) 

33,919 
(.0871) 

Newberry 



23,559 
(.0385) 

583,362 
(.2011) 

20,183 
(.0385) 

Archer 



21,182 
(.0324) 

483,274 
(.1498) 

18,094 
(.0324) 

Alachua 



48,008 
(.1215) 

1,016,294 
(.4792) 

41,920 
(.1215) 

Gainesville 
Urban  Area 

6,204,332 

(6.5446) 

697,737 
(9.0876) 

5,410,354 
(6.5446) 

652,094 
(9.0876) 

LaCrosse 

366,428 
(.0655) 

13,438 
(.0155) 

284,666 
(.0655) 

11,337 
(.0155) 

Micanopy 

520,765 
(.1161) 

6,441 
(.0047) 

410,530 
(.1161) 

5,324 
(.0074) 

Waldo 

533,672 
(.1208 

12,325 
(.0135) 

421,134 
(.1208) 

10,373 
(.0135) 

Hawthorne 

675,850 
(.1775) 

16,799 
(.0223) 

538,577 
(.1775) 

14,260 
(.0223) 

R(l,  2,  3,  4) 

3,588,569^ 
(1.2756) 







All  costs  are  in  September  1974  dollars. 

This  figure  includes  pumping  and  transmission  costs, 
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Table  5.6    Present  Value  of  Secondary  and  Tertiary  Treatment  Costs  for 
the  Minimum  Cost  Option  for  Alachua  County,  Florida,  1975 
to  1990^"^ 


Activated  Sludge 

High-rate 

Trickling 

City 

System 

Filter 

System 

1975 

1990 

1975 

1990 

High  Springs 

___ 

73,790 

2,248,431 

68,651 

(.0871) 

(.4455) 

(.0871) 

Newberry 



42,890 
(.0385) 



39,514 
(.0385) 

Archer 



38,313 
(.0324) 



35,225 
(.0324) 

Alachua 



92,373 
(.1215) 



86,285 
(.1215) 

Gainesville 

31,090,496 

1,535,004 

30,296,512 

1,489,361 

Urban  Area 

(16.0446) 

(4.5876) 

(16.0446) 

(4.5876) 

LaCrosse 

673,327 

23,738 

591,564 

21,637 

(.0655) 

(.0155) 

(.0655) 

(.0155) 

Micanopy 

986,469 

11,089 

876,234 

9,972 

(.1161) 

(.0047) 

(.1161) 

(.0047) 

Waldo 



21,690 

900,693 

19,739 

(0.35) 

(.1208) 

(.0135) 

Hawthorne 



29,997 

1,176,076 

27,458 

(.0223) 

(.1775) 

(.0223) 

R(8,  9) 

2,297,853^ 
(.2983) 

— 

— 



R(2,  3,  4) 





4,474,878^ 
(.8301) 

— 

R(l,  2,  3,  4) 

6,506,944^ 
(1.2756) 







All  costs  are  in  September  1974  dollars. 


This  figure  includes  pumping  and  transmission  costs. 
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Table  5.7   Present  Value  of  Wastewater  Treatment  Costs  for  Time 
Option  5  for  Alachua  County,  Florida,  1975  to  1990^ 


Activated  Sludge 

System 


High-rate  Trickling 
Filter  System 


City 


High  Springs 
Newberry 
Archer 
Alachua 


Gainesville 
Urban  Area 

LaCrosse 


Micanopy 
Waldo 
Hawthorne 
R(3,  9) 
R(2,  3,  4) 
R(l,  2,  3,  4) 


tertiary  processes 


with 


41,808,368 
(20.6322) 

775,558 
(.0810) 

1,013,231 
(.1208) 


2,478,977 
(.3341) 


7,474,954 
(1.5551) 


without 


10,602,925 
(15.6322) 

417,534 
(.0810) 

533,672 
(.1208) 

569,520 
(.1343) 

726,721 
(.1998) 


tertiary  processes 


with 

2,554,474 
(.5326) 


40,650,256 
(20.6322) 

684,164 
(.0810) 

900,693 
(.1208) 

969,109 
(.1343) 

1,277,120 
(.1998) 


5,185,699 
(1.0225) 


v;ithout 


1,087,224 
(.5326) 

652,435 
(.2396) 

547,699 
(.1822) 

1,174,284 
(.6007) 

9,444,815 
(15.6322) 

326,139 
(.0810) 

421,134 
(.1208) 

450.637 
(.1343) 

580,857 
(.1998) 


4,052,585 
(1.5551) 


All  costs  are  in  September  1974  dollars. 


b^,  . 


This    figure    Includes   puraping  and   transmission   costs. 
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are  different.   Again  the  figures  in  parentheses  under  each  set  of  cost 
figures  are  the  capacities  to  be  built. 

North  Central  Florida  Regional  Planning  Council  [1973]  has  six  al- 
ternatives toward  supplying  wastewater  treatment  facilities  in  Alachua 
County  and  these  alternatives  are  shown  in  Table  5.8.   There  are  differ- 
ences between  their  estimated  capacities  and  those  used  in  this  study. 
There  are  also  differences  in  basic  assumptions  used  in  consolidations 
over  space  and  time,  e.g.,  highway  milages  were  used  in  NCFRPC's  study 
instead  of  a  Euclidean  distance  measurement  which  was  used  in  this 
study.   Therefore,  a  meaningful  comparison  betvjeen  their  alternatives 
and  the  results  of  this  study  is  difficult. 
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Table  5.8   Alternative  Plans  of  Regional  Wastevjater  Utilities  in 
Alachua  County  as  Suggested  by  NCFRPC 

Alternative  Regional  Plants 

9 
I         I      R(i) 
i=l 

II        R(5,  7,  8,  9)  +  R(l,  2,  3,  4,  6) 

III        R(l,  2,  3,  A,  6)  +  R(8,  9)  +  R(5)  +  R(7) 

9 

IV        R(l,  4,  6)  +  R(2,  3)  +  R(5)  +  S   R(i) 

i=7 

9 

V        R(l,  2,  4)  +  R(3)  +  Z      R(i) 

i=5 

b     ^ 
VI        R(l,  4,  5,   6)  +  R(2)  +  R(3)  +  R(5  )  +  Z   R(l) 

i=7 


This  regional  treatment  plant  serves  only  part  of  the  NW  section 
of  Gainesville. 

This  regional  treatment  plant  serves  Gainesville  except  for  the 
part  of  NTJ  Gainesville. 

Source:   North  Central  Florida  Planning  Council  [1973]. 


CHAPTER  VI 

LANDS  PR  EiVDING  SECONDARY  EFFLUENT  AS  A 

SUBSTITUTE  FOR  TERTIARY  TREATMENT  PROCESSES 


The  technology  employed  in  v/astewater  treatment  has  not  advanced 
much  in  the  past  several  decades.   There  are  several  difficulties  one 
encounters  vjhen  attempting  to  remove  more  than  95%  of  the  BOD  and  SS 
with  the  standard  processes.   Greater  removals  require  very  large  in- 
creases in  detention  tdjnes  and  a  corresponding  increase  in  tank  sizes. 
A  second  treatment  problem  is  the  removal  of  nutrients.   Nitrogen  and 
phosphorous  compounds  are  more  likely  to  be  responsible  for  excess  weed 
growth  than  BOD.   Standard  processes  do  not  do  a  good  job  of  removing 
these  nutrients.   These  tv/o  problems  can  be  handled  by  tertiary  treat- 
ment processes;  however,  the  cost  for  these  treatment  processes  alone 
is  higher  than  the  cost  of  conventional  treatment  processes.   Land 
treatment  is  therefore  considered  as  an  alternative  to  tertiary  vjaste 
treatment. 

Briefly  stated,  land  treatment  involves  the  use  of  agricultural 
land  and  crops  or  forest  products  to  absorb  and  filter  nitrates,  phos- 
phates, and  other  elements  from  wastewater  that  has  undergone  primary 
and,  usually,  secondary  treatment.   Excess  "purified"  water  is  then 
returned  to  the  water  course.   The  methods  of  applying  wastewater  to 
the  land  can  be  identified  as  infiltration  systems,  crop  irrigation 
systems,  and  spray-runoff  systems  [Thomas,  1973;  Thomas  and  Lav?,  1963]. 
Infiltration  systems  are  usually  designed  to  prevent  surface  runoff. 
High  loading  rates  make  evaporative  losses  relatively  insignificant. 
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and  up  to  99%  of  the  applied  v/astewater  may  be  contributed  to  ground 
water  as  recharged  [Lavcrty  et  al. ,  1961].   Crop  irrigation  systems  may 
or  may  not  control  surface  runoff.   Lov;  loading  rates  allow  much  of  the 
applied  V7astewater  to  be  lost  through  evapotranspirat ion,  and  the  con- 
tribution to  ground  water  is  largely  dependent  on  evapotranspiration 
losses  [Dalton  and  Murphy,  1973;  Graveland  and  Vickerman,  1972;  Parizek 
et  al.  ,  1967;  Sprout  and  Hopkins,  1972;  Young  ct  al. ,  1972].   Spray 
runoff  systems  are  designed  to  return  50Z  or  more  of  the  applied  waste- 
water as  direct  surface  runoff,  evapotranspiration  losses  are  variable 
but  relevant,  and  the  selection  of  sites  with  impermeable  soils  re- 
stricts the  contribution  to  ground  v/ater  [Law  et  al. ,  1970].   The  appli- 
cation of  wastewater  to  land  serves  the  following  purposes:   promotes 
growth  of  crops,  conserves  vrater  and  nutrients  that  are  normally  wasted, 
provides  economical  treatm.ent  of  the  wastewater,  and  reduces  the  pollu- 
tion load  on  surface  water  supplies. 

Many  researchers  have  contributed  estimates  of  the  fertilizer  value 
in  treated  sewage  effluent,  and  the  values  vary  v/ith  the  study,  with  the 
type  of  treatment  given  to  the  sewage  effluent,  and  with  the  source  of 
effluent.   Schreiber  [1957]  reports  that  the  amount  of  the  principal 
fertilizer  elements  in  sewage  effluents  from  15  California  cities  stud- 
ied is  60  to  100  pounds  of  nitrogen  chiefly  in  the  form  of  ammonia;  20 
to  40  pounds  of  potassium  occuring  as  potassium;  and  60  to  100  pounds  of 
phosphate  occuring  as  phosphate,  per  acre-foot  of  effluent.   Schreiber 
does  not  assign  a  dollar  value  to  the  fertilizer  contained  in  an  acre- 
foot  of  effluent;  however,  Hirsch  [1959]  states  that  the  1969  value  of 
nitrogen,  phosphorus  and  potassium  fertilizer  in  an  acre-foot  of  typical 
San  Diego  wastewater  is  $18.   A  Pennsylvania  State  University  research 
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expariment  [Parizek  et  al. ,  1967]  attempted  to  relate  effluent-derived 
agricultural  nutrients  to  an  equivalent  commercial  fertilization  rate. 
In  this  case  the  investigators  determined  that  with  an  application  rate 
of  two  inches  of  effluent  per  week  for  periods  varying  from  24  to  33 
weeks  of  irrigation  in  three  years  of  study  (1963-1965) ,  fertilizer  ap- 
plication to  crops  and  wooded  plots  was  equal  to  100  to  280  pounds  of 
nitrogen,  180  to  300  pounds  of  available  phosphorus,  and  120  to  320 
pounds  of  v;ater  soluable  potash  per  acre. 

From  research  on  the  effect  of  sev;age  farming  on  crop  yields,  it 
is  evident  that  most  crops  produce  much  higher  yields  when  irrigated 
with  effluent  than  when  not  irrigated  or  when  irrigated  with  ordinary 
water  but  not  commercially  fertilized  [Parizek  et  al. ,  1967;  Day  and 
Tucker,  1959]. 

The  suitability  of  domestic  V7aste^^7ater  for  irrigation  depends  on 
the  chemical  composition  of  water  supplies  and  how  this  composition  is 
affected  by  the  use  cycle.   For  instance,  salinity  and  the  ratio  of 
sodium  to  calcium  and  magnesium  are  important  factors  in  judging  suita- 
bility of  water  from  saline  soil  solutions,  the  osmotic  pressure  of 
such  solutions  interferes  with  the  movement  of  water  from  the  solution 
into  the  plant  root,  and  the  plants  may  suffer  from  incipient  drought. 
Sodium  soils  are  relatively  impermeable  to  air  and  water;  they  are  hard 
when  dry,  difficult  to  till,  and  plastic  and  sticky  when  wet.   These 
adverse  physical  conditions  retard  or  prevent  germination  and  water  re- 
moval by  plants  and  are  generally  unfavorable  for  plant  growth. 

Most  of  the  public  reservations  concerning  the  use  of  treated 
sewage  effluent  for  irrigation  and  fertilization  stem  from  the  feeling 
that  the  effluent  is  unclean,  unaesthetic,  and  contaminated  with 
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disease  organisms.   Because  of  this  attitude,  the  current  general 
practice  in  the  United  States  is  to  provide  at  least  primary  treatment 
and  to  utilize  treated  sewage  for  irrigation  of  nonfood  crops  only. 
Secondary  treatment  and  chlorination  greatly  improved  the  quality  of 
sewage,  and  for  aesthetic  safety  reasons  these  conventional  wastewater 
treatment  processes  are  recommended  for  upgrading  municipal  wastewater 
before  irrigation  use. 

Possibly  the  most  significant  standard  for  irrigation  water,  from 
the  standpoint  of  public  interest,  is  that  for  bacterial  count.   Re- 
strictions on  coliform  counts,  the  primary  index  used  for  bacterial  con- 
tamination, are  relatively  rigorous  for  irrigation  v;ater.   But  it  is 
important  to  realize  that  the  presence  of  coliform  bacteria  in  \i7aste- 
water  is  only  an  indication  of  the  potential  presence  of  pathogenicity, 
as  coliform  bacteria  themselves  are  not  pathogenic. 

Another  most  pressing,  most  difficult  and  expensive  unknoxsfn  to  be 
determined  is  the  fate  of  viruses  present  in  effluent  used  for  spray 
irrigation.   Even  the  most  efficient  secondary  treatment  plants  fail  to 
produce  virus-free  effluents  [Kollins,  1956].   The  field  studies  at 
Stanford  University,  Hawaii,  and  MIT  [Drcwry,  1973]  have  negative  re- 
sults on  the  tests  for  virus  existance  in  the  test  wells.   However, 
Wellings  et  al.  [1974]  report  that  they  have  isolated  virus  from  the 
effluent  used  for  spray  irrigation  from  the  spray  just  before  it  enters 
the  soil,  and  from  veir  water  which  represents  effluents  percolated 
through  20  feed  of  sandy  soil  at  St.  Petersburg  Northwest  Wastev/ater 
Treatment  Plant,  Florida. 

Farmers  ordinarily  irrigate  the  best  soils  or  fields  on  their 
farms.   The  same  should  be  true  for  municipalities  when  irrigating 


93 


with  effluent.   Unfortunately,  too  many  make  the  sad  choice  of  acquir- 
ing the  cl'ieapest  land  available  which  often  is  too  wet,  too  steep,  or 
has  insufficient  soil  depth  to  provide  adequate  renovation. 

Nev;  concepts  are  required  v?hon  managing  irrigation  systems  for 
wastewater  disposal.   Whereas  farmers  ordinarily  irrigate  for  crop  re- 
quirements, in  wastewater  irrigation  the  need  is  to  manage  for  adequate 
water  renovation.   The  crop  is  secondary.   This  suggests  that  usually 
the  municipality  or  company  producing  the  wastewater  should  own  the 
land  and  not  the  farmer.   A  farmer  whose  income  is  from  the  sale  of 
crops  will  find  it  difficult  to  irrigate  a  crop  which  perhaps  already 
has  more  V7ater  than  it  needs.   In  v;astewater  disposal  management,  one 
normally  applies  maximum  amounts  of  water  to  land,  not  optimal  amounts 
for  crops.   Most  of  the  irrigation  research  and  experience  to  date  has 
been  oriented  toward  crop  objectives  and  not  wastewater  disposal. 

Monitoring  of  the  effluent  which  leaves  the  renovation  area  is  an 
essential  part  of  system  management.   Both  surface  and  subsurface  water 
discharge  should  be  monitored  in  each  dominate  direction  of  water  egress. 
Monitoring  facilities  usually  assay  only  the  quality  of  the  discharge 
from  the  renovation  area.   Under  proper  management,  however,  quantity 
also  should  be  considered. 

In  this  study,  ive  assume  that  secondary  treatment  and  chlorination 
are  required  for  the  effluent  to  be  used  in  land  treatment.   According 
to  the  results  in  Chapter  V,  if  only  the  secondary  treatment  is  required, 
there  is  no  need  to  build  any  regional  treatment  plant,  and  each  city 
can  satisfy  its  demands  from  1975  to  1990  by  the  second  option  mentioned 
in  that  chapter.   If  land  treatment  is  used  as  a  substitute  for  tertiary 
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treatment  and  added  to  the  secondary  treatment  plants,  the  feasibility 
of  regional  land  treatment  operations  comes  into  question. 

The  answer  to  this  question  will  depend  on  the  economies  of  scale 
in  land  treatm.ent  and  whether  the  gains  in  economies  of  scale  are  large 
enough  to  offset  the  increase  in  the  costs  of  transmitting  effluent  to 
regional  land  treatment  sites. 

Data  do  not  exist  to  investigate  the  economies  of  scale  in  land 
treatment  operations,  nor  are  there  definite  studies  that  can  be  used 
for  reference.   However,  there  are  studies  [Wills,  1955;  Madden,  1967] 
on  farm  size  which  suggest  that  economies  of  size  in  agricultural  pro- 
duction are  exhausted  by  family  size  farm  units,  and  that  cost  func- 
tions arc  homogeneous  of  degree  one  for  modest  size  units.   One  would 
expect  a  simijar  situation  for  land  treatment  operations. 

Meyers  [1973]  points  out  that  in  the  design  of  an  irrigation  system 
for  any  xcastewater  facility,  five  variables  must  be  considered.   These 
are  dndly  flow  from  the  facility,  v;eekly  loading  depth  on  the  land  area, 
hourly  application  rate,  sprinkler  spacing,  and  nozzle  operating 
pressure.   These  variables  are  interrelated. 

Given  these  five  variables  one  can  attempt  to  develop  an  engineer- 
ing design  of  the  land  spreading  system.   Design  parameters  to  develop 
such  a  system  have  not  been  developed  for  Alachua  County.   And,  rather 
than  construct  a  hypothetical  system  v/ith  data  that  have  little  or  no 
relevance  for  the  design  problem  here  in  Alachua  County,  this  study  pre- 
sents a  cost-minimization  model  which  will  determine  the  desired  balance 
among  the  multiple  objectives  of  economic  crop  production,  minim.um  cost 
waste  disposal,  and  environmental  quality  objectives. 
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In  general  form  the  model  may  be  presented  as  a  set  of  eight 
functions  as  follows: 

Y  =   y(E  ,  L,  X,  F)    crop  yield  response 

function  (6.1) 

F   =   f(E  )  nutrient  input  function      (6-2) 

A   =   a(E  )  land  requirement  function    (6.3) 

TC    =   c  (E  )  cost  of  disposing  of  E   on 

crop  land  (6.4) 

TC    =   c  (E  )  cost  of  treating  and  dis- 

posing of  E  (6.5) 

D  ■-=  d  (E  )  damage  function  for  E  (6.6) 

D  =  d  (E  )  damage  function  for  E  (6.7) 

C  =  c  (L)  cost  function  for  land  (6.8) 

where: 

E   =  quantity  of  secondary  effluent  to  be  applies  to 
land  in  crops 

E„   =   quantity  of  secondary  effluent  for  disposal  by  an 
alternative  method,  say,  tertiary  treatment 

E   =   total  quantity  of  secondary  effluent  for  disposal 

F   =   total  quantity  of  nutrients  in  E 

TC    -   cost  of  transporting  and  applying  waste,  E   to 
cropland 

I'C    =   cost  of  treating  and  disposing  of  secondary 
effluent,  E  ,  by  tertiary  treatment 

L   =   acres  of  cropland  used  for  waste  disposal 

X  =   quantity  of  resources  (other  than  E   and  L)  used 
in  crop  production 

Y  =   physical  quantity  of  crop  production 

C    =   cost  per  acre  to  gain  control  of  land  (maybe  rent, 
annual  cost  of  purchase,  etc.) 

C   =   cost  per  unit  of  non-land  and  non-v7aste  resources 

X 

used  in  crop  production 
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P   =   selling  price  per  unit  of  crop 

y 

D   =   dollars  of  damage  (considered  as  a  cost  to  the 

system)  to  environmental  community  caused  by  use 
of  E   on  crop  land 

D   =   dollars  of  damage  (considered  as  a  cost  of  the 

system)  to  environmental  community  caused  by  dis- 
charge tertiary  effluent  of  Ey    (alternative  to  the 
disposal-on-land  method  of  treatment  and  disposal) 

It  is  assumed  in  the  model  that  the  decision  unit  wishes  to  m.ini- 
mize  total  cost  of  disposing  of  a  given  amount  of  secondary  effluent  E 
after  allowing  for  damage  to  the  environment  and  community  and  any  off- 
setting return  from  crops.   The  decision  unit  may  purchase  land  di- 
rectly or  rent  land  by  the  acre  on  a  cash  basis  or  on  crop-share  basis. 
Preparation  of  the  land  for  crop  production — e.g.,  clearing  cut-over 
pine  land  and  establishing  pasture  may  also  be  included  in  the  cost. 
Further,  the  damage  to  the  environment  and  community  by  use  of  both  the 
land-disposal  process  and  the  alternative  tertiary  treatment  processes 
must  be  borne  by  the  decision  unit.   Thus,  externalities  have  been 
internalized. 

The  objective  of  the  decision  unit  is  to  minimize  total  cost,  TC, 
which  is  defined  as  follows: 

TC   =   C2(E2)  +  d2(E2)  +  c^(E^)  +  d^(E^)  -PyYCE^^,  L,  X)^ 
+  c  (L)L  +  C  X  (6.9) 

subject  to  E  =  E  +  E   or  the  disposal  of  all  waste  generated  by  the 
decision  unit.   Introduce  the  multiplier  and  construct  the  Lagrange 
function 

TC  +  X(E  -  E   -  E^)  (6.10) 


Since  F  is  a  function  of  E^,  equation  (6.2),  therefore  Y(E  ,  L,  X, 
F)  can  be  rewritten  as  Y(E  ,  L,  X) . 
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Taking  the  first  partial  derivatives  and  set  them  equal  to  zero,  one 
obtains 

8TC     3D       9Y 

i  +  — i_p  =   A  (6.11) 

8E     9E^    y  8E^ 

9TC     3D^ 
dE„     3E„ 

8Y  8C 

^aL  =  ^L  +  iT'-  <'•"' 

P^f  =   C,  (6,14) 

These  four  relations,  (6.11)  through  (6.14),  together  x;ith  E  =  E  +  E 
provide  a  solution  for  the  five  unknowns,  E  ,  E  ,  L,  X,  and  X. 

The  terras  on  the  left  of  equation  (6.11)  are  the  total  marginal 
cost  of  employing  a  unit  of  secondary  effluent  in  land  treatment,  which 
include  the  cost  of  transporting  and  applying  waste  E  ,  environmental 
damage,  and  return  from  crops.   The  term,  X,    on  the  right  is  the  inputcd 
marginal  revenue  product  of  the  secondary  effluent.   The  equation  states 
the  familar  proposition  that  the  marginal  cost  of  production  must  equal 
inputed  marginal  revenue  product.   Equation  (6.12)  states  the  same  pro- 
position for  tertiary  treatment  of  E  . 

Equations  (6.11)  and  (6.12)  indicate  that  the  inputed  marginal 
revenue  product  of  secondary  effluent  must  be  the  same  in  each  of  its 
uses,  and  the  total  marginal  cost  for  land  disposal  method  must  be 
equal  to  the  total  marginal  cost  of  disposal  by  tertiary  processes. 

To  show  hov/  the  rate  of  effluent  application  per  acre  fits  into  the 
system,  we  refer  to  Figure  6.1(a).   With  the  price  of  the  crop,  P  , 
fixed,  total  revenue  curve  of  crop,  P  Y,  has  the  same  shape  as  the  crop 

y 

9Y 
yield  response  curve.   The  marginal  revenue,  P  -r——  ,  indicates  the 

y  9E^ 
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(a) 


8TC     3D       8Y 

— L  4.  — L  _  p  — 

9E,     3E,     y  3E 


Figure  6.1    Cost  Curves  of  Disposal  of  Effluent  on  Land 
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addition  to  gross  revenue  if  the  effluent,  E  ,  vjere  costless.   Note 
that  if  we  were  to  viev/  profitable  crop  production  as  the  "sole"  objec- 
tive in  such  a  situation,  ob  is  the  optimal  amount  of  effluent  to  apply. 
If  E   is  not  costless,  how  should  it  be  prices?   To  answer  this,  we 
need  to  go  back  to  equations  (6.11)  and  (6.12). 

Equations  (6.11)  and  (6.12)  are  two  of  the  conditions  for  total 
cost  minimization  for  the  entire  system.   These  conditions  indicate 
that,  as  mentioned  before,  the  total  marginal  cost  for  the  land  disposal 
method  must  be  equal  to  the  total  m.arginal  cost  of  disposal  by  alterna- 
tive methods.   In  Figure  6.1(b),  the  left  side  of  equation  (6.11)  is 
depicted.   For  simplicity,  we  assume  that  the  functions  for  yield  re- 
sponse (6.1),  total  cost,  TC  ,  (6.4),  and  damage  (6.6),  are  quadratic 
in  E  ,  thus  giving  linear  derivatives. 

As  in  Figure  6.1(a),  ob  represents  the  amount  of  effluent,  E  ,  to 
be  applied  for  maximum  crop  yield  per  acre.   Suppose  that  the  total 

marginal  cost  for  the  alternative  method,  ^-r h  ^^   ,  is  oc.   This 

would  require  an  application  of  effluent  in  excess  of  maximum  crop  yield 
per  acre.   If  the  total  marginal  cost  of  the  alternative  method  v;ere  to 
decrease  sufficiently,  the  optimal  application  rate  v;ould  drop  to  less 
than  ob.   Optimal  application  rates  less  than  oa  would  mean  that  dis- 
posal on  land  created,  in  total,  an  income  rather  than  a  cost. 

However,  this  equilibrium  may  not  be  reached  in  some  cases.   For 
example,  in  Figure  6.2,  the  two  vertical  lines  represent  total  marginal 
costs  for  E   and  E  ,  respectively,  and  the  horizontal  line  represents 
the  total  amounts  of  secondary  effluent  disposed  of.   In  panel  (a),  the 
equilibrium  is  reached,  and  an  amount  od  of  effluent  will  be  disposed  of 
by  land  treatment  operation.   The  rest  of  E,  E-od  will  be  disposed  of  by 
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(a)    Example  for  Equilibriuin 


(b)    Example  for  Disequilibrium 


Figure  6.2    Examples  of  Equilibrium  and  Disequilibrium  Conditions 
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the  alternative  method.   In  panel  (b) ,  for  a  given  amount  of  E,  there 
is  no  intersection  for  these  two  marginal  cost  curves.   There  is  no 
equilibrium  involving  both  laethods,  and  land  treatment  provides  the 
least  cost  method  of  disposing  of  all  the  effluent.   Of  course,  the  re- 
verse situation  may  also  occur  if  the  marginal  cost  of  the  land  treat- 
ment operation  is  dominated  by  the  marginal  cost  of  the  alternative 
method. 

Equations  (6.13)  and  (6.14)  state  that,  under  equilibrium  condi- 
tions, the  value  of  marginal  product  of  a  factor  equals  the  marginal 
cost  of  this  factor.   As  indicated  by  equation  (6.8),  the  price  of  land 
m.ay  be  a  function  of  the  amount  of  land  required  for  land  treatment 
operations.   The  right  side  terms  of  equation  (6.13)  shov;  that  the 
marginal  cost  of  land  acquisition  is  a  function  of  the  amount  of  land 
acquired,  and  equation  (6.14)  shows  the  marginal  cost  for  other  inputs 
to  be  a  constant,  C  . 

X 

Of  course,  other  variables  besides  the  marginal  cost  of  alterna- 
tive disposal  methods  could  affect  the  optimal  application  rate;  for 

example,  the  other  two  conditions  sho\-7n  by  equations  (5.13)  and  (6.14). 

3P  Y 

The  marginal  value  product  of  effluent,  — ;r —  ,  also  has  important  im- 

8E^ 

plication  for  the  price  a  farmer  can  afford  to  pay  for  varying  amounts 
of  effluent,  E  .   If  crop  production  is  viewed  as  a  separate  decision- 
making unit,  the  am.ounts  of  effluent  less  than  ob  have  a  positive  value 
and  could  command  a  price;  application  in  excess  of  ob  would  cause  fi- 
nancial loss  to  the  crop  enterprise  and  v.'ould  require  compensation. 

This  model  can,  of  course,  be  expanded  to  include  more  alternative 
treatments,  crop  systems,  and  constraints.   Hovrever ,  in  any  cost  mini- 
mization model  the  optimal  rate  of  application  of  effluent  to  cropland. 


102 


E  /L,  depends  on  knowledge  of  the  functional  relationships  (6.1)  through 
(6.8)  and  prices  and  costs. 

Although  in  Florida  a  number  of  experiments  have  been  conducted  to 
show  the  effect  of  irrigating  wastev;ater  on  pasture  and  selected  crops, 
very  few  have  made  estijiates  of  the  yield  response  to  v;astewater  and 
fewer  experiments  have  been  conducted  to  develop  new  management  tech- 
niques for  a  land  treatment  system  [Alexander,  1972;  Hor tenstine,  1973; 
Institute  of  Food  and  Agricultural  Sciences,  1973;  Overman,  1971; 
Overman  and  Smith,  1973].   Such  management  techniques  may  include  the 
following:   what  crop  should  be  grown,  what  type  of  irrigation  system 
is  best  for  the  desired  agricultural  enterprise,  what  is  the  best  ir- 
rigation frequency  for  disposing  of  v.-astewater,  what  is  the  harvesting 
requirement  for  forage  crops  (grasses)  and  how  often  is  removal  of  the 
forage  required.   The  use  of  crops  and  land  to  renovate  v;astev;ater  and 
thus  eliminate  unwanted  materials  through  so  called  "living  filters" 
is  different  from  the  conventional  problems  in  agricultural  production. 
In  a  conventional  agricultural  production  problem  one  tries  to  maximize 
net  returns  subject  to  limited  resources  or,  on  the  other  hand,  to  use 
minimum,  amounts  of  resources  to  produce  a  certain  level  of  returns.   The 
emphasis  has  been  placed  on  agricultural  production  research.   While  in 
land  treatment  system.s,  the  emphasis  is  placed  on  the  disposal  of  waste- 
water and  agricultural  production  plays  a  secondary  role.   As  a  result 
there  is  a  serious  lack  of  knowledge  of  management  techniques  in  the 
areas  of  negative  marginal  productivities  of  resources. 

Based  upon  the  input-output  relationships  developed  by  the  IFAS 
Irrigation  Task  Force,  Holt  [1972]  has  provided  estimates  of  the  produc- 
tion practice  necessary  to  produce  crops  under  irrigation  on  formerly 


103 


unf armed  deep  sand  soils  in  North  and  West  Florida.   A  center  pivot 
irrigation  system  is  used  to  cover  a  138~acre  area  in  his  study,  the 

initial  investment  estimate  is  $200  per  acre,  and  the  variable  cost 

2 
estimate  is  $1.25  per  acre  per  inch  of  water  applied.    Table  6.1  shov7S 

the  estimates  for  six  selected  irrigated  crops.   Column  two  shovjs  that 

for  six  crops,  the  returns  over  specified  expenses  are  all  negative 

except  sorghum  for  silage  using  the  1971  unit  prices  given  in  column 

two.   With  later  prices  the  crops  would  be  more  profitable.   Column 

four  shov/s  the  cash  expenses  on  fertilizer  for  the  crops. 

Overman  and  Smith  [1973]  have  studied  the  nutrient  composition  of 
secondary  effluent  at  the  southwest  treatment  plant  in  Tallahassee, 
Florida.   Table  6.2  shows  the  nutrient  composition  of  secondary  effluent 
in  their  study  and  the  amounts  of  nutrient  content  in  the  effluent  if 
the  irrigation  rate  is  tv7o  inches  per  week  for  twenty-five  v;eeks.   If 
these  nutrients  can  be  used  to  eliminate  the  costs  of  fertilizer,  and  if 
there  are  no  adverse  effects,  then  it  would  be  profitable  for  all  the 
six  crops  mentioned  in  Table  6.1  except  soybeans  for  beans  (see  Column 
5).   Further,  if  calcium  is  obtained  from  secondary  effluent,  all  the 
six  crops  would  be  profitable.   However,  the  possibility  of  having  ad- 
verse effects  by  the  application  of  excess  v/ater  and  nutrients  were  not 
considered  by  Overman. 

Table  6.3  shows  the  estimated  total  expenses  per  acre  of  the  irri- 
gated crops  mentioned  above  by  Holt  [1972].   Table  6.4  shows  the  land 
preparation  and  irrigation  system  costs  from  three  studies.   Based  on  a 


2 
Based  on  an  irrigation  frequency  of  5  days,  6  inches  annually,  and 

total  initial  investment  includes  deep  v/ell,  pump,  power  unit  and  system, 

but  no  depreciation. 
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Table  6.2   Nutrient  Composition  of  Secondary  Effluent  at  Southwest 
Treatment  Plant  in  Tallahassee 


Element  Ib./acre-in 


.  a  Irrigation  Rate  (2  in/wk.) 


(lb . /acre) 


N  8.2  410 

P  2.7  135 

K  1.4  70 

Ca  7.1  355 

Mg  2.2  110 

Na  8.9  445 
Cu              less  than  0.02 

Fe  0.11  5.5 

Zn  0.033  1.6 


From  composite  samples  for  period  April  5,  1972  to  September  27, 
1972. 

Based  on  twenty-five  weeks  of  irrigation. 

Source:   Overman  and  Smith  [1973]. 
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640-acre  irrigation  site,  Muraro  [1972]  has  estimated  the  costs  for  a 
land-spreading  systein.   A.nderson  and  Hipp  [1973]  have  estimated  the 
costs  of  land  preparation  for  selected  pastures.   And  Holt  [1972],  as 
mentioned  before  has  estiraates  for  central  pivot  irrigation  system. 

For  a  140-acre  field,  Westberry  [1974]  has  estimated  the  costs  of 
ovming  and  operating  a  central  pivot  irrigation  system,  and  the  cost  of 
corn  production  v/ith  and  without  irrigation.   According  to  his  estima- 
tion, the  cost  per  bushel  of  corn  is  $1.70  v^ith  irrigation,  and  $1.75 
without  irrigation  with  estimated  productions  of  155  bushels  par  acre 
and  60  bushels  per  acre,  respectively.   His  results  show  that  a  corn 
grov7er  receives  a  higher  expected  revenue  with  an  irrigation  system. 

The  results  from  Holt  [1972]  and  Westberry  [1974]  show  that  irri- 
gation has  the  potential  not  only  to  increase  yield  but  to  provide 
insurance  against  the  possibility  of  serious  yield  reduction  due  to 
prolonged  dry  periods. 

The  most  important  unknov/n  is  whether  an  agricultural  land  treat- 
ment system  will  adequately  purify  water  to  meet  acceptable  standards. 
If  so,  the  cost  of  such  treatment  is  considerably  belov?  that  for  terti- 
ary treatment.   If  one  assumes  that  on  site  a  land  treatment  operation 
can  break-even,  one  can  estimate  how  far  secondary  effluent  can  be 
pumped  to  a  land  treatment  site  and  still  compete  with  tertiary 
treatment. 

Westberry  [1974]  has  estimated  the  fixed  cost  of  a  center  pivot 
irrigation  system  as  $51.07  per  acre  annually,  and  a  variable  cost  of 
$2.06  per  acre  inch  of  water  applied.   A  center  pivot  system  operates 
in  a  160-acre  block  and  irrigates  a  circular  area  of  138  acres. 
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Baldwin  [1975]  has  estimated  a  fixed  cost  of  $96  per  acre  annually  and 

a  variable  cost  of  $1.65  per  acre  inch  of  water  applied  for  such  a 

3 

system. 

Table  6.5  shows  pipeline  and  pum.ping  station  costs  for  time  options 
2  and  5  (Table  5.1)  and  the  estimated  irrigation  system  cost  from 
Westberry  [1974]  and  Baldwin  [1975]. 

From  these  cost  figures,  one  finds  that  time  option  2  provides  a 
minimum  cost  of  pipeline  per  mile  for  the  cities  of  High  Springs, 
Newberry,  Archer,  Alachua,  and  the  Gainesville  Urban  Area,  and  time  op- 
tion 5  provides  a  minimum  cost  for  LaCrosse,  Micanopy,  Waldo,  and 
Hawthorne.   The  different  time  option  for  the  latter  four  cities  is 
probably  due  to  their  small  population  changes  betvjeen  1975  and  1990. 

Wastev/ater  estimates  for  all  cities  in  Alachua  County  except  those 
for  the  Gainesville  Urban  Area  were  small.   Therefore,  time  option  5  vjas 
used  to  calculate  Irrigation  system  costs  for  all  communities  except  for 
the  Gainesville  Urban  Area.   If  other  time  options  had  been  used,  some 
flows  V70uld  have  been  too  small  to  be  realistic  for  the  construction  of 
an  irrigation  system.   Time  option  2  was  used  for  the  Gainesville  Urban 
Area.   All  systems  costs  include  construction  and  operation  and  mainte- 
nance costs. 

The  difference  between  the  cost  of  tertiary  treatment  and  site  costs 
for  land  spreading  were  used  to  estimate  break-even  distances  over  which 
effluent  could  be  moved  for  land  spreading  purposes  (Table  6.6). 
Baldv.'iu's  [1975]  cost  estimates  (column  6  in  Table  6.5)  and  a  high-rate 
trickling  filter  system  were  used  for  the  calculation.   Break-even 


3 

These  cost  figures  were  derived  from  Harrison's  [1975]  work. 
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distances  show  hov/  far  one  could  transport  the  secondary  effluent  from 
each  source  for  land  spreading  and  be  indifferent  between  land  spreading 
and  tertiary  treatment. 

Table  6.6  shows  the  break-even  distances  for  tir.ie  options  2  and  5. 
These  results  show  that  secondary  effluent  can  be  transported  considera- 
ble distances  for  land  treatment  and  still  compete  v/ith  tertiary 
processes. 

In  estimating  break-even  distances,  it  "was  assumed  that  Newberry, 
Archer,  and  Alachua  shared  regional  tertiary  treatment  costs  in  propor- 
tion to  their  v/astewater  flov7S.   The  shares  of  their  wastev7ater  flow  are 
0.2423,  0.1805,  and  0.5773,  and  their  corresponding  assessment  in  1975 
was  assum.ed  to  be  $579,569,  $432,747,  and  $1,380,872  respectively. 

If  these  estimates  are  cox'rect,  and  if  soil  and  crops  could  elim.i- 
nate  nutrients  and  unwanted  materials  from  secondary  effluent  effec- 
tively, then  the  cost  for  land  treatment  V70uld  be  the  difference  between 
the  site  preparation,  transmission  and  pumping  costs,  and  the  net  re- 
turn from  crop  productions,  which  is  considerably  less  expensive  than 
tertiary  treatment. 

Relationships  (6.2),  (6.4),  and  (6.5)  involve  basically  engineering 
estimates  and  it  would  be  expected  that  the  parameters  in  these  equa- 
tions could  be  estim.ated  with  relatively  greater  precision  than  (6.1) 
which  is  an  essentially  biological  relationship.   If  (6.5)  refers  to  a 
tertiary  system  of  treatment  and  disposal,  one  would  expect  that  the 
estimates  in  previous  chapters,  together  with  the  updating  of  costs, 
would  give  rather  precise  knowledge  of  this  relationship. 

Knowledge  is  most  meager  in  the  case  of  damage  functions,  (6.6)  and 
(6.7).   Although  these  damages  are  expressed  in  dollars,  in  the  case  of 
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Table  6.6   Break-even  Length  of  Transmission  Pipeline  for  Land  Treat- 
ment in  Alachua  County,  Florida,  1975  to  1990 


City 

Tertiary 

Process  Cost 

Breal. 
Distance 

c-even 
(miles) 

Option  2 

Option  5 

Option  2 

Option  5 

High  Springs 

1,313,235 

1,467,250 

32.68 

36.94 

Nev/berry 

598,880 

658,683 

20.40 

22.55 

Archer 

448,878 

500,970 

16.82 

18.91 

Alachua 

1,425,237 

1,651,628 

33.69 

39.78 

Gainesville 

Urban  Area 

25,723,425 

31,205,441 

100.39 

117.04 

LaCrosse 

317,198 

358,025 

19.45 

22.04 

Micanopy 

470,352 

479,559 

24.23 

24.72 

Waldo 

488,925 

518,472 

23.94 

25.46 

Hawthorne 

650,697 

696,263 

26.35 

28.31 

Calculated  from  Tables  5.5,  5.6,  and  5.7. 


Calculated  from  Table  6.5.   These  figures  vzere  obtained  by  sub- 
tracting costs  of  pumping  station  and  Baldwin's  irrigation  system  cost 
from  tertiary  process  cost  then  divided  by  the  cost  of  pipeline  per 
mile. 
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many  of  the  potenLially  harmful  aspects  of  applications  of  wastewater 
the  fundamental  physical  relationships  are  not  knovra.   Damages  may  also 
be  difficult  to  measure  in  dollars.   For  example,  the  life  of  a  person 
is  difficult  to  express  in  dollars,  and  the  transformation  in  space  and 
time  of  heavy  metals  represents  processes  that  are  not  well  understood. 
It  is  for  such  reasons  that  policy  normally  attempts  to  reduce  potential 
adverse  environmental  impacts  by  setting  rather  subjective  standards. 

Functions  (6.6)  and  (6.7)  may  not  be  linear,  since  these  damages 
depend  not  only  on  the  amount  of  effluent  disposed  but  also  on  the  rate 
of  application.   For  example,  if  tertiary  treatment  is  used,  the  damages 
shown  by  (6.6)  are  functions  of  the  amount  of  effluent  disposed  as  well 
as  the  effluent  disposal  rate.   The  ability  of  the  environment  to  re- 
store itself  is  limited  and  more  or  less  constant  over  a  given  period  of 
time.   If  a  regional  treatment  plant  treats  wastewater  for  the  whole 
county  and  releases  its  tertiary  effluent  at  one  location,  the  rate  of 
pollutants  released  at  this  location  is  still  tremendous  even  though 
the  pollutant  contents  of  the  effluent  are  low,  and  hence  create  an 
environmental  hot  spot.   The  same  result  may  happen  if  land  treatment 
is  used.   This  damage  to  the  environment  may  discourage  the  region- 
alization  of  wastewater  treatment  facilities. 

The  management  and  operation  of  wastewater  treatment  facilities 
add   to  county  incomes  both  directly  and  indirectly.   The  proportions 
of  these  expenses  that  go  into  the  county's  business  and  household 
sector,  respectively,  determines  the  final  effects  of  treatment  on 
county  income.   The  indirect  effects  of  treatment  on  county  income  are 
important  and  should  not  be  ignored.   Since  agricultural  production 
normally  creates  considerable  indirect  activity  one  V70uld  expect  that 
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tlie  land  treatment  alternative  would  have  greater  beneficial  impact  on 
the  community  than  v.'ould  tertiary  treatment. 

The  state  of  the  arts  in  advanced  wastewater  treatmeiit  is  stil]  in 
its  developmental  stage.   There  remains  considerable  uncertainty  as  to 
the  effectiveness  of  both  chemical  and  land  spreading  m.ethods  for  ad- 
vanced treatment.   The  effectiveness  of  land  spreading  see'ms  to  be  site 
specific  and  a  function  of  such  factors  as  type  of  soil  and  slope  of 
land.   The  effectiveness  of  tertiary  processes  in  a  real  vjcrld  treatment 
situation  also  remains  in  question. 


CHAPTER  VII 
SUMMARY  AND  DISCUSSION 


Summary 


In  the  planning  and  design  of  a  wastewater  treatment  facility  for 
a  city,  town,  or  an  area  v/here  the  requirements  for  treatment  are  ex- 
pected to  increase  V7ith  time,  the  initial  size  of  the  treatment  plant 
and  the  timing  of  capacity  additions  and/or  replacements  over  some  time 
horizon  need  to  be  answered  in  the  context  of  an  optimal  staging  policy. 

The  economies  of  scale  provide  the  main  incentive  for  regionaliza- 
tion  for  wastevjater  treatment  facilities  as  well  as  for  constructing 
larger  initial  treatment  plant  sizes  and  the  timing  of  capacity  addi- 
tions and/or  replacements  over  some  time  horizon.   Several  studies  have 
been  made  that  have  discovered  the  existence  of  economies  of  scale  in 
both  the  construction  and  operation  and  maintenance  costs  of  w^astev/ater 
treatment . 

These  econom.ies  of  scale  come  mainly  from  the  following  sources. 
In  the  process  industries,  such  as  petroleum  refining  and  chemical  pro- 
duction, the  output  of  a  processing  unit  tends  v^^ithin  certain  physical 
limits  to  be  roughly  proportional  to  the  volume  of  the  unit,  other 
things  being  equal.   The  amount  of  materials  and  fabrication  effort  re- 
quired to  construct  the  unit  is  more  apt  to  be  proportional  to  the  sur- 
face area  of  the  unit's  reaction  chambers,  storage  tanks,  connecting 
pipes,  etc.   Since  the  area  of  a  sphere  or  cylinder  varies  as  the 
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two-thirds  power  (or  six-tenths)  of  volume,  the  cost  of  constructing 
process  industry  plants  can  be  expected  to  rise  as  the  tvjo-thirds  power 
of  their  output  capacity.   By  the  same  token,  there  are  economies  of 
scale,  in  operation  and  maintenance  costs.   In  addition,  a  host  of  other 
advantages  such  as  more  qualified  operating  personnel,  and  higher 
degrees  of  automation  nay  be  gained  through  the  utilization  of  this 
concept . 

It  is  quite  clear  that  economies  of  scale  do  exist  in  vjastewater 
treatment.   However,  there  are  several  reasons  why  unit  costs  do  not 
decline  indef initel^'.   These  reasons  include  the  difficulty  of  managing 
and  operating  an  enterprise  of  an  increasingly  large  scale  volume  of 
production,  and  the  increases  in  the  transportation  costs  of  raw  sewage 
for  large  scale  treatment  plants.   In  this  study,  the  emphasis  has  been 
placed  on  the  trade-offs  between  the  gains  in  economies  of  scale  in 
treatment  and  the  increase  in  transportation  costs  of  raw  sewage  over 
space,  and  the  trade-offs  between  the  gains  of  economies  of  scale  in 
treatment  aiid  the  increase  in  construction  costs  over  time. 

Alachua  County,  Florida  vras  chosen  as  the  study  area.   Population 
estimates  of  years  1975,  1980,  and  1990  for  each  of  the  nine  incorpo- 
rated cities  in  the  county  were  obtained  from  the  estimates  provided  by 
the  North  Central  Florida  Regional  Planning  Council  [1973].   The  re- 
quired treatment  capacities  for  1975,  1980,  and  1990  for  the  nine 
cities  were  estimated  by  multiplying  population  projections  by  135 
gallons  per  capita  per  day — an  average  sewage  flow  provided  by  Lochr 
[1968]. 

Two  types  of  secondary  treatment  and  three  tertiary  treatment  pro- 
cesses were  considered  throughout  the  study:   high-rate  trickling  filter 
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system,  activated  sludge  system,  and  tv/o  clarifier  lirae  clarification 
processes,  lime  recalcination,  and  ammonia  stripping  process, 
respectively.   The  relationsliips  between  capacity  and  construction  and 
operation  and  maintenance  costs  were  obtained  from  EPA  studies  [EPA, 
1973a;  Michel  and  Johnson,  1970;  Smith  and  McMichael,  1969;  Smith  and 
Eilers,  1971]. 

The  concept  of  regionalization  has  two  aspects:   one  is  the  number 
of  separate  sewerage  s^'stems  operated  by  a  single  authority,  i.e.,  the 
centrali/^ation  of  administration;  and  the  other  is  a  single  connected 
sewer  system  serving  several  municipalities  with  one  central  plant. 

An  optimization  model  for  the  search  of  minimum  cost  combination 
of  regional  wastewater  treatment  plan  over  space  and  time  was  developed. 
This  model  provided  a  guideline  for  the  development  of  a  heuristic 
method  of  developing  an  optimal  staging  process  for  wastev;ater  treat- 
ment in  Alachua  County. 

The  space  and  tim.e  dimensions  of  regionalization  were  considered. 
In  order  to  find  the  potential  sites  of  regional  treatment  plants,  a 
cost  minimization  model  was  used.   For  the  consolidation  over  time  an 
exponential  objective  function  was  minimized  subject  to  some  discrete 
capacity  requirements. 

The  rates  of  increase  in  cost  over  time  were  estimated  from  EPA 
cost  indexes  as  .12,  .10,  and  .07  per  year  for  force  main  construction 
cost,  pumping  stations  and  sewage  treatment  plant  construction  costs, 
and  O&M  costs  for  both  force  main  and  sewage  treatment  plant, 
respectively.   Service  lives  for  pumping  station  and  wastewater  treat- 
ment plants  were  taken  as  25  years,  and  for  pipelines  as  50  years.   An 
interest  rate  of  5%  per  year  v/as  assumed. 
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A  least  cost  staging  policy  to  satisfy  the  wastewater  treatment 
requirements  of  the  nine  cities  in  Alachua  County,  Florida  for  the 
years  1975,  1980,  and  1990  Vi/as  determined  with  the  model.   The  results 
were  sensitive  to  the  required  quality  of  secondary  effluents. 

If  the  required  quality  of  secondary  effluent  is  set  at  above  80% 
removal,  activated  sludge  plants  were  required  to  provide  the  least 
cost  method  of  treating  wastexijater.   The  results  show  that  secondary 
treatraent  cost  w^ill  be  minimized  if  High  Springs,  Alachua,  Archer,  and 
Newberry  cooperate  to  build  a  regional  treatment  plant  in  1975  to  sat- 
isfy their  1980  demands,  and  then  construct  individual  plants  in  1990 
to  satisfy  the  additional  demands  at  that  time.   The  other  cities  could 
satisfy  their  demands  by  building  individual  activated  sludge  plants  in 
1975  to  satisfy  their  1980  demands,  v/ith  additions  being  made  in  1990. 
The  minimum  tertiary  treatment  cost  combination  v/ould  require  the  con- 
struction of  two  regional  treatment  plants,  one  for  High  Springs, 
Alachua,  Archer,  and  Newberry,  and  one  for  Waldo  and  Hawthorne  in  1975 
to  satisfy  the  1980  demands  of  these  six  cities.   Their  additional  de- 
mands in  1990  would  be  met  with  plants  in  each  city.   Gainesville, 
LaCrosse,  and  Micanopy  v/ould  satisfy  their  tertiary  treatraent  demands 
vjith  individual  plants  constructed  in  1975,  with  additions  in  1990. 

If  the  required  quality  of  secondary  effluent  is  set  at  80%  re- 
moval, the  results  shov;  that  secondary  treatment  cost  for  the  county 
will  be  minimized  if  each  city  builds  a  trickling  filter  plant  which 
satisfy  1980' s  wastewater  treatment  demands  in  1975,  with  additions  to 
capacity  being  made  in  1990.  The  minimum  tertiary  treatment  cost  com- 
bination suggests  that  there  v/ill  be  regionalization  of  treatment  in 
1975  to  satisfy  the  1980  demands  of  three  cities — Newberry,  Archer, 
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and  Alachua.   Additional  demands  in  1990  will  be  satisfied  by  building 
individual  treatment  plants  in  these  three  cities.   The  other  six 
cities  would  satisfy  their  tertiary  treatment  with  individual  plants 
constructed  on  the  same  schedule  as  the  trickling  filter  plants. 

One  of  the  objectives  of  this  study  has  been  the  study  of  the 
trade-offs  between  the  construction  of  tertiary  treatment  processes  and 
the  usage  of  secondary  effluent  as  irrigation  water.   Under  conventional 
secondary  treatment  about  90%  of  BOD  is  removed.   Almost  no  nitrogen  and 
phosphorous  are  removed.   Thus,  secondary  effluent  may  cause  problems  of 
overgrov;th  of  vegetation  in  receiving  waters.   Cost  of  removing  wastes 
above  95%  increase  rapidly.   Land  treatment  by  irrigation,  through  the 
uptake  of  nutrients  by  crops  or  grass  may  attain  a  higher  removal  of 
BOD  and  COD  and  remove  up  to  98%  of  the  nutrients  in  secondary  effluent. 
However,  precise  knowledge  of  the  land  treatment  process  is  not 
available. 

The  possibilities  of  substituting  land  for  tertiary  treatment  vjas 
discussed,  and  the  factors  which  should  be  considered  in  land  treatment 
were  discussed  with  the  aid  of  an  optimization  model.   The  break-even 
distances  of  transmission  pipelines  for  land  treatment  as  a  substitute 
for  tertiary  treatment  were  calculated  and  presented  for  the  nine 
cities  in  Alachua  County,  Florida. 

Discussion 

Factors  such  as  arrangement  of  outfall  structure,  cost  allocation 
among  participants,  methods  of  financing,  and  methods  of  implementation 
have  not  been  considered  in  this  study. 
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The  main  purpose  of  this  study  was  to  consider  the  time  and  space 
dimension  of  regional  wastev/ater  treatment  planning.   Because  of  dis- 
tances between  cities,  the  gains  in  economies  of  scale  in  the  construc- 
tion of  regional  treatment  plants  to  serve  more  than  one  city  are  more 
than  offset  by  transmission  costs  of  rav.'  sev/age  from  each  source  to  po- 
tential treatment  sites.   Therefore,  the  minimum  cost  solution  of  this 
study  includes  tv/o  regional  treatment  plants  for  activated  sludge  sys- 
tem, and  one  regional  treatment  plant  for  trickling  filter  system. 

The  distance  between  each  tv/o  cities  v;ere  such  that  only  ten  com- 
binations of  cities  v;ere  selected  by  the  cost  minimization  model.   In 
studies  vTliere  the  cities  are  closer  together,  the  scanning  process  used 
by  this  study  may  create  a  larger  number  of  potential  sites.   If  so,  the 
search  to  find  the  optimum  sites  for  treatment  plants  vrill  be  more  in- 
volved and  may  require  the  imposition  of  other  constraints. 

Requirements  for  the  transmission  of  vjastewater  depend  on  the  topo- 
graphic condition  of  the  area  studied.   In  the  design  of  pipelines  and 
pumping  stations,  topographic  conditions  should  be  taken  into  considera- 
tion.  In  this  study,  it  was  assumed  that  a  pumping  station  was  required 
for  each  regional  plant;  of  course,  more  than  one  pumping  station  may  be 
needed  in  some  cases.   The  process  of  scanning  all  the  possible  sites 
makes  it  difficult  to  take  the  topographic  factor  into  consideration  for 
each  site.   However,  if  a  site  is  excluded  from  the  solution  by  the 
assumption  of  only  one  pumping  station,  it  v/ould  not  appear  in  the  solu- 
tion if  more  than  one  pumping  station  were  required.   Therefore,  the 
assumption  is  not  critical  to  the  results. 

A  cost  minimization  model  for  a  land  treatment  system  was  presented. 
And  the  break-even  distances  over  which  one  can  transport  secondary 
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effluent  from  sources  to  land  treatment  sites  as  an  alternative  to 
tertiary  treatment  v/ere  estimated.   An  important  study  topic  is  the 
income  generated  by  crop  production  through  the  disposal  of  secondary 
effluent,  and  how  this  income  is  distributed.   This  topic  has  important 
implications  for  rural  and  urban  developm.ent. 

T.and,  including  easement  costs,  was  not  estimated  in  this  study. 
This  is  £1  very  important  variable  in  the  actual  implementation  of  a 
treatment  program--particularly  for  land  treatment. 

The  impacts  of  energy  crisis  on  wastewater  treatment  as  well  as  on 
land  treatment  were  not  discussed.   However,  one  may  suspect  that  the 
effects  of  the  energy  crisis  v;ould  increase  all  costs  proportionally, 
and  hence,  the  solution  of  this  study  may  not  be  affected. 

Another  aspect  of  regionalization,  centralization  of  administra- 
tion, was  not  considered  for  the  situation  in  Alachua  County.   Since 
all  the  cities  are  not  far  from  each  other  and  their  wastewater  treat- 
ment demands  are  low,  the  centralization  of  management  and  administra- 
tion may  be  feasible  to  utilize  personnel  m.ore  thoroughly. 
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COMPUTER  PROGRAM  FOR  THIS  STUDY 


123 


13 


10 


72 


73 


75 


74 


75 


77 


COf'^MCN    r'M 

COf-'NGK    CG«(2,9),WT{3,9),CAH(3,9),CI{?,     3,9)»P0P 
♦TRC(3,9),TC(3,9),RTRT{3).RCP(3),RTC(3),RPP(3),R 

*  TC { 3 , 9 ) , RT RT ( 3 ) t  RCP ( 3 ) , R T C ( 3  ) , RPP { 3  ) , R  1  R ( 3 )  , 

*  P  .V  P  (  3  ,  9  )  ,  C  X  (  9  )  ,  C  Y  (  9  )  ,  P  P  {  3  , 9  )  ,  W  T  T  (  3  ,  9  )  ,  D  I  S  (  9  )  ,  C 

*  CCl{2,3,10),COMl(2,3,10),TCCI{2,3,9),TCnfn(2,3 

*  PC0fn(2,3,9)  ,ACaP(3,9),ACP{3,9),RACP(3)        ,A;.T( 
:^^    TAl^T{3,9),TRT2(2,3,lC),TRT3(2,3,10),FX2(2,3),F 

*  TRC2{3t9),PMP2{3,9)  , RPP2 ( 3 )  , P?2 ( 3 , 9  )  , FX3{ 2 , 3 ) , 

^    2,3,9) ,A^TCM(3,10) 
COtVMON    WT2  I  3  ,  9  )  ,  ViTT2  (  3  ,  9  )  ,  RT  C  2  (  3  )  ,  TC  2  [  3  ,  9  } 

DIMENSIC\  IGU(9) 

READ (5,1)  < (CCR( I vJ) t 1=1«2), J=lr9) 

READ  (5,2)  (  (PCP{  I  ,J)  ^1  =  1.3),  J  =  l-,  9) 

FGRf'5AT(2F5.3) 

FORMAT(3F6oO) 

DO  13  1=1,2 

DO  13  J-1,9 

CCR(l,J)=CCR( I,J)/2.625 

00  39  r-:M  =  l,4 

DO  10  1=1,3 

DO  10  J=l,9 

ACAP( I, J)=PCP( I, J )*13  5.0/1COCCCO. 

CAP{ 1, J)=PCP( I,J)*135.0/1CG0CC0. 

GO  TC  (72, 73, 74, 74), PM 

CAP(3,5)-CAP(3,5)-5.C 

CAP(2,5)=CAP(2,5)-9.5 

CAP(  1,5)--^CAP(  l,5)-9.5 

GC  TC  77 

DC  75  J=l,9 

ACAP{3,J)=ACAP(3,J )-ACAP{ 1, J ) 

ACAP(2,J )  =  ACAP(2, J)-ACAP(  1,J) 

CAP(3,J)=CAP(3,J)-CAP{ 1, J) 

CAP(2,J)=CAP( 2,J)-CAP(1,J) 

CAP( 1,5)=CAP( l,5)-5, :-4.5 

CAP(3,5)=CAP{ 3,5)+4.5 

GO  TC  7  7 

DC  76  J=l,9 

ACAP(3,J)=ACAP(3, J)-ACAP{2,J ) 

ACAP(2,J)=ACAP{2, J)-ACAP{ I, J) 

CAP(3,J)=CAP( 3,J)-CAP(2,  J) 

CAP(2,J)=CAP{2,J)-CAP(1,  J) 

CAP(1,5)=CAP( 1,5) -5. ^-4, 5 

CAP(3,5)  =  CAP( 3,5)+4=  5 

GC  TC  7  7 

DC  11  1=1,3 

DC  11  J=ly9 

TRC2( Ir J )=CAP ( I, J )-3.5 

PKP2{  I,  J)-Pr:G(TRC2(  I  ,J)  ,  I,J,  I) 

WT2(  i,J)  =  TRf-UTRC2(l  ,  J)  ,  I  7  Jt  n 

TRC( I , J)  =  ACAP (  I, J  )-3.5 


(379)  , 
TR(3) , 

P  (  3  ,  9  )  , 
,9)  , 
3,10), 
X(2,3), 
PCCK 


124 


11 

12 
3 

41 
42 


81 

142 

46 


48 
82 
85 
49 
50 
44 
22 
23 
26 
51 
52 
45 
24 
25 


PMPl I , J) 
WT( I , J)^ 
DO  12  K= 
DO  12  1= 
DO  12  J= 
CALL  TRT 
WRITE (6, 
FORMAT { • 
WRITE (6, 
FORMAT ( ' 
'^^■RITE(6, 

for:-;  AT  { • 
■-    ) 

WRITE (6, 

FORKAT( • 

=  ) 

W  K IT  e ( 6  » 
FORMAT! • 
WRITE (6 
FORMAT ( 
WRITE (6 
FORMAT ( 

=  2)  ) 
WRITE (6 
FORMAT! 
WRITE (6 
FORKATC 

WRI rE(6 
FORMAT ( 
WRITE (6 
FORKAT( 
'WRITE  (6 
FORMAT ( 
WRITE  (6 
FORN'AT! 
WRITE (6 
FORKAT( 
Ij  R  I  T  E  I  6 
FORMATt 
WRITE (6 
FORKAT( 
'WRITE  (6 
FOR['AT( 
WRITE (& 
FORN.ATC 
V'JRITE{6 
FORN'ATl 
'rtRIT£{6 
FORMAT ( 
w  R I  T  E  {  6 
FURf''AT( 


=  rr'G(TRC(  I,J)  ,  I,  J,2) 

TUi';  (TRC{  I  J  J)  ,  I,  J,  2) 

1,2 

It3 

1»9 

(K,  I, CAP (  I,  J  ),ACAP( I  5  J) J  J) 

3) 

1'  ) 

4) 


(  (  POP(  I,J)  ,J-U9)  ,  1-1,3) 
t///'     PCPULATIC;,^'/'     S3(/- 
)     (  (COiU  It  J)  »J-lt9)  ,  1-1,2) 
///•    CORDliVATG     IN    FILES',/ 


42)  {  {CAP (  I, J), J  =  1,9)  , 1=1,3) 
'///'  TRT  CAPACITY  IM  MGC'/ 


'  ,  9  F  1 4  ,  0  )  ) 
',2(/'  ',9F14.4) 

•  ,3{9F14.4,/'  '} 


8  3  )  (  (  ACAP  (I  ,  J  )  ,  J  -^  1 ,  9  )  ,  I  =  1 ,  3  ) 
'///•  AWT   CAPACITY  S3(/'  «,9F]4.A)) 

Bl)  ( ( (TRT2(K, I,J),J=l,9),I=lt3) ,K=1,2) 
'///•  INDIVIDUAL  TRT2  CGST»t6(/'  «»9rl4,2)) 

A6)((CC1(K,I,J),J=1,9),I=1,3) 
'///•  PRESE;a  TRT  CCNST  COST'/'  ',3{/«  'v9F14 


a)  (  ( 
'/// 
Z)     ( 

v// 

B5)  I 

-:;)(( 

0)  (( 

•/// 

4)  ( 
'/// 

2)  (  ( 
'/// 

3)  {  ( 
•/// 

b)  (  ( 
'/// 

1)  (  ( 
•/// 

2)  (  ( 
'/// 

5)  ( 
'/// 
't)  {  ( 
•/// 
*;>)(( 


CGMl 
'  PR 
{  (  T  R 

•  IN 
IPMP 
'  CO 
TCCl 
»  PR 
TCOM 

•  P  R 
{  WT  ( 

•  TR 
TCCl 
'  TR 
TCOM 
'  TR 

pyp2 

•  Tk 
PCCl 
'  PR 
PCOM 

•  PR 
(PHP 
«  PU 
PCCl 
'  TR 
PCOM 


(K,I 
ESEN 
T  3  (  K 
CIVl 
2(1, 
NVEK 
(1,1 
ESEN 
1(1, 

ESE^ 
I, J) 

A\'sr- 

'«  2  y  1 

T3TR 
1(2, 
T3TR 
(It  J 
T?  P 
(1,  I 
£SE,\ 
Kit 
ESEK 
(  I,  J 
KPIN 
(2,1 
T3  P 
1(2, 
T3  P 


,J)  , 
T  TR 
,  I,  J 
DUAL 
J), J 
TION 
,J), 
T  TR 
I, J) 
T  TR 
,J=1 
ISS  I 
,J)  , 
M  CO 
I, J) 
M  0  + 
),J  = 
KG  C 
,J), 

r  Piv 
I,  J) 

T  P?-' 
),J- 
G  CO 
tJ), 
MP  C 
I, J) 


J  =  l, 
T  0  + 
),J  = 
TRT 
=  1,9 
AL  F 
J  =  l, 

^^  CO 
,J  =  i 

M  0  + 
,9), 
ON  C 
J  =  1t 
\ST 
,J=1 
K  CG 
1,9) 
CST  • 
J-I, 
P    CO 

,J=1 

P  0  + 
1,9) 
ST  ' 
J---1, 
C^;ST 
,J=1 


9),  I  =  lv 
N  CCST' 
1,9), 1  = 
3  COST' 
)tl=l,3 
f'iP  COST 
9),  1  =  1, 
NST  COS 
,9)  ,1  =  1 
M  CGST' 
1=1,3) 
DST«,3( 
9),  1  =  1, 
COST' ,3 
,9), 1=1 
ST'  ,3(/ 
,1=1,3) 
,3(/'  • 
9  )  ,  I  =  1 , 
MST  COS 
,9), 1=1 
M  CGST' 
,1=1,3) 
,3( /'  ' 
9  )  ,  I  =  1  , 

COST' , 
,9), 1=1 

COST', 


,  3(/'  ♦,->'F14.2)) 

1,3)  ,K=1,2) 

,6{/'  •,9F14.2)) 

) 

«,3(/'  ',9F14o2)) 

3) 

T  •  ,  3 ( /  '  •  ,  9F  1 4 . 2  )  ) 

,3) 

, 3(/'  •,9F14.2)) 

/•  ',9F14.2)) 

3) 

(/'  ',9ri4.0)) 

,3) 

'  •,9F14.0)) 

, 9F14«0)  ) 

3) 

T',3(/'  ',9^14.2)) 

,3) 

,3(/'  »,9F14.2)) 

,  9  F  1  4  .  2  )  ) 

3) 

3(/'   ',9F14.0)) 

,3) 

J(  /•  '  ,9F1^.0) ) 
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5  3 

5  5 
56 


32 


33 


34 


WRIT 
rCRM 
WRIT 
FORM 
DO  5 
DG  5 
TAVJT 
WRIT 
FORM 
WRIT 
CALL 
FORiV 

DC  3 
CALL 
LI  =  I 
DC  3 
CALL 
WRIT 
CALL 
CALL 
DO  3 
CALL 
LI  =  I 
DG  3 
CALL 
LJ=I 
DC  3 
CALL 
WRIT 
CALL 
CALL 
DO  3 
CALL 
LI-I 
DC  3 
CALL 
LJ  =  I 
DC  3 
CALL 
LK-I 
00  3 
CALL 
WRIT 
CALL 
CALL 

no  3 

CALL 
LI  =  I 
00  3 
CALL 
LJ=I 


E(6j  5 
AT(  « 
£{6,5 
AT(  • 
5  I--1 
5  J-1 
(  ItJ) 
E(6,5 
AT(  » 
E{6,3 

ZERO 
AT(  • 
2  I-l 

COST 

2  II- 
COST 

£(6,5 
GRAV 
ZERC 

3  1=1 
COST 

+  1 

3  Ii  = 
CCST 
I+l 

3  IJ  = 
COST 

E  (  6  ,  5 
GRAV 
ZERC 

4  I-l 
CCST 

+  1 

'->  11  = 

CCST 
I  +  l 
4  IJ  = 

CCST 
JH  1 

4  IK  = 
COST 

E(6,5 
OR  A '7 
ZERl 

5  1  =  1 
CCST 

+  1 

5  II- 
CCS1 
HI 


i)  { 

'/// 

4)  ( 

«/// 

,3 

,9 

=  AV,T 

6)  I 

«//■/ 

) 


(  A..jT(  I,  J),J  =  I,9)  ,  1  =  1,3) 

•  PRESENI  AV.T  CCNST  C0ST«,3(/'  'j9F14e2)) 

( AWTCM(  I  , J)  , J  =  l,g),  1  =  1,3) 

'     PRESENT    AV,T    0  +  M    COSTS  3(/'     SvF14o?)) 


(  I,  J  )  +  AV.TON{  I,  J  ) 

(     TAVyT{  I,J)  ,J  =  1,9)  ,1=1,3) 

'  PRESENT  A'rtT  C0ST«,3(/»  ',9F14.2)) 


',///'  CC*-Pi:JATIGN«  ,4X,9{  I6,4X)  ) 

(1,1) 

LI  t9 
(11,2) 
)  I,  II 
TY{2) 

,7 
(1,1) 

LI,R 
(11,2) 

LJ,9 
(IJ,3) 

)  }  ,  I  I , IJ 
TY(  3) 

,6 
(1,1  ) 

LIt7 
(11,2) 

LJtS 
(  IJ,3) 

LK,9 

(  IK,4) 

)  I ,  II  ,  IJ  ,  I  .< 

TY(4) 

,5 
(1,1) 

L  I  V  ': 

(  n,  2) 
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36 


DC  3^;  IJ  =  LJ,7 

CALL  CCST( IJ, 3) 

LK=IJ+1 

DO  35  IK-LK,8 

CALL  CCST(iK,4) 

LL=IK+1 

DO  3d  IL=LLr9 

CALL  C0ST(IL,5) 

WRITE(6t5)  I,II,IJjIKtIL 

CALL  GRAVTY(5) 

CALL  ZERC 

DO  36  I=lt4 

CALL  CCSTC 1,1) 

LI=I+l 

DC  36  II^LI;5 

CALL  COST( 11,2 ) 

LJ=II+1 

00  36  IJ=LJ56 

CALL  CCST(IJ,3) 

LK=IJ+1 

DC  36  IK=LK,7 

CALL  CCST{ IK, 4) 

LL=IK+l 

DO  36  IL=LLvS 

CALL  C0ST(IL,5) 

LM=IL+1 

DO  3  6  If—L.",-:* 

CALL  CCST  (  IM,6) 

vm  I  T  E  {  6  ,  b  )  I ,  I  I ,  I  J  ,  I  K  ,  I  L  »  I  ,V; 

CALL  GRAVTY(6) 

CALL  ZERC 

DO  37  I=L,3 

CALL  CCST( 1,1) 

LI=I+1 

DO  3  7  II=LI,4 

CALL  CCST( 11,2) 

LJ=I  I  +  l 

uL  37  1J=LJ,5 

CALL  CCSK  IJ,3) 

LK  =  IJ-H 

DC  37  IK=LK,6 

CALL  CCST(  IK,'t) 

LL=IK+1 

DC  37  IL=LL,7 

CALL  CCST ( IL,5) 

LM^lL+1 

DC  37  iy  =  Lf',8 

CALL  COST  (  1M,6) 

Lt>i=Ii'  +  l 

00  3  7  lrs-LN,9 

CALL  CCST  ( I  i,l) 
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v,SlTii(6,-3)  It  I  I,  IJt  [y^y  IL,IM,  IN 
CALL  GRAVTYi?) 

37  CALL  ZEaC 
OD  38  1=1,2 
CALL  CGST(i,l) 
LI=I+1 

DO  38  I  I  =  LI»3 

CALL  CCST( 11,2) 

LJ=II+1 

DC  38  IJ=LJ,4 

CALL  CCST(IJ,3) 

LK=IJ+1 

DO  3S  IK=LK,5 

CALL  CCST{IK,4) 

LL=IK+I 

DC  38  IL=LL,6 

CALL  C0Si(ILt5) 

LM=IL+1 

DC  38  IF^LN-T? 

CALL  COST  {  i:;,6) 

LN^liX  +  1 

DO  38  IN^L.M,8 

CALL  COST  (IN, 7) 

LG=IN+1 

DO  38  IC=LG,9 

CALL  COST  (10,8) 

WRITE  (6,5)  U  I  I,  IJt  IK,  I L,  If-- ,  IN,  IC 

CALL  GKAVTYO) 

38  CALL  ZERC 
DC  39  1=1,1 
CALL  CGSKI,!) 
LI=I+1 

DO  39  I1=LI,2 

CALL  C0ST(II,2) 

LJ=II+1 

DO  39  IJ=LJ,3 

CALL  CGST(IJ,3) 

LK=IJ+1 

DO  39  IK=LKt4 

CALL  CCSTdKi-^t) 

LL=IK+1 

DC  3  9  IL=LL,5 

CALL  CG5T  (  U.,5  ) 

LK=IL+1 

DC  '3J    1M=LK,6 

CALL  COST  {If'i,6) 

LN=If'  +  l 

00  39  Ivi  =  L>N,7 

CALL  COST  (1^,7) 

LG=IN+1 

DO  39  IC=LC,3 
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39 

30 


CALL 
LP=I 
DO  3 
CALL 

WRIT 
CALL 
CALL 
CONT 
STOP 
END 


COST  (10,8) 
C  +  l 
9  I  P  ^  L  P  ,  -^ 

COST  (IP, 9) 
E(6,5)  I,Il,IJ,IK,IL 

GRAVTY(9) 

ZERO 
IKUE 


IF»  IK,IG,  IP 


12 


SUBRGUTI 

COKHGN  F 

COf^MON  C 

*TRC{3»9) 

*  TC(3,9) 

*  PMP(3»9 

*  CCl(2t3 

*  PCCN'KZ 

*  TAWT(3, 
«  TRC2{3t 

*  2,3,9), 
COKHON  ^; 
DIFEKSIC 
DO  101  H 
DO  101  K 
SX=0.0 
SY=0.0 

1-1  =  0.0 
DC  2  J 
SX=SX+WT 
SY=SY+WT 
W=W+kTT ( 
X  =  SX/l-, 
Y=SY/W 
LSD  =  0 
F-  =  0.0 
0=0.0 
H  =  0.0 
WRITE (6, 
FORMAT ( • 
DO  1  1=1 
0D( I  )  =  (  ( 
WX( I )=WT 
WY(  n=WT 
F  =  F  +  V-iX(  i 
G=G+WY( I 
H  =  H  +  .iTr  ( 
X1=F/H 
Y1=G/H 


NE  ORAVTY(N) 


(2,9 
C(3, 
TRT  ( 
CX{9 
0)  ,C 
»9)  » 
jTRT 
,PMP 
TOM( 
(3,9 
DD{ 


),WT 
9)  ,R 
3)  rK 
)tCY 
OMU 
ACAP 
2(2, 
2(3, 
3,10 
),WT 
9 )  ,  K 


(3,9 

TRT( 

CP(3 

(9), 

2,3, 

(3,9 

3,10 

9)  ,R 

) 

T2(3 

X{9) 


) jCAP{ 3,9),CI(2,  3,9) , POP (3, 9) , 
3) ,RCP (3),KTC( 3),RPP(3) ,RTR(3)  , 
),RTC(3) ,RPP(3),RTR(3), 
PP(3,9) ,WTT(3,9),DIS(9),CP{3,9) , 
10),TCCl(2O»9)  ,TCCm(  2,3,9)  , 
),ACP{3,9) ,RACP(3)   , AWT (3, 10)  , 
),TRT3(2,3,10),FX2(2,3),FX(2,3), 
PP2(3)  f PP2(359)  ,FX3(2,3)  ,PCC1( 

,9) ,RTC213  ),TC2{3,9) 
,WY(9) 


1,>\ 
T(K, J)*CX( J) 
T(K,J  )';CY(  J) 
i<  r  J  ) 


12)  X,Y 

•  ,10X,2{ F10.6,5X)  ) 

,^ 

CX(  1  )-X)>.--.=  2+(CY(  I  )-Y  )**2  ) 

T(K,I )*CX{ I ) 

T(K,1 )*CY{  I) 

)/CD(  I  ) 

)/CD(  I  ) 

K , I) / C  C (  I  ) 


■r-r.5 


129 


5 
11 

20 

30 


31 


lOC- 


71 


72 


73 


70 


IF(AB5(X 
IF  (AeS( 
GG  TO  30 
LSD=LSD+ 
X  =  X1 

IF  (LSD. 
GO  TO  6 
WRITt(6, 
FORMAT { » 
FORMAT ( 

■•  ,2{H  Y  = 
00  31  L= 
RTRT{L)= 
RACP(L)= 
RCP{L)=C 
RTC(L)-0 
RTC2{L)= 
RPP(L)=0 
RPP2(L)= 
RTR(L)=G 
DO  100  J 
DIS( J)=S 
TC(K,  J)  = 
TC2{K,J) 
RACP(K)= 
RCP{K)=R 
RTC(K)=R 
RTC2(K>= 
RPP(K)=P 
RPP2{K)= 
CONTIKUE 
CALL  TRT 
LK=0.0 
RTR(K)=T 
IF(RTR(K 
WRITE{6r 
FORr^AK  » 

:^  •^X,'  M  = 
GO  TC  4 
RTR{K)=T 
IF(RTR{K 
WRITt{6, 
FCR/^AT(  ' 

^  ''t  X  7  •  M  = 
GO  TC  4 
LK=LK+1 
GO  TO  72 
v!RITt(6, 
WRIT  1:16, 
FORN-.AT{  ' 


l-X).GT.n.COCl)  GG  TO  3 
Yl-Y)  .GT.  COCCI)  GO  TO  3 


;T.50)  GO  TC  5 


11) 


•  tFl 
1,3 
0»0 
G.O 
.0 
.0 
0.0 
,0 
CO 
.0 

=  1 ,  H 
QRT(  { 
CIS(J 
=  CIS( 
ACP{K 
CP{K) 
TC(K) 
RTC2( 
P  {  K  ,  J 
RPP2( 


GREATER  THAN  50  ITERATIONS') 
•  CENTER  OF  GRAVITY'     ,20X,'  X=  •,F10.4 
0.4)  ) 


CX{ J )-X)**2+ 
)*WTT(K, J) 
J)*WTT2{K,J  ) 
, J)+RACP{K) 
+  CP(K,  J  ) 
+  TC(I<,J  ) 
K)+TC2( K 
)+RPP{K) 
K)+PP2(K 


{CY( J)-Y)-*2) 


J) 


J) 


(M,K,RCP(  K)  ,  -!aCP{K),  10) 

RT2(;''UK,  10)+RTC2(K)+RPP2(K) 
)  .GT.FX2(M,K)  )  GO  TC  70 

71)  F  X  2  (  M  ,  K  )  ,  -v  ,  l< 
'///'    TRT 2  INCiVICUAL  SUM'/'  ',F16.0, 
'  ,  12,  •       K=  • ,  12) 

RT3(M,K, 10)+RTC(K)>RPP(K ) 
)  .GT.FX3(  f-',K)  )  GO  TO  101 
73)  FX3(N-,K)  ,M,K 

«///«   TRT3  INDIVIDUAL  SU:^i«/»  •,F16,0, 

♦,12,'       K  =  ',12) 


20)  X,Y 

2  '3  )  (  T  C  (  K  ,  I)  ,  I  -  1  ,  N  ) 

ICFAL  TR^'  CCST'  /'  •  ,9ri4 


2) 


130 


26 


22 


27 


21 


23 


101 


FOR 
ING 
*CAP 
k'RI 
FOR 

*  • 
CC  = 
AA  = 
ViRl 

*  TOM 
FOR 

*  AA 
RTR 

ir-lRI 

FOR 
*ST, 

*  /t 

kRI 

FOR 

LK= 

IF( 

GO 

CON 

RET 

tND 


T£{6, 

MAT( 

COST 
ACITY 
TE{6t 
MAT{  ' 
•  t9Fl 
CCKH 
AWT{K 
TE(6, 
{K,10 
MAT(  • 
,AV.T, 
T(K)  = 
TE{6» 
MAT  ( 
TYPE 

S3{ 
TE{6, 
.VAT( 
LK  +  1 
LK.EG 
TC  72 
TINUt 
URN 


26)  (CX(  JJ)  tCY{  JJ)  yV,rT(K,JJ)  ,PP{K,  JJ)7CP(K,JJ  ) 

•  CCROIWATE  (X,Y)  , TRANSMISSION  COST^PUMP 
,  CAPACITY«/9(  «  •,5FI6.2)  ) 
'/^{/*     •,5F16.2)) 

22)  (  CIS  (  I),  1  =  1,:^) 

CISTANCE  FROM  SOURCE  TO  CENTER'/ 
4,4) 

,K,10)+CGf';l(f-'vK.  10) 
, 10)+AWT0M{K,10) 

27  )  CC  ,  CCl  {  M,  K  ,  10  )  ,  CO.Vl  {  M,  K  ,  1  0  )  ,  AA  ,  AWT  (  K?  1  0  )  ,  Aw 
) 

CONVENT  I  CNAL  TRT — CC  ,  CC  1,  CO.'-'l   ADVANCED  TRT — 
AUTOM'/'  ' ,6(F16cO,4X)  ) 
AA  +  CC 
21)  RTC  {  K  )  ,  R  PP  (  K  )  ,  RTRT  (  K)  ,  r. 

•  REGIONAL  TRANS  COST, PUMPING  COST, TRT  CO 
r',/'  « ,3(F16.2»5X)  ,  •  M=  ',12) 
F16.2,5X),»  /■'=  ',12) 

23)  RTR{K) 

'  REGID^AL  TCTAL'/'  «,F16.2) 

.2)  GC  TO  iOl 


17 
16 


18 
1 


FUNC 

COKM 

CONiV 

*TRC{ 

*  TC( 

*  PHP 

*  CCl 

*  PCC 

*  TAW 

*  TRC 

*  2,3 
GC  T 
GC  T 
t<  =  l 
L=16 
GC  T 
GC  T 
K=l 
L=16 
GO  T 
K  =  6 

L  =  ll 


TICN  TRv(  CI  ,F,i\) 

CN  KM 

CN  CCR(2,9),WT(3,^),CAP( 3,9) ,CI( ?,  3 , 9 ) , POP( 3 , 9 ) , 

3,9),TC(3,9) ,RTRT(3),RCP(3),RrC{3),RPP(3)  ,RTR(3) , 

3,9) ,RTRT{3) , RCP { 3 ) , RTC ( 3  )  , RPP ( 3 ) , RTR { 3 )  , 

(3,9) ,CX( 9 ),CY( 9), PP( 3,9)  ,WTT( 3,9), DIS( 9) ,CP{ 3,9) , 

(2,3,10), com ( 2, 3, 10), TC 01(2,3,9), TCOMl (2,3,9), 

rn(2,3,9) ,ACAP(3,9) ,ACP(3,9)  ,RACP(3)   , AWT (3, 10), 

T(3,9),TRT2(2,3,10),TRT3(2,3,10),FX2(2,3),FX{2,3), 

2(3,9),PMP2{3,9),RPP2(3),PP2(3,9),FX3{2,3),PCC1( 

,9) ,ArtT0M(3, IC) 

0  (16,17,  18) , KM 

C  (1,2,2)  ,  I 


(1, 


,3)  ,1 


C  4 
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GC  TC 
K=L6 
L=l 
t\  =  8.5 
B  =  201 
CC  =  15 
CCM=4 
R  =  0.0 
Fl  = 
F2=  . 
TCOMl 
TCCK 
DO  10 
TCOKl 
i=  M  ,  I  , 
TRM  =  T 
RETUR 
END 


5*G**0. 463-5-2. 0436 

.69/154.38 

24.7*A**1.37949*B 

C.0*(144.9/122.8) 

5 

.12  04  44 

0747619 

(  N  f  I  ,  .*-'  )  =  C  .  0 

N  ,  I  t  F  )  -C  C*  (  1  +  F  1  )  **  {  K+  1 )  /  (  1  +  rs )  «*K*L  /  50 

J  =  1,L 

(  N,  I  ,  K  )  -COM*  (  1+  F2  )  **  {  J  +  K  )   /  (  l-i-fl  )  **  {  J  vK-1  )  HCOiMl  ( 

}Ji) 

CCKKj  I»i'-  )+TCC;"l(N,  I,f;) 

N 


17 

16 


18 
1 


SUBROUTINE  TRT {K ,  I  , Q ♦ AQ , J ) 

COKf-^CN  Kf^ 

COi-'MGN  C0R(2,9),WT(3,9),CAF'(3t9),CI(2,  3  »  9)  ,  POP  (  3  ,  9  )  , 
^TRC{3,9) ,TC{3t9),RTRT{3) , RCP ( 3  )  , KTC ( 3 ) f R^P 1 3 ) , RTR ( 3 )  , 
^  TC(3,9) ,KTRT(3) ,RCPl3>,RTC(3)t«PP(3),RTR{3) , 
:^  PKP  (3,9),CX{9)tCY{9)vPP{3,9),WTT(3,9),DIS(9),CP(3,9), 
:=  CC1{2,3,  10)tC0Ml(  2,3,  10  )  ,  TCC  1  (  2  ,  3,9  )  ,  TCOMl  (  2,  3,9  )  , 
•~    PC0m(2,3,9)  ,aCaP{3,9),  ACPI  3,9)  ,RACr'{3)   ,  AWT  (3,  10), 
^    TAWT(3,9},TRT2(2,3,10),TRT3(2,3,10),FX2(2,3),FX(2,3), 
i^  TRC2(3,9),PI'.P2(3,9),RPP2(3)vPP2{3»9),FX3(2,3),PCCl( 
^  2,3,9)  ,Air.TnM(3,10) 

A=184. 56/129. S4 

8=144.9/110.8 

T=184. 56/173. 38 

S=144.9/1C2.5 

{  1  6  ,  1 7  ,  1 8  )  ,  ^■  M 
(1,2,2) 


GC  TG 
GC  TG 

L  =  16 

GC  TO 
GC  TG 

N!=l 

L  =  16 
GC  TC 
M  =  6 
L=ll 
GG  TC 
K=16 
L=l 
GC  TC 


4 
(  1 


I 


)  ,1 


(5,6) ,K 


CC  =  EXP(13.655  4)*Q**(  «62b39)-^=T 
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CGf'i= 
GC  T 
CC  =  E 
COM  = 
R=0. 
TSC  = 
TSCG 
RC  =  n 
RCO= 
AS=E 
ASG= 
AW  =  T 
AWTO 
Fl  = 
F2  = 
AViTO 
CO.Ml 
AWT  ( 
CC1{ 
DO  1 
AV-;TG 
^  -1) 

co^:l 
^  v-i 

TRT2 
TRT3 
KETU 
END 


ICO* 
0  7 
XP(13 
10.0* 
05 

EXF{- 
=  EXP( 
XP{-1 
cXPO 
XP(-0 
tXP{3 
SC+RC 
=  TSCO 
.10^ 
.0/4 
K  {  I  ,  J 
(Ktl, 
ItJ)  = 
K,I,  J 
0  N=l 

M  n ,  J 


(4.5046  +  0.6A96*ALOG10li^)  )*S 


.4586)*Q^ 
*{'.-.6720- 


-( .6456)*T 

0.6C23=.MLOG10{ 


)  ) 


05283)*AC**.7  761C2*10.*=^5*A 

14078) *AQ**. 69292  8*365*8 

2467)*AQ**.508C07*10.**6*A 

994  )*AQ**. 728564*365* B 

12  259)*AQ**1.1C74  4*10.**6*A 

12  21) *AQ**. 77  8987*365*8 

S 

CO+ASO 

4 

19 

0.0 

=  0.0 

*(  l  +  Fl)**{iV  +  i)     /(1+R)**M    *L/25 

CC*(1  +  F1)  **(Mfl  )     /  (  l  +  R)**iM*L/25 


AWTONK  I  ,  J  )  ^-AWTa*(  lvF2  )**  (  J+M       )  /  (  1  +  R)  **(N  +  iM 

(  K  ,  I  ,  J  )  -CCHl  (  K  ,  I  ,  J  )  +CCi;*  {  1  +  F2  )  **  {  N  +  N       )  /  il^R)  **  {  i\  + 

) 

(,<,  1  ,  J)=CC1{K,  I  ,  D+COi-'KKj  I,  J) 

(K,  I  ,  J)==TkT2(K,  I  ,  J)+A'-.T{  I  ,  J)-:A'/;TC;'(  I,  J) 

RN 


17 
16 


FUNG 

CONM 
*TRC( 

*  TC( 

*  Pf-iP 

*  CCl 

*  PCO 

*  TAh 

*  TRC 

*  2,3 
GO  T 
GC  T 
K  =  l 

L  =  16 
GC  T 
GC  T 
K=l 
L  =  16 
GC  T 
K  =  6 


TICN    PHG{G,I,N^,N) 

CN    Nf^ 

GN    CCR(2,9),WT(3,g),CAP(3,9) ,CI (2,    3 , 9 ) » POP ( 3 , 9) , 

3,9) ,TC{3,9),RTRT(3),RCP(3),RTC{3) .R^P ( 3 ) , RTR ( 3 ) , 

3,9) ,RTRT(3) ,RCP(3) , RTC ( 3  )  , RPP ( 3  )  ,RTR { 3 ) , 

(3,9) ,CX(9),CY(g),PP(3,9),wTT(3,9),DIS(9),CP(3,9) 

( 2 , 3 , 10 ) , CGM  U  2 , 3 , 10 ) , TCCl ( 2 , 3 , 9 ) , TCOMl { 2 , 3 , 9 )  , 

f'-l  (2,3,9)  ,  AC  AP(  3,9)  ,  ACP{  3,9)  ,RACP(3)        ,AWT{3,10)  , 

T(3,9),TRT2{2,3, 10 ) , TRT3 ( 2, 3 , 10 ) , FX2 ( 2 , 3 ) , FX ( 2 , 3 ) 

2 ( 3 , 9 ) , PMP2{ 3 , 9 ) , RPP2 ( 3 ) , PP2 { 3 , 9  )  , FX  3 { 2 , 3  )  , PCC 1 ( 

,9) ,AWTGH(3,10) 

0    {16,17,18) ,KV 

C     (  1  ,  2  ,  2  )  ,  1 


C    4 

C     (  1,2,3)  ,  I 
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L  =  ll 
GO    TG    A 
K  =  16 

L  =  l 

CC  =  76  300.  0=:^C*-(.  7682  )*(  201.69/15^1.  38) 

COM- 1 .6 10  3*Q** ( . 7  37A  3 ) *  10 . 0* { 144  o9/ 1 22. 8 ) *365 

R=C.05 

Fl=    .104214 

F2=    .0747619 

PCONUN,  I,f-')-0.0 

PCCl  {  N,  I  ,  K  )-CC*(  l-^F  U  ^*  ( i<  +  1)     /  (  1  +  R  )*"K':-L  /25 

00    10    J=1tL 

PCO^:i{N,  l,l')=    PCG?n{'^,  I,M)+CCM^={l  +  F2)**(  Ji-!s 

5={  J  +  K-1) 

Pi^iG^-PCCl  (N,  I, 

RETURN 

tKO 


)/( 1+R)** 


)4-pcg;-;i(n,  i,m) 


SUBROUTINE    ZERO 

COr'MCN  VH  ^    ^^^^ 

COPr'CN    CCR(2,9),,iT(3,g),CAP(3,9),Cn2,     3, 9), POP 

*TRC(3,9) tTC(3,9),RTRTI3),RCP{3),HTC(3),RPP(3)tR 

*  TC(3»9),RTRT(3),KCP(3),RTC{ 3),RPP(3),RTR(3), 

*  PMP  {  3,9  }  »CX(  9  )  ,CY  (  9  )  ,  PP  (  3,9  )  vVlTT  (  3,  9  )  ,  D  IS(  9  )  ,C 

*  CC  I  ( 2 , 3 , 10 ) , coy  1 ( 2 , 3 , 10 ) , TCC 1 ( 2 , 3, 9 ) , TCOMl ( 2 , 3 

*  PC0Fa(2,3,9) ,ACAP{3,9),ACP(3,9),RACP(3)   ,AWT( 

*  T  AWT ( 3 , 9 ) ,TRT2( 2, 3 , 10 ) , TRT3 ( 2, 3 , 10 )  , FX2 ( 2, 3 ) , F 

*  TRC2(3,9),pr--P2{3,9)  ,  RPP2  (  3  )  ,  PP2  {  3  ,  9  )  ,  FX3  (  2  ,  3  )  , 

*  2,3,9),A/.TGr:  (3,10) 
DG  1  1=1,2 

DO  1  J=l,3 

FX2{ I,J)=0.0 
FX3{ I,J)=C.0 
RETURN 
END 


13,9), 
TR(3) , 

P(3,9) , 
,9), 
3,10) , 
X  (  2  ,  3  )  , 
PCCU 


SUBROUTINE  CaST(  I  ,  J) 

C0Mf-:CN  ^^' 

CO^'FG^  CGR{2,9),WT(3,9),CAP(3,9),CI(2,  3  ,  9  )  ,  PGP  (  3  ,  9  )  , 
*TRC13,9) ,TC(3,9) , RTRT { 3 ) , RCP ( 3 ) , RTC ( 3 ) , RPP ( 3 ) , KTR ( 3 ) , 

*  TC(3,9),RTRT(3),KCP(3),RTC(3),RPP(3),RTR(3), 

*  PKP(3,9),CX(9),CV{9),PP13,9),WTT{3,9),DIS{9)  ,CP(3,9)  , 

*  C  C 1  {  2  ,  3  ,  1 0  )  ,  C  0  M  U  2  ,  3  ,  1 0  )  ,  T  C  C 1  {  2  ,  3  ,  9  )  ,  T  C  0  M 1  (  2  ,  3  ,  9  )  , 

*  PCGK1{2,3,9) , ACAP( 3,9) ,ACP( 3,9) ,RACH{3)   ,AWT(3,10), 

*  TAV.T(3,9),TRT2{2,3,10),TRT3(2,3,10)  ,  FX2  (  2  ,  3  )  ,  FX  (  2  ,  3  )  , 

*  TRC2(3,9),Pr-p2(3,9)  ,  RPP2  (  3  )  ,  PP2  (  3  ,  9  )  ,  FX3  (  2  ,  3  )  ,PCL  1  ( 

*  2,3,9)  ,AUTC."-'(3,10) 

CG(^^'CN  »^T2(3,9),WTT?(  3,9),RTC2{3  ),TC2(3,9) 

DG  10  K=l,3 
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/* 


CX(  J)=CC;M  1,1  ) 

V.TT{K,J)=wT(Kt  I) 

ViTT2{K,J)=WT2{Kt  I  ) 

ACPCK,  J)=/iCAP{K,  I  ) 

PP2(K,J)=PMP2(K, I ) 

CP(K, J)=CAP(K, I) 

PP(K»J)=PFP{K, I) 

CYC J)=CCR{2,I ) 

DO  10  M=l,2 

FX2(M,K)=rX2(  H,K)-i-TRT2  (f!,K,  I  ) 

FX3 ( M 1 K )  =  FX3 ( H ,K ) +TKT3 ( M , K,  I  ) 

RETURN 

END 
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